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MOTION OF CONDUCTING BODIES IN A MAGNETIC FIELD 


This book represents a collection of articles on the 
theory of induction, magnetohydrodynamic (MHD) machines 
with a liquid metal working medium, the electromagnetic 
processes in an ideal, conducting MHD machine, and the 
higher spatial harmonics of the magnetic field of an 
induction MHD machine. It also discusses the transverse 
edge effect in plane, induction MHD machines, the longi- 
tudinal edge effect in linear MHD machines, the pondero- 
motive forces influencing conducting media in the traveling 
magnetic field of a cylindrical inductor, and a theory for 
the propagation of pulsed electromagnetic fields in moving, 
conductive media. 


From the Editorial Board [5* 


The book Dvizheniye provodyashchikh tel v_magnitnom pole (Motion of Con- 
ductive Bodies in a Magnetic Field) which we offer to the reader is a survey of 
work on electromagnetic processes in magnetohydrodynamic (MHD) induction ma- 
chines (generators, pumps) with working media of liquid metal. 


The topics in this collection touch on the theory of MHD machines examined 
from the so-called electrodynamic approach, i.e., disregarding the magnetohy- 
drodynamic effects in the machine channel. 


The surveys include most of the papers published on the problem in ques- 
tion and several which were in the process of publication as the manuscript 
was being prepared for press; some original findings are also presented. 


It is proposed in the immediate future to devote a separate work to a 
survey of magnetohydrodynamic phenomena in MHD induction machines. 


Please send all requests and criticisms about the book to the Institute of 
Physics of the Academy of Sciences, Latvian SSR at 19 Turgenev St., Riga. 


STATE OF THE THEORY OF MAGNETOHYDRODYNAMIC INDUCTION MACHINES 
WITH WORKING MEDIA OF LIQUID METAL 


Ya. Ya.” Liyelpeter 


Magnetohydrodynamic (MHD) induction machines with a working media of 
liquid metal may be used as converters of the mechanical energy in the liquid 
metal flow into electrical energy (generators) or, on the contrary, as con- 
verters of electrical into mechanical energy (motors, pumps, brakes). Up till 
now liquid-metal MHD machines have found their chief use as pumps for trans- 
ferring liquid metal. The use of these machines for generating electrical 
energy involves the development of methods for converting thermal energy into 


* Ninbade in tie margin indicate pagination in the original foreign text. 


energy of liquid metal flow. Elliot's proposal of a two-phase conversion 

cycle is noteworthy (Ref. 1); other systems are also well known. Comparatively 
few papers (Refs. 2, 3) have been published on liquid-metal MHD generators in 
contrast to plasma generators. 


Approximate numerical evaluations of the basic characteristics of MHD 
induction generators indicate that at powers of the order of thousands of 
kilowatts their efficiency reaches 50-702. 


Papers by A. I. Vol'dek (Ref. 4) and N. M. Okhremenko (Ref. 5) give a 
synopsis of the theory of induction pumps. Below we shall scrutinize several 
aspects of the state of the theory, which are explained in less detail in the 
articles mentioned. 


A ianndy beter rictaral diagrams of MHD induction machines are known. The 
design of a specific MHD machine must always solve two problems: (1) selection 
of the optimum structural plan, and (2) finding the optimum relationships be- 
tween the design variables for the given system. In optimizing the design of 
a specific machine, we must know the numerical relationships between the phys- 
ical quantities determining its properties. In other words, we must analyze 
the characteristics of the specific machine with prescribed specific loads and 
geometrical dimensions. 


From this it follows that the basic tasks in MHD induction machine 
theory are to analyze the properties of a set of structural diagrams at pre- /6 
set specific loads and geometrical dimensions and to elaborate methods of 
optimizing dimensions and specific loads at the prescribed useful power. The 
first part of the task is by now in a comparatively well-developed state. 
The author gives a survey of the articles devoted to matters of optimization 
(Ref. 6). In the following we will touch only upon the first part of the 
theory of MHD induction machines. 


The design of an MHD machine with liquid-metal working medium must 
unite the properties of electrical and hydraulic machines, and to the maximum 
degree must satisfy the specific demands made on both types of machine. The 
result of the latter is obviously that the channel of the optimum machine will 
have a cross-section varying over its length and a variable mean flow velocity 
even in the machine's active zone, while the electromagnetic and hydraulic 
processes in the channel will be complexly interconnected, making it impossible 
to study them separately. In this form the theory is at present completely 
undeveloped. In all the problems investigated, the channel cross-section 
is assumed to be constant along its length. 


A concurrent solution of the electromagnetic-field and hydrodynamic 
equations in the last case is exceptionally difficult. Therefore, in practice, 
recourse is had to various sorts of approximations, e.g., the electromagnetic 
processes and the hydraulic phenomena are examined separately. Therefore, we 
may conventionally speak of the electromagnetic theory and the hydraulic theory 
of MHD machines, although they are necessarily interconnected. The first por- 
tion has, at present, been developed to the greatest extent. The second is 
for the time being in the initial stage of elaboration. 


In electromagnetic theory, a liquid metal is regarded as a solid body, 
e.g., a strip, cylinder, or the like, which moves at a velocity equal to the 
average speed of the liquid. In this case, we must determine the electro- 
dynamic force field, integral force, and energy (power) transmitted from the 
winding to the liquid metal (in the motor and brake regime), or in the other 
direction (in the generator regime). From these data we may compile the 
equivalent electrical circuit for the machine and its parameters, and may 
analyze its external characteristics and operating conditions. 


Consideration of hydrodynamic effects in the simplest case is reduced 
to the fact that part of the electromagnetic pressure is expended on equa- 
lizing internal frictional forces in the liquid. Hence, the employed pressure 
is less than the electromagnetic pressure. Hydraulic losses in friction in 
this case are determined for the prescribed channel shape, due consideration 
being given to the effect of the magnetic field on hydraulic resistance. 


In the approximation approach to MHD machine theory, one problem of /7 
hydraulic theory is to define certain corrections to electromagnetic theory. 
In principle, this route may be justified to some degree, for we may assume 
that in the optimum machine the distribution of electromagnetic forces over 
the channel cross-section should not be very uneven, Nor will the distribu- 
tion of averaged velocities substantially differ from the configuration during 
flow because of external pressure forces. The need for more detailed study 
of magnetohydrodynamic characteristics is, however, quite obvious. 


Let us examine in somewhat more detail the state of MHD induction machine 
electromagnetic theory, about which the largest number of papers has been 
published. This part of the theory comprises the study of electromagnetic 
field structure in the working gap of the machine, as well as of processes in 
the magnetic and electrical circuits which determine the effective resistance 
and reactance of the windings. Magnetic and electrical field structure in 
the working gap depends on the structural shapes of the stators (inductors) 
and the channel configuration. Every structural layout of the machine has its 
own specific phenomena which have a substantial effect on the machine character- 
istics, e.g., in a plane induction pump an important role is played by a trans- 
verse edge effect which is lacking in a cylindrical pump. 


Assuming the magnetic permeability of the magnetic circuit to be constant, 
we may describe the electromagnetic processes by a system of linear differen- 
tial equations. When studying field structure in the working gap, we usually 
assume that one of the components of magnetic induction on the inductor sur- 
face is given, and we shall relate this component with the current in the mach- 
ine windings. 


The plane linear MHD machine has been studied in greatest detail. Ref- 
erences 7-11 have discussed electromagnetic processes while neglecting the 
longitudinal and transverse edge effects, on the assumption that only the fun- 
damental harmonics of the magnetic field is present. In this case, the fields 
depend only on a single coordinate directed across the channel. Cases are ana- 
lyzed with conductive channel walls (Refs. 8, 9) and without them (Refs. 7, 10, 


11). References 16-22 study the effect of the higher spatial harmonics of 
the magnetic field, which are caused by discrete spacing of the winding con- 
ductors and by serration of the inductor surface. This subject has also been 
touched upon in (Ref. 23). The theory of MHD induction machines may likewise 
partially use detailed studies of the higher harmonics in ordinary electrical 
machines. A comprehensive survey of this subject is found in Geller and 
Gamata's work (Ref. 40). 


A comparatively large number of studies deal with the transverse edge 
effect in a plane MHD induction machine (Refs. 24-39). They solve prob- /8 
lems which are stated in various ways and which examine the motion of a 
conduction band of finite width in a traveling magnetic field. The basic find- 
ing of these studies may be reduced to computing the coefficient of pressure 
attenuation and analyzing the influence of the transverse edge effect on 
characteristics of the machine. A new principle for simulating vortex fields 
in the conduction band is elaborated in (Ref. 33). References 34, 35, and 38 
examine the transverse edge effect when there are shortcircuited buses in the 
channel of a plane MHD machine. 


Despite the great number of published works, the question of transverse 
edge effect cannot be assumed to be thoroughly studied since many solutions 
(Refs. 26-28, 30, 31) are obtained for various initial conditions, while no 
general analysis or sufficiently comprehensive comparison of results with 
experimentation has been conducted. 


A number of works have studied magnetic field structure of an inductor 
of finite length and the problems associated therewith of phase asymmetry in 
the winding (Refs. 41-46). The totality of these phenomena is usually called 
the longitudinal edge effect in the primary circuit. It has been ascertained 
that a break in the magnetic circuit leads to the generation of pulsed com- 
ponents of the magnetic field along with the traveling component. Methods 
have been proposed for balancing the pulsating fields. Distortion of field 
structure associated with motion of the conducting medium through the final 
zone of propagation of the traveling magnetic field (longitudinal edge effect 
in the secondary circuit) has been partially studied (Ref. 47). (Ref. 48) 
examines the edge effect in the primary and in the secondary circuit together. 


A survey of longitudinal edge effect is found in (ref. 4). 


The theory of cylindrical MHD induction machines is treated in fewer 
works than is that of plane machines, although this theory is more promising 
from the electromagnetic viewpoint (there is no transverse edge effect) since 
in practical computations the effect of channel curvature on field distribution 
may in the vast majority of cases be neglected, i.e., the problem may be re- 
garded as a plane one. (Refs. 7, 49-54) scrutinize various problems for cy- 
lindrical MHD machines with constant channel cross-section without regard to 
longitudinal edge effect. (Ref. 49) computes the ponderomotive force acting 
on a conductive cylinder of finite length in an infinite cylindrical inductor 
and confirms the fact (which was observed experimentally previously) that force 
density peaks at a certain relative cylinder length. 


The cited list of works on various problems of the motion of conducting 
bodies in a traveling magnetic field is not exhaustive. There are sev- /9 
eral experimental works studying magnetic field structure of an actual in- 
ductor in a MHD machine, the ponderomotive forces acting on solid conductive 
bodies, operating characteristics of the machine, and other effects (Refs. 55- 
59), but they are considerably fewer in number than the theoretical works. 


When we speak of the MHD machine electromagnetic theory on the whole, we 
may note that there is a substantial number of works on this theory, but that 
by no means all matters have been investigated with adequate care. While the 
solutions themselves of electromagnetic field equations have been obtained for 
a large number of different cases, there is usually inadequate analysis of 
these solutions and comparison of them with other similar problems. The find- 
ings of the studies are rarely reduced to convenient computational coefficients 
and formulas defining the parameters of equivalent electrical circuits or con- 
veying the basic energy characteristics of the machine - power, efficiency, 
power factor, etc. All this impedes practical utilization of research find- 
ings in the design of specific MHD machines. 


One purpose of the present collection is to make a critical survey and 
compare the findings of papers on the electromagnetic theory of MHD induction 
machines in order to facilitate their use in design work. The collection 
adduces surveys of individual problems in the theory, e.g., on longitudinal 
and transverse edge effect, the influence of higher spatial field harmonics, 
etc. A separate range of questions is considered in the studies dealing with 
motion of conducting bodies in pulsed electromagnetic fields. A survey of 
these works is given in connection with their possible use in metering tech- 
nology, e.g., in measuring the flow rate of a conducting liquid (Refs. 59-60). 
Various devices with traveling and pulsed magnetic fields are also applied 
for the same purpose (Refs. 61-63). 


In conclusion let us briefly touch upon the state of the hydraulic theory 
of MHD induction machines. 


Accurate solutions to problems of conductive liquid motion in a traveling 
magnetic field may be obtained only for the simplest cases of laminar flow, 
e.g., where velocity has only a single component, and magnetic field amplitude 
in the gap is either constant or varies in accord with a given law. Several 
such solutions have been published (Refs. 8, 64-69), but their results may be 
used merely to give a qualitative representation of the nature of the flow, 
since laminar flow conditions are very infrequently observed. 


Of definite practical interest are the results of an experimental study 
of turbulent flow of a conducting fluid in a constant transverse magnetic /10 
field (Refs. 70-73). These experiments have been used to derive empirical 
relationships for the coefficient of hydraulic resistance in smooth and rough 
channels of constant cross-section, while taking into account the influence of 
the magnetic field. To design MHD induction machines, however, these findings 
must be utilized while taking into consideration the fact that the drag co- 
efficient may be different in a traveling field. In this connection the results 


of N. M. Okhremenko's article (Ref. 74) on the effect of a traveling magnetic 
field on the drag coefficient are very interesting. 


Solutions are known for the distribution of electromagnetic fields in a 
conducting fluid moving with the velocity distribution given by a power law 
(Ref. 75). 


Finally, there are certain experimental studies on fluid flow structure 
in the flow-through parts (channels) of MHD machines disregarding the magnetic 
field effect (Ref. 76). The chief purpose of these studies is to define the 
hydraulically optimum channel shapes - in particular their intake and outlet 


zones. 
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ELECTROMAGNETIC PROCESSES IN AN IDEAL, INDUCTION MHD MACHINE 


A. K. Veze, L. Ya. Ulmanis 


1. Introduction 


The basic theory of induction MHD machines is concerned with the cal- 
culation of electromagnetic processes in an infinitely wide conductive band 
located in the traveling electromagnetic field of a flat inductor. If it is 
assumed that the conductive layer is infinitely long and that the sinusoidal 
traveling wave produced by the idealized inductor has only a basic harmonic 
component, such a problem may be readily solved analytically. The first prob- 
lems of this type were solved in conjunction with the development of a theory 
for asynchronous engines (Ref. 1-2). Subsequently, such problems were solved 
for MHD machines. The results obtained were utilized for designing induction 
pumps (Ref. 3-6), but the theoretical formulas had to be refined by the intro- 
duction of certain empirical coefficients. These coefficients were employed 
to provide an approximate determination of the difference between the assump- 
tion advanced in the theoretical formulation of the problem and actual practice 
(Ref. 8, 18, 26). 


The results obtained when such problems were solved were also applied 
when designing MHD generators, brakes, flow meters, electromagnetic mixers of 
metals in melting furnaces, chutes for transporting molten ferrous metals, and 
other equipment. 


This article presents a summary of the solutions for electromagnetic 
processes in a conductive band located in a traveling magnetic field of a flat 
inductor, under the assumption that the dimensions of the device are infinite 
in the direction of motion of the field and in the direction in which the elec- 
tric current passes. 


2. Formulation of the Problem 


A different number of layers with differing electroconductivity may be 
located in the operational zone of MHD machines. For example, an insulation 
layer, the channel wall, and in the middle a layer of electroconductive 
liquid are usually located on each side of the inductor in the clearance of 
flat induction pumps operating symmetrically, i.e., five layers in all. 
Pumps are also being designed for pumping two liquids at once; in such /16 
devices, the number of layers is greater. Therefore, it is advantageous 
to examine the problem in the general case, assuming the presence of many 
layers in the inductor. 


A flat ideal inductor, producing a traveling magnetic field, may be 
represented in the form of one or two infinitely thin layers of current, whose 
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linear density (linear current load) 
has only one tangential component. 


When there are two layers of 
current, we shall call the inductor a 
two-sided inductor, and when there is 
one layer we shall it a one-sided in- 
ductor. 


Our problem may be formulated as 
follows. 


There are n parallel layers (Fig- 


Figure 1 ure 1) in the traveling magnetic field 
of a flat inductor. The specific con- 
Arrangement of Conductive ductivity of these layers is Ojos Soo os 


Layers in a Two-sided Inductor 
s +> Oo» and the thickness is b> bos si 


Es ba? respectively. We shall assume that all the layers move at constant vel- 
ocities with respect to the inductor equalling Vy> Vgs tree Va In order to 


compute the velocity of each layer, let us introduce the following notation: 


U,—U. 
s)=— ; 
Vs 
Uz — U. 
S= 2 2 ; 
Us 
Cak 
$n = oF 3, 


where v, = 2tf is the velocity of the traveling magnetic field with respect to 


the inductor (t - polar division of the inductor, f- frequency of the current 
supplying the inductor). We may assume that the magnetic permeability of /17_ 
all the layers, “ioi the magnetic circuit of the inductor, is constant, and 
equals U= 4r *°10742. It equals u = Const. in the magnetic circuit of 


the inductor. The sone division of the inductor, and also velocity of the 
field along both of its sides are assumed to be the same. 


The magnetic field H the vector potential A and the electric field 
E have the following components in the given case: 


H (Hx; 0; Ha); 
A {0; Ay; 0}; 
E (0; Ex; 0}. 
From this point on, we shall designate Ay simply by A and E7 E. 
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The following differential equation may be solved for each of the k 
layers 
OPA, | Ar OA, 
ax? tz? Hoe og (1) 


In the case under consideration, all the electromagnetic quantities 
have a sinusoidal dependence on wst ~ az, and they may be expressed by means 
of functions having the following form 


x= Xo (x) ellast—az) | 


where x is the corresponding electromagnetic quantity. Equation (1) then 
assumes the following form 


a? 
Fok — Br 2o =0, (2) 


where 


Ba= J a?+ ipon@sSn= pir + ipar; 


(3) 
2 2 a 
yum ]/ anosa Funcosetor a /Yatti +1+1 S (4) 
CETT F UROROSk — a? oe see i- 1 — (5) 
rn = Qrpar; 
_ BROROSR (6) 
= ee, ae . 


a 


On surfaces dividing adjacent layers, the magnetic and electric fields 
must satisfy the following boundary conditions: 
j /18 
Byu=By, Hy=H» for x=b\— 55 
Bo=Ba; Ho= $; 
x2=Bx; Ha=H3 for x=bi+b:- 5; 
(7) 


Bar =B rtn; A, =H, zakt) for *=h4+.. tomt 


(where 6 is the distance between both sides of the inductor). 


We shall define one of the components of the traveling magnetic field on 
the inductor surface ~- i.e., in the case of x = +.. Practical computations 


have shown that it is more advantageous to define the tangential component of 
the field, since - if the magnetic permeability of the inductor magnetic circuit 
is u = @ n numerically it equals the linear current load of the inductor, and 
does not depend on processes occurring in the electroconductive layers. 


Thus, in the general case the boundary conditions on the inductor surface 
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may be formulated as follows: 


Hz = Hoye(ost-22) for x=, 

(8) 
H,= — Hoge!lost-az+9) for xed. 
Here ¢ is the angle at which the phase shifts between the magnetic fields on 
both sides of the inductor, and it assumes any value in the general case. In 
practice, » = 0 customarily holds - i.e., the windings on both sides of the 
inductor are switched onconcurrently. In the literature, with the exception of 
several articles (Ref. 9, 10), it is this case which is primarily investigated. 


In all, there are 2n boundary conditions which are requisite in order 
to determine 2n integration constants. If it is assumed that in the case of 
xe È the linear current density equals zero, we have a one-sided inductor. 


If the problem is symmetrical with respect to the x = 0 plane - i.e., 
Er Þa- (k-1)° o = Ta- (k-1) and ẹ¢ = 0 - the number of layers n is odd, since 


for even n there would have to be two adjacent layers with different physical 
parameters in the middle of the clearance. However, in this case the problem 
is non-symmetrical. The solution of the symmetrical problem in the mean band 
is expressed by the even or odd function of x depending on whether the /19 
desired quantity has even or odd symmetry. However, solutions in symmetrical 
bands, which do not come in contact with each other, are characterized by a 
sign before the coordinate x for quantities with even symmetry, and are char- 
acterized by a sign before the coordinate x and before every expression for 
quantities with odd symmetry. 


b 


There is no necessity of finding a solution for every layer when solving 
the symmetrical problem. We may find a solution for only one half of the 
symmetrical region, assuming that there is a tangential component which equals 
zero in the middle of the clearance. The solution for the other half of the 
symmetrical region may be found by using the relationships presented above 
between the solutions in symmetrical layers. It is simpler to solve the problem 
by employing symmetry, since we must determine only n + 1 constant, instead of 
2n constants. 


The force lines of the magnetic field may be determined by solving the 
differential equation 


dx = Refs 
dz ReH, (9) 


The force acting upon the conductive band is composed of a constant and 
pulsating components. The latter pulsates with the double frequency of the 
traveling magnetic field. The force may be computed analytically by multiplying 
the instantaneous values of current density and magnetic induction. Pulsations 
of the force density may have an influence upon the hydrodynamic processes in 
liquids located in a traveling magnetic field. When the ponderomotive forces 
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in induction pumps and other MHD machines are investigated, the pulsation of 
the force component is usually disregarded, and its constant component is com- 
puted, which equals the force averaged over time. 


In the problem under consideration, there are two components of the 
force density - the tangential component T: and the normal component fz Their 


constant components may be computed according to the following formulas 
1 j 
hay Re(j,8.*) (10) 
and 
1 
z= — 7 Re(/yBz*). m 


They depend only on the x-coordinate. The electromagnetic pressure per unit 
of length P3 in the k-th layer, which is produced by the traveling magnetic 


field - i.e., the density of the tangential force component, is /20 
xotd, 
pa=~oh- | Re(jsBs*)dx, (12) 
26, 
Xa 
and the mean density of the normal force component is 
Xotb, 
I ; 
pam aif Re(j,B,*) dx. (13) 
Xo 


In the symmetrical case, the tangential component of the force density 
has even symmetry, and the normal component has odd symmetry. This means that 
the normal component of the resulting: force equals zero in the symmetrical 
case. It only contracts the body, but does not displace it. 


The power transmitted to the conductive band from the inductor by means 
of the traveling magnetic field may be expressed by means of the normal com- 
ponent of the Poynting vector, which expresses the density of the electromag- 
netic energy flux per unit of time: 


S= FEV") (14) 
The active power transmitted to the body through a unit of surface area 


equals the real part of the the complex Sa and the reactive power equals the 


e 
imaginary part of the complex 5° 


If the coordinate system is related to the moving conductive band, then 
the Poynting vector may express that portion of the energy which is converted 
into Joule heat losses. If the coordinate system is related to the inductor, 
S expresses that portion of the electromagnetic energy which is converted into 


mechanical energy. 
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The studies (Ref. 1-25) solved different problems of the type examined 
above. Let us analyze certain cases, which are of the greatest interest. 


3. Electromagnetic Phenomena in a Traveling Field of a One-Sided 
Inductor 


(a) Conductive Half-Space 


In a study by I. M. Kirko (Ref. 11) the simplest problem of this type 
was investigated - the motion of conductive half-space in a traveling magnetic 


filed. 


A conductive medium with the specific electric conductivity o fills /21 
the half-space x > 0, and moves at a constant velocity v along the direction 
of motion of a traveling magnetic field (Figure 2), which is produced by an 


infinitely thin current layer located on the surface x = 0. In the case of 
x < 0, it is assumed that the half-space is filled by an ideal ferromagnet 


with up = œ and o = 0. 


X van The solution of Maxwell equations 
provides the following values for the 
components of the electromagnetic field 


a in conductive half-space: 
VIVVOTTTITTE Ot j e 
Gs Boum nes, a9 
Figure 2 Bim= woAce®*, (16) 
Half-Space Above a One~-Sided E,,= foposAo e- br, (17) 
Inductor B 


je iw toasAg e-bx, (18) 


Thus, all of the field components and the current density decrease 
according to an exponential law as one recedes from the surface x = 0. 


At the distance + from the surface x = 0, the magnetic field decreases 
1. 
by a factor of e - i.e., this quantity represents the depth to which the travel- 
ing magnetic field penetrates (Ref. 3) 


1 1 2 
h= = = |/=. (19) 
= Je 
In the special case, when e > 0 - i.e., when the conductive medium is 
absent or fs = 0 - we have 1 
ho=hr=—- 
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In this case, the pole division of the inductor t - the quantity 
which is contained in the dimensionless frequency e - determines the depth to 
which the traveling magnetic field penetrates. This phenomenon is arbitrarily 
designated as the geometric surface effect. 


The surface effect in a traveling magnetic field is characterized /22 
by the fact that the planes of the same phase of the field are not perpen- 
dicular to the x > 0 plane, but are inclined toward it at the following angle 


1 
P= tant. 
2 
We obtain the amplitude of the normal component of magnetic induction 
from expression (15) 


Vite (20) 


and the tangent of the angle at which the phase shifts x < 0, as compared with 
the phase surface density of the current Ae: 


Yı cos pox— po Sin pox 


tan Ọx F — à 
pi Sin p2x + wpe COS tpox 


(21) 


The ratio of the amplitudes of the tangential and normal induction com- 
ponents increases with an increase in €; 


4 
B.T l+e’. (22) 


The angle at which the phase of the tangential induction component shifts 
changes proportionally to the x coordinate; 


P= — peor. (23) 


All of the quantities investigated are determined by the parameter e. 


Force Lines of the Traveling Magnetic Field Above a One-Sided 
Inductor in the Case of e = 0 


Figure 3 


19 


Figure 4 


Force Lines of a Traveling Magnetic Field Above a One-Sided 
Inductor in the Case of te = 50 


Figures 3 and 4 illustrate the magnetic force lines in an infinite half- 
inductor for two values of €, which were calculated 


space above a one-sided 

by Yu K. Krumin. Similar problems were investigated by Schilder in (Ref. 23). 
The parameter € also determines the density of the pondermotive force 

averaged over time and the electromagnetic power; 


1 € 
=— 2 — — fix 
fz 2 apolo yita e , (24) 


1 e 
=- 2 -4x 
fx p) p2poAo Vite ent, (25) 
or f 123 
a= =p; 
n (26) 
1 w ‘pe 
Pa=x* g Hoo yI + 7o’ 
1 w » 
agg wae ae (27) 


The power factor of an idealized energy convertor, without allowance 


for the leakage flux of the winding, is 
de 
cos Q= =>. 
=g (28) 
The dimensionless quantity € is an important characteristic of the /24 


electromagnetic processes in problems of the type which we are considering. 
It was called electromagnetic slipping in a study by Kh. I. Yanes (Ref. 19). 
This quantity may be also regarded as the magnetic Reynolds number 
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Regents 


Yin 


in which the velocity u is represented by the relative velocity of motion of 
the traveling magnetic field and of the conductive medium as - The distance 


2 a is assumed to be a characteristic linear dimension, and the magnetic 
viscosity equals v_ = 1 
m Hod 


The quantity € changes sign when the slipping sign changes: in the 
braking and pumping regimes £ > 0, and in the generator regime of the MHD 
machine € < 0. 


(b) Conductive Layer in the Field of a One-Sided Inductor 


The articles (Ref. 12, 13) investigated the following problem concerning 
three layers: the conductive layer which moves at a constant velocity 2tf 
(1 - s) with respect to the inductor is located at the distance 6 from the one- 
sided inductor of a traveling magnetic field (thickness of the conductive lay- 
er, b; specific conductivity, o; specific conductivity of the regions I and 
III equals zero). The non-conductive layer behind the conductor is not limited 
(Figure 5). 


Let us analyze certain re- /25 
lationships obtained when the Max- 
well equations are solved for the given 
problem. 


In order to simplify the analysis, 
let us introduce the following dimen- 
sionless parameters, except for e: 


b=6a, (29) 
Fi 5 
igure 3=dc, (30) 
Solution of the Problem Regarding Seya (31) 
Three Layers in the Field of a iig (32) 
One-Sided Inductor 


@ = poowsb?. 


= The dependence of the induction components on the dimensionless coordinate 
x is shown in Figures 6 and 7. Bọ is the tangential induction component on 
the inductor surface, and ¢ is the angle at which the phase shifts between both 
components with respect to the phase of the tangential component on the inductor 
surface. The solid curves refer to the normal induction component; the broken 
curves refer to the tangential component. Graphs are given for different values 
of the parameter €, and the value of £ = 4.5 corresponds to the force maximum 
for a given thickness b = 0.35. 
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Figure 6 Figure 7 


Damping of the Traveling Change in the Phase of the 
Magnetic Field Amplitude Traveling Magnetic Field in 
in the Conductive Layer a Conductive Layer Over a 
Above a One-Sided In- One-Sided Inductor in the 
ductor in the Case of Case of 6 = 0 

6 =0 


It may be seen from the graphs that both the amplitude and the argument 
¢ of the normal induction component decrease as € increases. The phase /26 
of the normal induction component sharply changes based on the thickness of 
the metal layer during this surface effect. When passing through it, a 


rapidly decreases, and the phase of the normal induction component lags behind 
the tangential component on the surface of the inductor. 


In the case of u = œ, H, equals the linear current load Ay: If the 


0 
space is filled by a substance with o' # 0 and y' # œ in the case of x < 0, 
then the following relationship holds between Ho and Ay: 


H, PI 
“ge ~ — 7S [(a? sh a6 + Bch a8) sh Bb + apert ch Bb), ro 


h 
pie d= (a+B)[uB’(ach ad+8 sh aô) + p'a (a sh aô+ $ ch ad) Jes + 


+ (B—a)[yp’ (ach að— B sh ad) + p'a (a shað— pB ch að)]e-8t; 
p’=Ya?+p’0'os. 
The reaction of secondary currents in the conductive layer is character- 
ized by a change in the vector of the normal induction component B,/Bq on the 


inductor surface for 8= 0, as a function of €, and in the relative 
thickness of the conductive layer b (Figure 8). In this case ọọ characterizes 
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the difference between the angle at which the phase shifts between the normal 


induction component for different values of b and e and the tangential com- 
ponent on the surface of the inductor. If e + 0, then Be zi aad 6.4 z . 
Bo 


Both the amplitude of the vector and ọ decrease with an increase in e€. 
e 


In the case of | B > 0, the angle at which the phase shifts between 


Xx 
Bo : 
both components strives to the limiting value r Thus, for b >1 XxX 
B 
0 


strives to the limit from the left side of the bisector, and in the case of 

= 0.1 - 0.2 it strives to the limit from the right side. However, for 

= 0.25 ~- 0.75 the geometric location of the end of the vector has a point 

of inflection approaching the limit first from the right side, and then from 

the left side. In the case of b > 0, the geometric position of the end of the 

vector represents a semicircle. In the case of b = 2, the vector dia- 
Bo 

gram practically coincides with the curve in the case of b = œ. 


N, 


A 
t 
` 


a éso `e 107$ 10° 10°? 10T 10° 10" 10 
Figure 8 Figure 9 
Geometric Location of Dependence of Fa on w in the 


the End of the Vector 7 Case of a0 
By as a Function of e and b 


Bo 
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Let us present the relationships for the components of the ponderomotive 
force in terms of the dimensionless criteria employed in the work (Ref. 10): 


gabsh2y _ sin 2v2 vty (st 2y; ņ sin 2y 
yı Y2 Yı Y2 
+2 (ch 2y; — cos 2y2); 


2 2) 2 ce N a 
R’ = (sh?-y, +sin? y2) [ee ch?3-+sh 26+ sh? i + 


+ 


= 2a? 1 (34) 
+e? $ sh 2y, + 2 sin 2y2+ viv (ch? yı — sin? v2)| ae 


2 e ; 
vv ch de® (yz sh 2yı — yı sin 2y2); 


iare. 1/0, 
aa 9 > X2 y p] 
The transverse component of the force has the following form 


pz _ Exa _ o v 
F= Hf = SR 3 (ch 2yı — cos 2y2) +2 (y2 sh 2yı —yı sin 2y)| 


+ 


£28 


(35) 
where (Fet and Prag are the averaged transverse and longitudinal components 


of the force, respectively, acting upon a volume of substance 2b m - i.e., on 
a column with unit area eut out of the plate . 


Figure 10 


Dependence of R on w in the Case of $ = 0 


Figures 9 and 10 present the dependences Fa and F, according to (34) 


and (35) as a function of w for different values of the parameter b' = b in 
the case of 6 = 0, 


The dependence of the ratio of the maximum values of the force components 
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F and F 


zmax on b/is given in Figure 11. The limiting value of the force component 


Fe in the case of €e > © can be used as an approximate estimate of the maximum 


value of the function F (e). 


If the thickness of the conductive layer is small 


- - i.e., the depth to which the traveling magnetic field penetrates is larger 
than b -- and if the conductive layer is close to the inductor, then F x 


lim F_. 
x 
c> 


Fai 
DLANA 
œg a~ 


OOOO! 
nnen 


Figure 11 


Dependence of F on b 
z max 


x max 


max 


The flux of electromagnet- /29 
ic power through a unit of surface 
of the one-sided inductor also has 
two components LO and Tiz: 


; wsBo? (&, oat a 
Pa=Re M= 50 OR (Ba + p17 +p?) sh 20,+ 


+ h Gtp- 1) sin 2v2 +p (ch 20, — cos 2v); (37) 


2c > _ 3B? PH+) 
P,=Im Th. fer | g 2 sh 2v, — 


_ Wo (2t tho? — , — ra 
tba (tpi? ty- ate A) sin 2v2+ p,? ch 2v, + wz? cos Quef. 


The tangential component of the Poynting vector has only a real part: 


wsB,? — ree 
Il,= anA [1 +p, ch 20; + apo sin 2v2 + 


(38) 


+ (pr?-+ the? + 1) sh? v, + (12+ tp22— 1) cos? va]. 


Here 


vi=apı(6+b—x); ve=arpo(6+b—x); 
R= (sh byi + sin? bye) [p12 + apo)? ch2S + sh 25+ sh? 6]+ 
+e%yp, sh bpi + po sin Wipe + (,2-+ yy?) (ch? by; ~ sin? bye) ]+ 
+ 2pitpe ch 8 eò (pz sh 2p, —H, sin 2bypo). 


I, expresses that portion of the electromagnetic energy which is propa- 


gated in the direction of motion of the traveling magnetic field. 
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The results derived from solving the problems described above may be 
applied in designing electromagnetic chutes for transporting molten ferrous 
metals. Similar problems were solved for electromagnetic mixing of a liquid 
metal in melting furnaces (Ref. 14 - 17, Ref. 20 - 22) and for non-contact flow 
meters for molten metal (Ref. 24, 25). For electromagnetic mixers, the prob- 
lems were solved with many layers, with allowance for an insulation £30 
layer above the inductor and the screen. Since the depth to which the travel- 
ing magnetic field penetrates in these cases is less than the thickness of the 
molten metal layer, it is usually assumed that its thickness is infinite. In 
view of the fact that the nature of electromagnetic processes in the electro- 
conductive layer does not depend on other layers located between the layer 
under consideration and the inductor, the nature of these processes in a molten 
metal is the same as in conductive half-space. 
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Figure 12 Figure 13 
Dependence of the Phase Shift Dependence of the Ratio of the Mag- 
on e€ netic Induction Amplitudes on e€ 


The problem was solved (Ref. 24) for a non-contact electroconduc- /31 
tive medium velocity meter for a system consisting of an inductor, 
an electroconductive layer, and a ferromagnetic packet. It was assumed that 
u = © for the inductor and the ferromagnetic packet 


p= — arctg (th pb tg 2b) (39) 


and 


4/1 40 
x= (ch 2p16+cos 22b. Sa 
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¢ is the phase lag of the normal component of the traveling magnetic 
field when it passes through a layer of the electroconductive medium, and 
x is the ratio of the amplitudes of the normal induction component. 


Figures 12 and 13 present the dependence of the angle at which the phase 
shifts $ and the ratio of the amplitudes upon the parameters £ for the three 
values of b = ba. 


4. Electromagnetic Processes in the Traveling Magnetic Field 
of a Two-Sided Inductor 


Several authors have investigated the problems related to the phenomena 
in the field of a two-sided inductor. The problem was solved in the works 
(Ref. 9, 10) for any value of pọ in a general form. These problems are of 
interest in designing equipment which employs the normal component of the 
electromagnetic force, and for equipment operating in the case of p = 180°. 
They are also of interest for determining the influence upon the tangential 
force component when ọ deviates from 0. 


The problem was solved in (Ref. 10) for one layer (the tangential com- 
ponent of the magnetic field was defined on the layer surface). In this case, 
the following values of the fields hold within the layer: 


He= x SEO [ch B (b +x) +e’? ch B (b — x)] eltwst-az), 


(41) 
= RD pS [sh B (b +x) —e!? sh B (b — x) ] ellost—az), (42) 
E= — e [ch B (b +x) +9 ch B(b—x)] el(wst-a2), (43) 


In this case, the thickness of the conductive layer equals 2b. 


Relationships (41) and (42) may be employed to obtain the equation /32 
for the force lines of the magnetic field, which has the following form for 
no-load operation: 

dx _ cha(6+x)sinaz+cha(b—x)sin(az—wp) 


dz sh a (b +x)cos az—sh a (b — x) cos (az=) a 


In contrast to the angle ¢ which was introduced previously, here Y 
characterizes the phase shift between currents flowing in the direction of 
the y-axis, since it is more advantageous to measure the current instead of 
the magnetic field in practice. Only in two cases does ọ = ¥ hold: when 
they both equal 0, or 180°. Figures 14 and 15 present a picture of the force 
lines of the magnetic field during no-load operation (calculated by Yu. K. 
Krumin) for values of ¥ equalling 0 and 90°, in the case of 2b = 0.5 - i.e., 


T 
when the thickness of the plate is two times smaller than t, since the depend- 
ence on ¥ is more apparent for thin plates than it is for thick plates. 
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In this case, it was found that two components (which are different 
from zero) of the ponderomotive force which is averaged over time and space 
may exist - the longitudinal component which is in operation in the direction 
of motion of the field, and the transverse component which operates in 133 


the direction of the x-axis. 
è 


onho 


Figure 14 Figure 15 
Force Lines of the Traveling Mag- Force Lines of the Traveling Mag- 
netic Field in the Clearance of a netic Field in the Clearance of a 
Two-Sided Inductor in the Case of Two-Sided Inductor in the Case of 
y=0 y = 90° 


The force components may be expressed as follows in dimensionless form: 


_ F 1 = a 
t= -rH = gilt cos 9) Fo+ (1— cos o) Fio), (45) 
where 
F= yosh2yi-+yi sin 2 yz 
= Oa) (ch2y,—cos2 Ya)’ (46) 
Fieo=b —¥2Sh2yi—yisin2 yo l (47) 
=t We y2) (ch 2 yi t cos 2y 
Fa Da (vi2—yo")sh2yisin2y. (48) 


= — 5 = sin p — 
* pHo ? (v1? + y?) (sh? 2 yi tsin? 2y) 
where Eck and Fane designate the averaged transverse and longitudinal force 


components, respectively, acting upon a column cut out of a plate with a 
cross section of 1 unit of area. 


Figures 16-18 graphically illustrate the dependence of (46) - (48) on 
w, but, instead of b, a more convenient criterion is introduced 


ee 

T 
In the case of w «<b, a directly proportional relationship always exists 
between the magnitude of the force and the value of the dimensionless criterion 


W. 
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Figure 16 


Dependence of Fo on w 


In the case of w >b and w >> 1 - i.e., when there is a strongly expressed /34 
surface effect - Fo x Figo v a » where the z-component of force strives 

20 
to zero in inverse proportion to ¥2w . For the z-force component, when the 
inductor is switched on concurrently with Fo there is one characteristic region -- 


in the case of w >> b - when ġ = 0°, but in the case of w «1 the force changes 


in inverse proportion to the first power of w . For Figo? all the maxima /35 


of the curves occur in the case of w >> 2.5. There is no such limit for Fas 


and the maximum value of the force increases with a decrease in b, shifting 
toward smaller values of w simultaneously. When there is a sharply expressed 
skin effect, the values of Fo and Figo coincide, and the maxima strive to the 


value of 0.354. 


The transverse force component Fo also has maxima with respect to w. 
When there is a decrease in b, the maximum value of Fe strives to the value 
sin ọ 
of 0.25. In contrast to Fo Fa is an alternating function of the condition 


w. However, the absolute values of its extremums decrease so rapidly that only 
the region of the first maximum is of practical importance. In the case of 
w >b, we have 


oshV20_ (49) 


If the surface currents are given on the surface of the inductor 
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Figure 17 


Figure 18 
Dependence of Figg On u Dependence of Fe on w 
sin 6 
A=A, ellws!—az) for x= +b; 


A=Ap eflost-az+¥) for ty 


as well as the space behind the surface currents - i.e., if it is occupied by 
a substance with the specific conductivity o' and the magnetic permeability u' 
in the case of |x|>b - then the field strength H, = H_, in the case of x = b 


0 z0 
and Ay are related by the following relationship : 
Ho _ w’BI(w’B+ up’)e™® — (w’B— wp’) e286 — 2up']er (50) 
A, (w’B+ up’)? e — (w’B— pp’)? e289 


(u is the magnetic permeability of the conductive band). 


The angles at which the phase shifts between the fields ọ and the currents 
Y in the inductor are related as follows: 


Im Lin B+ HB’) e — (wB uP’) e~285] ef — 2B’ 
t [(W P+ wp’) e766 — (u’/B— pp’) e-t] — 2 pp" e (51) 
Pe arene L(n’B+ HB’) e — (y’B—uB’)e-**] e9 —2 up’ ” 


(WBF up’) ee — (wP up) e] —2 up’ eit 


In the case of ¥ = 0°, p = Y and 
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Ho 1 


A ae 
° 1422 oth go 
wP (52) 
İn the case of ¥ = 180° ọọ = ¥ and 36 
H, l 
ae a o» 
o 14 EP th po 
uP 
where 
aes: isc, 
(a to’ E Bo 


I. A. Tyutin and E. K. Yankop (Ref. 6) have solved the problem for 
¢ = 0° and the number of layers n = 5, taking into account both the influence 
of the immobile channel walls, and the non-conductive heat-insulating layers. 


Then S4} = Sy = 8, = Se = Ls 
01 =05=0; 02=04= 01; 03=0; 
Bi=Ps=a; Ba=Bo= Bi; B3=B; 
bo=b4= 61; b3=b; bi=bs=b;; 
g=0°; Ha=Ho=H; p. =œ. 
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The normal component of magnetic induction on the inductor surface is assumed 
to be given. 


R. A. Petrovich solved a similar problem, defining the tangential com- 
ponent on the inductor surface. Since it is more advantageous to define the 
tangential component, let us only investigate the solution of Petrovich. On 
the other hand, if we know one solution, we may readily obtain another solu- 
tion by means of the following relationship 


H9=iHgcth a b,, (54) 


where Heo is the amplitude of the normal component of the magnetic field on 


the inductor surface. 


Let us write the amplitudes of the induction components in separate 
regions of the clearance: 


Bzy = Ho Aoi (F chax—Lshax), 


R (55) 


; b b b 
Bzry= bo foi |e shẹ È sh Pe (4-2) +B: ch By ch n(2—a)], (56) 


Berry MOAB at ch px, (57) 
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The density of the induced currents in regions III and IV is: 


Here 
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Avs wo fo [F shax—L chax], 


Brv = pofo Be le sh Bech (s-er $) + 


+B ch BS sh Bi (=-4)]. 


j=- Hgo p: ch Bx, 
x — — Aio b b 
jiv =R [pshp J sh fe (--ż) + 


+ßech Boch Be (:-4)| : 


R=api sh a bich Bibs ch B È +Bê cha bish Biby ch BS + 


b 
+BB: ch a b; ch Beb: sh By tap sha b; sh Bebe sh B t, 


1 
paz | (a+ Be) (Be+B)ch(p Ž + Biba. ~ab) + 


+ (a—B) (B+ Bi)ch(PS + Bib ta +ab,)+ 


+ (api) (Bi—B)ch( -BF +Bbitay +aby) 


+ (a+B1) (B:—B)ch(—B-5 +da $ —aby) 


l l 
L= thers (B+B:)sh(a E +abı—ß ż — Bebe) + 


+ (a—Br) (B+B) sh (BÈ +Bibi+a È +ab) + 


+ 


| 


(58) 


(59) 


(60) 


(61) 


(62) 


(63) 


(64) 


+ (a+r) (Be— B) shp -Bbita $ +ab:) + 
138 


+ (Bi) (B-P sh (Bbp y +a tabs) |. (65) 


The electromagnetic pressure acting upon the conductive band, which is 
averaged over time and over layer thickness and which is referred to a unit of 
layer length, equals 


Pe 


= He? Ac? cows [ef (shee are singt ) (66) 


4b [Ri pı the 
The normal component of the Poynting vector on the surface x = $ has 


the following form 


a=} Eu A = 
(67) 


; 242 
= fous fot (F* chax—L* sh ax) (F shax— L ch ax). 


By computing the active power, we may find the power which is liberated 
in molten metal, and we may determine the losses at the channel walls in the 
form of the difference between the total flux of active power passing through 
the channel wall and the power liberated in them. 


Let us investigate certain simplified cases (Ref. 7). 


If 5, = 0, the electromagnetic pressure of the pump will be 
P MoA lF, Fe 
e 6(Kii+K2)’ (68) 


where & = 2TP a is the length of the pump active zor ; Pa is the number of 


pairs of winding poles; 
1 - 1 f 
Ki=F,; shy a(é—6) +F3ch a a(s—6); 


K2=F.ch ta(ð—b); 
F,=ch 1b + cos 2b; 


oi 
F=- sin pa2b + te sh yb; 
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Fyn tt sh pib — $z Sin 2b. 139 


æ a 


The active power transmitted to the molten metal and the electromagnetic 
pressure are related by the following simple relationship 


Pa= pe Qs, (69) 
where Qs is the output of the pump corresponding to the synchronous velocity 


of the moltem metal 

Q,=2tf ab (70) 
(where a is the width of the pump channel in the direction of the y-axis). 

The electromagnetic power Pa represents the sum of two components: the 

power of the Joule heat losses 

Py=pesQs (71) 
and the mechanical power 

P2=pe(1—s)Qs=peQ, (72) 

where Q is the real output of the pump. 


The reactive power consumed when the magnetic flux passes through the 
clearance is. 


2 1 
P,= BaP ahale, F3cha(é—b) +5 (F2+FiyTF 2) sh a(ò—b)]. (73) 


Let us examine the special case s «<1. In physical terms, this means 
that the magnetic induction produced by the induction currents in the con- 
ductive band is considerably less than the primary (exciting) induction. The 
expressions for the electromagnetic pressure and powers may then be con- /40 
siderably simplified: 


pem Pf otalkr ë (74) 
2_ 29. 
2sh 2 
Pam healclalh eab? B (75) 
aoe 
2sh J 
P,=2rfuAoal cth Š, (76) 


Here we have 
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1 ab+shab 
as ier a a) 


As may be seen from expressions (74 and 75), the electromagnetic press- 
ure and the active power are directly proportional to e in the case of e <<1. 


If t >> b, then the tangential component of the magnetic induction is 
small as compared with the normal component. The latter is almost uniform 
with respect to the clearance height. 


In this case, we have 


— _2poleAo? 
Pe= 0 (8?-+ 626%)’ G8) 
= 2¥0Ao*2rfable (79) 
a(ð? Feb?) ’ 


P,= HoAo?2tfalbs 


f2 2 3) m2 (80) 
=a (82-4 eb) 2e— St (Ite +eVI +e”) ab (6— on}. 


The studies (Ref. 18, 19) solved the problems with respect to electro- 
magnetic processes in a conductive layer having the finite width 2a (i.e., with 
allowance for the transverse edge effect). It is assumed that the magnetic 
field in the clearance of the inductor was plane-parallel (see the article by 
Yu. A. Mikel'son in the present collection). In particular, this work obtained 
the expressions for the electromagnetic pressure p, and the power P,. Assum- 
ing that a = ~ in these formulas, i.e., assuming that the width of the plate 
was infinitely large, they were transformed into 


p= = ufoB? (81) 


l+ė 
and [41 


Po= aon 4ablxof B?, Gy 


i.e., they correspond to the approximation of the solution of R. A. Petrovich 
in the case of 1 >> given above. Here e' =œ b ; B is the induction in the 


S 
clearance during no-load operation. 


The theory for electromagnetic processes in an infinitely wide layer, 
located in a traveling magnetic field, has been developed in greater detail 
than in layers of finite width, since a mathematical analysis is slmpler. 

REFERENCES 


l. Ollendorff, F. Einheitliche Theorie der Drehfeldmaschinen (Unified Theory 
of Rotating Field Machines). Archive f. Elektrotechnik, 24, 2, 1930. 


35 


10. 


11. 


12. 


13. 


14. 


36 


Lopukhina, Ye.M. Issledovaniye asinkhronnogo dvigatelya s rotorom v 
vide pologo tsilindra (Study of the Asynchronous Engine with a Rotor in 
the Form of a Complete Cylinder). Elektrichestvo 5, 26, 1950. 

Tyutin, I.A. Mekhanicheskiye sily v begushchem magnitnom pole (Mechanical 
Forces in a Traveling Magnetic Field). V kn: Voprosy energetiki (In 
the Book: Problems of Energetics), 3. Izdatel'stvo AN Latv. SSR, 1955. 

Tyutin, I.A. Vvedeniye v teoriyu induktsionnykh nasosov (Introduction to 
the Theory of Induction Pumps). Trudy Instituta Fiziki, AN Latv. SSR, 
8, 1956. 

Tyutin, I.A. Elektromagnitnyye nasosy dlya zhidkikh metallov (Electro- 
magnetic Pumps for Molten Metals). Riga, 1959. 

Tyutin, I.A., Yankop, E.K. Elektromagnitnyye protsessy v induktsionnykh 
nasosakh dlya zhidkikh metallov (Electromagnetic Processes in Induction 
Pumps for Molten Metals). Trudy Instituta Fiziki AN Latv. SSR, 8, 1956. 

Liyelpeter, Ya.Ya., Petrovich, R.A. K teorii ploskikh induktsionnykh 
nasosov (Theory of Plane Induction Pumps). Izvestiya AN Latv. SSR, 
Seriya Fizicheskikh i Tekhnicheskikh Nauk, 1, 1964. 

Liyelpeter, Ya.Ya., Tyutin, 1.A. Metodika rascheta induktsionnykh nasosov 
dlya zhidkogo metalla (Method for Designing Induction Pumps for Molten 
Metal). Trudy Instituta Fiziki AN Latv. SSR, 8, 1956. 

Mikel'son, Yu.A. Provodyashchiy sloy v begushchem elektromagnitnom 
pole dvukhstoronnego nesimmetrichnogo induktora (Conductive Layer in a 
Traveling Electromagnetic Field of a Two-Sided Non-symmetrical Inductor). 
Izvestiya AN Latv. SSR, Seriya Fizich.i Tekhnicheskikh Nauk, 2, 1965. 

Veze, A.K., Krumin', Yu.K. Ob elektromagnitnoy sile, deystvuyushchey 
na beskonechno shirokiy provodyashchiy sloy v begushchem magnitnom pole 
ploskikh induktorov (Electromagnetic Force Acting Upon an Infinitely 
Wide Conductive Layer in a Traveling Magnetic Field of Plane Inductors). 
Magnitnaya Gidrodinamika, 4, 1965. 

Kirko, I.M. Kriterii podobiya elektrodinamicheskikh yavleniy pri otnos- 
itel'nom dvizhenii magnitnogo polya i provodyashchey sredy (Criteria 
for Similarity in Elektrodynamic Phenomena During the Relative Motion 
of a Magnetic Field and a Conductive Medium). V kn: Voprosy energetiki 
(In the Book: Problems of Energetics), 3, 1955. 

Veze, A.K., Liyelausis, 0.A., Petrovich, R.A., Ulmanis, L.A. Provodyashiy 
sloy v begushchem elektromagnitnom pole odnostoronnego induktora (Con- 
ductive Layer in a Traveling Electromagnetic Field of a One-Sided In- 
ductor). V kn: Voprosy magnitnoy gidrodinamiki (In the Book: Problems 
of Hydrodynamics), III. Izdatel'stvo AN Latv. SSR, 1963. 

Veze, A.K., Ulmanis, L.Ya. Raspredeleniye elektromagnitnogo polya i 
ponderomotornykh sil v beskonechnoy provodyashchey polose, pomeshchennoy 
v begushchem magnitnom pole odnostoronnego induktora (Distribution of 
the Electromagnetic Field and Ponderomotive Forces in an Infinitely 
Conductive ‘Band Placed in a Traveling Magnetic Field of a One-sided 
Inductor). In press. 

Vaynberg, G.S. K teorii ustroystva dlya elektromagnitnogo peremeshivaniya 
rasplavlennogo metalla v dugovykh elektropechakh (Theory of Equipment 
for the Electromagnetic Mixing of Molten Metal in Arc Electric Furnaces). 
Elektrichestvo, 2, 1958. 


15. 


16. 


17. 


18. 


19. 


20. 


21. 


22. 


23. 


24. 


25. 


26. 


27. 


28. 


Vaynberg, G.S. 0O vybore chastoty ustroystv dlya elektromagnitnogo per /42 
emeshivaniya metallov v elektropechi (Selection of Equipment Frequency 
for Electromagnetic Mixing of Metals in an Electric Furnace). Elektri- 
chestvo 5, 1958. 

Kochnev, E.K. K teorii ustroystv dlya elektromagnitnogo peremeshivaniya 
rasplaviennogo metalls (Theory of Equipment for Electromagnetic Mixing of 
Molten Metal). Elektrichestvo, 7, 1959. 

Kochnev, E.K., Rezin, M.G. Issledovaniye ustroystva po elektromagnitnomu 
peremeshivaniyu rasplavlennogo metalla (Study of Equipment for Electro- 
magnetic Mixing of Molten Metal). Izvestiya Vysshikh Uchebnykh Zaved- 
eniy, Elektrotekhnika, 9, 1962. 

Vol'dek, A.I. Toki i usiliya v sloye metalla ploskikh induktsionnykh nas- 
osov (Currents and Stresses in the Layer of Metal of Plane Induction 
Pumps). Izvestiya Vysshikh Uchebnykh Zavedeniy, Elektromekhanika, 1,1959, 

Yanes, Kh. I. Uchet vliyaniya vtorichnoy sistemy v lineynoy ploskoy mag- 
nitogidrodinamicheskoy mashine (Influence of the Secondary System in the 
Linear Plane Magnetohydrodynamic Machine). Trudy Tallinskogo Politekhn 
Instituta, Seriya A, 197, 1962. 

Vortnichuk, N.I., Krutyanskiy, M.M. 0O vybore optimal'noy chastoty toka 
statora pri peremeshivanii zhidkoy stali s pomoshch'yu begushchego mag- 
nitnogo polya (Selection of the Optimum Frequency of the Stator Current 
When Molten Steel is Mixed by Means of a Traveling Magnetic Field). V 
kn: Voprosy magnitnoy gidrodinamiki i dinamiki plazmy (Problems of Mag- 
netic Hydrodynamics and Plazma Dynamics). Izdatel'stvo AN Latv. SSR,1959. 

Ostroumov, G.A. 0O peremeshivanii rasplavlennykh metallov begushchim magn- 
itnym polem (Mixing of Molten Metals by Means of a Traveling Magnetic 
Field). V kn: Voprosy magnitnoy gidrodinamiki i dinamiki plazmy (Prob- 
lems of Magnetic Hydrodynamics and Plazma Dynamics). Izdatel'stvo AN 
Latv. SSR, 1959. 

Ostroumov, G.A. Fiziko-Matematicheskiye osnovy magnitnogo peremeshivaniya 
rasplavov (Physico-Mathematical Bases of Magnetic Mixing of Melts.). 
Metallurgizdat, 1961. 

Schilder, J. Pohybujici se elektromagnetické pole ve vodivém prostredí 

Pe ce romane i Field Motion in a Conductive Medium). Elektrotechnicky 

Ulmanis? L?ya. Beskontaktnyy raskhodomer dlya zhidkikh metallov (Non-Con- 
tact Flow Meter for Molten Metals). Byulleten' Izobreteniy, 19, 1962. 

Ul'manis, L.Ya. Fizicheskiye yavleniya pri induktsionnom vozdeystvii beg- 
ushchego magnitnogo polya na sloy zhidkogo metalla (Physical Phenomena 
During Induction Interaction of a Traveling Magnetic Field on a Lay- 
er of Liquid Metal). Avtoref. Kanddiss. (Author's Abstract of Candi- 
date's Dissertation). Riga, 1963. 

Ul'manis, L.Ya. K voprosu o krayevykh effektakh v lineynykh induktsionnylkh 
(Problems of Edge Effects in Linear Induction Pumps). Trudy Instituta 
Nasosakh Fiziki, AN Latv. SSR, 8, 1956. 

Watt, D.A. A Study in Design of Traveling Field Electromagnetic Pumps 
for Liquid Metals. Harwell, 1955. 

Valdmanis, Ya.Ya. Elektromagnitnyye sily deystvuyushchiye na beskonechn- 
uyu provodyashchuyu polosu v pole odnostoronnego induktora trekhfaznogo 
toka (Electromagnetic Forces Acting Upon an Infinite Conductive Band in 
the Field of a One-Sided Inductor of Tri-Phase Current.) Izvestiya AN 
Latv. SSR, Seriya Fizich. i Tekhnicheskikh Nauk, 1, 1965. 


37 


HIGHER SPATIAL HARMONICS OF THE MAGNETIC FIELD OF AN 143 
INDUCTION MHD MACHINE 


Yu. Ya. Mikel'son 


1. Introduction 


In terms of their principle of operation, induction MHD machines are 
similar to asynchronous electric machines. Just as in asynchronous engines, 
the winding of the stator produces a traveling magnetic field, under the in- 
fluence of which currents are induced in the molten metal of the MHD machine 
(in the rotor of the asynchronous machine). The interaction of these currents 
with the magnetic field of the stator leads to the formation of ponderomotive 
forces. Therefore, there is a certain similarity between certain problems in- 
volved in the theory of asynchronous engines and induction MHD machines. One 
of these problems is the deviation of the traveling magnetic field from a 
sinusoidal one. This deviation is caused by several factors. The most im- 
portant factor is the spatial distribution of the stator winding and the non- 
uniformity of the air gap (projections and grooves on the steel surface of the 
stator). These factors occur both in asynchronous engines and in MHD machines. 
However, they may have a different influence upon the operation of the mach- 
ine due to differences in both devices. We would like to emphasize the follow- 
ing differences: 


(1) The "rotor" of the MHD machine represents a solid, liquid conductive 
medium, in contrast to the discrete conductors of the rotor winding of the 
asynchronous machine; 


(2) The clearance between the rotor and the stator of the asynchronous 
machine may differ greatly from this clearance in the induction MHD machine; 


(3) The stator of the induction MHD machine may be both one-sided and 
two-sided; 


(4) The operational regimes of the machines under consideration are 
different. As a rule, the asynchronous engines operate with slipping which is 
close to zero. Not all of these differences play the same role. An increase 
in the clearance between the rotor and the stator in a MHD machine, as compared 
with an asynchronous machine, reduces the role of the higher spatial harmonics 
(h.s.h.), while the wide region of slippings and the solid "rotor" of 144 
the MHD machine may lead to an increase in the influence of the h.s.h. upon 
it. 


Many authors have investigated the harmonics produced due to the spatial 
distribution of the multiphase stator winding, and have determined their in- 
fluence upon the motion of the rotor. The results derived from a large number 
of these works have been generalized in (Ref. 1). 
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Figure 1 


Infinite Layer of a Conductor (Above) and a Magnetic Current (Below). The 
Linear Current Load (Ref. 12) is Given on the Surface of the Magnetic Circuit. 


The spatial distribution of a winding of multiphase current leads to 
the fact that a magnetic field is produced which consists of an infinite series 
of harmonics traveling in the opposite directions. In particular, harmonics 
of the orders 6k + 1 occur for a symmetrical, three-phase winding, where 
k = 0, + 1, +2,... The positive and negative harmonics have different directions 
of motion. The number of grooves per pole and per phase and the contrac- 
tion of the winding step are described by means of the so-called winding co- 
efficients kyin . In the majority of cases, the h.s.h. are disregarded when 

v 
the electromagnetic processes in a MHD machine are investigated (Ref. 2-5). 
In many cases, this is valid - for example, in the case of large non-magnetic 
clearances, and a sufficiently thick conductor layer. Due to the smallness 
of the clearance between the rotor and the stator, the h.s.h. play a signifi- 
cant role in asynchronous motors. Therefore, special measures should be taken 
to suppress them. With respect to MHD machines, the influence of the h.s.h. 
in a molten metal and in a stator upon electrodynamic force density and energy 
losses has not been studied sufficiently (Ref. 6-11). 

I. M. Postnikov (Ref. 12) has presented a method for computing the losses 
from h.s.h. in the rotor of an asynchronous machine. His method may be employ- 
ed to calculate the losses in a MHD machine with a one-sided inductor for a 
sufficiently thick layer of molten metal (theoretically infinite) (Figure 1). 
However, this method is not applicable for calculating losses in a MHD 145 
machine with a two-sided inductor. In the case of the two-sided inductor, 
the layer of molten metal has a finite thickness in contrast to the infinitely 
thick layer (Ref. 12). The windings of the stator with magnetic circuits are 
located on both sides of the molten metal layer. 


met 
The electromagnetic fields in the nolemjot a MHD machine may differ 
greatly in the case of a two-sided and a one-sided inductor of the magnetic 
field (Ref. 6). In addition, the determination of the expansion coefficients 
(Ref. 12) of linear current loading in series with respect to individual har- 
monics, when there are grooves on the surface of the stator, is an independent 
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problem (Ref. 18). These coefficients were assumed to be known on the basis 
of the theory of electric machines in (Ref. 12). 


The influence of h.s.h. on the density of the electrodynamic forces and 
Joule heat losses in the model of a MHD machine with a two-sided symmetrical 
inductor was investigated in (Ref. 10). It was shown in this report that in 
several cases the h.s.h. makes a significant contribution to the enérgy losses 
and to the density of the electrodynamic forces. 


The influence of a non-uniform air gap upon the motion of the rotor of 
an asynchronous machine was also investigated in (Ref. 1). This problem may 
be reduced to determining the induction of the magnetic field on the surface 
of a smooth medium (rotor) with a magnetic permeability of u = © and electric 
conductivity of o = 0, when there is one infinitely deep groove (Ref. 13) or 
an infinite series of grooves (Ref. 14) in another medium (in the steel of the 
stator) also with uy = œ and o = O (Figures 2 and 3). 
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Figure 2 


Infinite Groove in a Magnetic Circuit Opposite a Smooth Magnetic Circuit 
(Ref. 1) 


For the theoretical computations, it is assumed that the dimensions /46 
of the portions of the systems shown in Figures 2 and 3 along the z-and y- 
axes are infinite. 


The difference of the magnetic potentials between a smooth and 
serrated medium with = œ and o = 0 is assumed to be constant. The problem 
may be solved by the method of conformal mapping on the basis of the theory of 
the scalar magnetic potential. The presence of grooves leads to a decrease in 
the average magnitude of the induction on the smooth surface of the "rotor". 
This is taken into account by the Carter coefficient in the first approximation 
in practical applications (Ref. 1). The exact distribution of the field in the 
clearance is determined in (Ref. 14). However, in view of the cumbersome nature 
of the final results, the amplitude of separate harmonics corresponding to the 
projections was calculated only in this special case by Freeman (Ref. 16). 
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Figure 3 


Infinite System of Grooves and Projections Opposite a Smooth 
Magnetic Circuit (Ref. 14) 


A.I. Vol'dek (Ref. 15) has provided a simplified method for calculating these 
harmonics. For example, S. P. Pecheritsa (Ref. 17) has employed his methods. 
However, even for asynchronous motors it is impossible to assume that the prob- 
lem of serrated pulsations has been studied exhaustively. The problem remains open 
regarding the distribution of the field within the clearance, and not only on 
the surface of the "rotor", regarding the influence of the depth and width of 
the grooves upon this distribution, as well as other problems. A conductive 
medium is located within the clearance in a MHD machine. The distribution of 
the field in this medium, taking the non-uniformity of the clearance into ac- 
count, is also of interest. This distribution influences the density of elec- 
trodynamic forces and energy losses in an induction MHD machine. A theoretical 
investigation of this problem was performed in (Ref. 18), and some of the re- 
sults derived in this study will be investigated below. 


2. H.s.h. From the Spatial Distribution of the Winding with [47 
a Smooth Steel Surface of the Stator 


There are several factors which complicate a theoretical investigation 
of the influence of h.s.h. upon the electromagnetic processes in an induction 
MHD machine. Any real MHD machine has finite dimensions, which completely de- 
termine the properties of the materials from which the channel of the MHD mach- 
ine, the stator, and the winding are made. The velocity of the molten metal in 
the channel of the MHD machine is not constant over the cross-section. It is im- 
possible to make strict allowance for all of these factors. Therefore, several 
simplifications must be employed in theoretical computations. The finite dim- 


ensions of the MHD machine are disregarded in an investigation of the role of 
h.s.h. 
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Let us investigate the simplified model shown in Figure 4. 


Figure 4 


Non-Symmetrical Model of a MHD Transformer with Conductive 
Channel Walls 


Here, region I is the working component of the MHD machine with con- 
ductivity o and magnetic permeability Ho: The regions II and III represent the 


non-magnetic clearance. 


The winding of the stator is usually located in special grooves of the 
stator steel. The grooves make the non-magnetic clearance of the MHD machine 
non-uniform, but we shall first investigate the smooth steel surface of the 
stator (the influence of the grooves will be investigated in section three). 
In this connection, let us place the winding of the model of the stator (see 
Figure 4) in the clearance between the channel of the MHD machine and the sta- 
tor. By changing the distances dy and d,» we may obtain any arrangement of 


the winding in the clearance. Thus, in the case of di = 2i and d, = Lo the 


winding lies on the steel surface of the stator. The winding (only one phase 
is shown) consists of q linear plates having the thickness h (q grooves per 
pole and per phase). The distance between the plates in one phase is ti - h 


(the width of the "projection" is t, - h). The distance between the for- /48 
ward and inverse currents of one phase is à; the polar step is t; (the ratio 


à characterizes the contraction of the winding step). 
T 


Regions IV and V in Figure 4 are the channel walls of an induction MHD 
machine with conductivity ou magnetic permeability uo? and thickness 85 and 


6,. 
2 
every point and equals v(0,0,v). It is assumed that the current is a sin- 
usoidal current with the angular frequency w, and is uniformly distributed 
over all plates of one phase. The current amplitude of one phase of the wind- 
ing in region II is Io]; in region III it is Ig9. The phase shift between 


It is assumed that velocity of the medium of region I is identical at 
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the currents in the upper and lower inductors is ¢. A medium with yu = œ and 
o = 0 is located behine the winding. The other geometric dimensions are 
shown in the figure. 


After a theoretical investigation of the model of the MHD machine shown 
in Figure 4 (see appendix 1), we may find the z-component (which is averaged 
over time and the coordinates) of the force density 


R=) fev (1) 
Vv 
the x-component (which is averaged over time and z-coordinate) of the force 
density 


(2) 
where v = 6k + 1, k = 0, +1, +2, ...3 f, gives the force density of the yth 
harmonics of the traveling magnetic field. oy may be positive or negative, 


fz= $ fev (x), 


depending on the slipping. The mean Joule heat losses Qyoule in the chan- 


nel of the MHD machine may be obtained by means of the x-component of the 
Poynting vector: 


Re S,= + fx. 
‘0 


(3) 
The energy losses in the molten metal* are 
(4) 
Qio 1 =Re [Szi lx=5— xl |==] 
in the channel walls 
Q J% =Re [Sziv le=ot+01— xIV zap), (5) 
Q jy =Re [Sev l=- — Sev Iz=), (6) 


where Re designates the real part. 


In numerical calculations of the force density and Joule heat loss- /49 
es in an induction MHD machine, it is advantageous to introduce the relative 
force density and Joule losses: 


26 
I a IE’ (7) 
b 
T =Q, TA (8) 


The contribution made by the h.s.h. to the force density and Joule heat 
losses is the difference between the sum of the series and the first term of 


* Energy losses in a vertical column with unit transverse cross-section, 
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the series (1), (2). 


The computations of the force density and Joule heat losses according to 
formulas (7) and (8) are cumbersome. The computations were performed on a 
BESM-2 in the Computing Center of Lenningrad State University imenii P. Stuchok 


for the special case when 
6:=b2=0, d;=d2=l=h=d, =q, h=0, 


=0, Io=1, =I, bat =}: 2 
F a=lo=lo t 3g’ q=1; 3; 5, 
and it was found that the contribution made by the h.s.h. to the force den- 
sity and Joule heat losses may be quite significant for specific values of the 
parameters. In this special case, i and Qo were expressed by means of the 


following dimensionless parameters: 


(9) 
UT 
s=] -iz and. 


The computations of E and q were performed for a constant linear cur- 

rent density Jo = 319 `Z Y . We shall present some results derived from the numeri- 
eae 

cal calculation. Thus, Table 1 presents the relative force density in the case 
of € = 1, € = 20, q = 1 and different values of ¢ and s. Table 2 presents the 
relative Joule heat losses in the channel of a MHD machine for the same values 
of the parameters. As may be seen, the contribution made by the h.s.h. to the 
force density and the Joule heat losses depends on the paramteter ~. The role 
of the h.s.h. is greatest for small values of the parameter ¢. Thus, for ex- 
ample, in the case of g = 107 3, & = 20, e€=1,q=1,s=1, Ê zl comprises 


only 72% of fo The situation is the same for Qs and Os 1 In the case of 


& > 0.1, it is practically impossible to take into account the influence of 
the h.s.h.. It must be noted that in the case of s = 1 - i.e., when the molten 
metal is immobile - the contribution made by the h.s.h. is minimal. 


It thus follows that it is impossible to estimate the influence of /51 
the h.s.h. in the case of s = 1; it will be too low. The tables present 


fz and I in the case of q = 1, which corresponds to the case in which the 


contribution made by the h.s.h. is greatest. For small ¢, q > 1 in real MHD 
machines. An increase in q leads to a reduction in the role of the h.s.h.. 
This may be readily seen from the graphs shown in Figures 5 and 6. 


Thus, for example, Figure 5 presents a graph showing the dependence of 
the dimensionless force density fo on the slipping (sum of the series and first 


term of the series) for the following values of the dimensionless parameters 
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1072, e = 3.72, Ẹ = 10, q = 1; 335. 


(Ref. 10) 


Figure 6 shows the dependence of the Joule heat losses on the slipping 


It should again be noted 


that the first term of the series of the force density and the Joule heat loss- 


es differs theleast from the sum of the entire series in the case of s= 1. 


for the same values of the dimensionless parameters. 


1 the entire force density practically co- 


incides with the first term of the series (1), whereas in the case of s = -2 


the case of s 


in 


Thus, for example, 
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Figure 5 Figure 6 

1 - First Term of the Series 1 - First Term of the Series Q 
f | in the Case of q = 1; as Mae 
zl 3 in the Case of q = 1; 2 - Q, in the 
2--£ in the Case of q = 1; = Jo 

“21 z Case of q = 1; 3 - Qo in the Case 
7 fa a Ames) of q= 3; 4- Q50 in the Case of 
4 - f in the Case of q = 5. TE 


(generator regime) and in the case of q = 1 the difference is a factor of four. 
In the case of q = 3 the difference is a factor of 1.62, and in the case of 
q = 5 -- 1.24 (see Figure 5). The situation is the same for the Joule heat 


losses (see Figure 6). 


It follows from the tables and graphs presented above that allowance /52 
for only the first harmonics of a traveling magnetic field does not always 
lead to a correct result. A rigorous distinction must be drawn between the 
region of the parameters (¢, S, q, etc.), where it is possible to disregard 
the influence of the h.s.h. and where it is impossible to do so. 


3. Influence of an Inductor Surface Having Projections on the 


Electromagnetic Field Distribution in a Conductive Band 


The grooves and projections on the steel surface of the stator deform 
the electromagnetic field in the molten metal of a MHD machine. This distorted 
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field may be expanded in series and we thus obtain the so-called harmonics 
“corresponding to the projections" (Ref. 15). However, this expansion is 
rendered more difficult by the fact that the "distorted field" is unknown, 
which must be expanded in series. In the air gap of asynchronous machines, 
this field may be found under the following assumptions (Ref. 14): 


(1) The rotor and the stator have the conductivity o = 0 and magnetic 
permeability u = œ (see Figure 3); 


(2) The grooves on the surface of the stator are infinite or finite in 
depth; 


(3) The rotor may be smooth or may have grooves; 


(4) A constant difference of the magnetic potentials is given between 
the rotor and the stator; 


(5) The dimensions of the system (Figure 3) along the y- and z-axes 
are infinite. 


The definitive results derived from the study by Coe and Taylor (Ref. 14) 
are not suitable for numerical calculations. In practice, the method of Vol'dek 
(Ref. 15) is employed more frequently to determine the amplitudes of the har- 
monics corresponding to the projections. In this method, the grooves on the 
rotor surface influence the magnetic permeability of the air gap. The expansion 
in series of the approximate periodic curve of the magnetic conductivity leads 
to the harmonics corresponding to the projections. If allowance is only made 
for the first term of the expansion, this produces a decrease in the average 
induction in the gap due to the non-uniformity of the air gap by a factor of 


ks -~ where ks is the Carter coefficient (Ref. 1). For the most part, in prac- 
tice the non-uniformity of the air gap is taken into account by means of this 
coefficient -- i.e., it is assumed that it is equivalent to a certain increase 


in the uniform gap. 


There is a conductive medium in the air gap in MHD machines. The electro- 
magnetic field in this medium differes greatly from the electromagnetic field 
in the air gap of an asynchronous machine. The depth of the grooves on the 
steel surface of the stator has a finite value. It is not sufficient to take 
into account the non-uniformity of the air gap in the MHD machine by means of 
the Carter coefficient. 


In this connection, the study (Ref. 18) discussed an approximate 153 
model of the MHD machine (Figure 7) in order to allow for the influence 
of the non-uniformity of the air gap upon electromagnetic processes in the MHD 
machine. This model has infinite dimensions along the y- and z-axes. Regions 
III and V have the conductivity o = 0, and the magnetic permeability pu = uo’ 


Region IV is the conductor with the`conductivity o and magnetic permeability 
Ups Region I is a groove occupied by an alternating current with constant den- 


sity amplitude jg and with angular frequency w; regions IT and II' are the 
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empty grooves. The vertical lines z = - 5 and z= = divide half of the period 


(the system has a period along the z-axis equalling 21). The geometric dim- 
ensions of the regions are shown in the figure. 


Figure 7 


Model of the MHD Transformer with a One~Sided Inductor, with 
Allowance for Non-uniform Gap 


The empty grooves are left for loops of the winding for the second and 
third phases. This system corresponds to a one-sided inductor having a tri- 
phase current with the number of grooves per pole and phase equalling unity. 
The solutions obtained make it possible to write the expression of the vector 
potential for a two-sided symmetrical inductor (Ref. 18) and in other special 


cases. 


According to (Ref. 18), the vector potential A(0,A,0)* in different 
regions (see Figure 7) has the following form. 


an chAm(x+En) 
= -7 9 m nj 4 
A 2 2 (2Kxprcimn + pojoln Omo) ehh hate COS Amz (10) 
1, 
—gHofo(*tln)?, 
54 
ch > hn (2+ Ln) 
Qe n 
Au= KrqiBm (ee ~3z_1 
u Dey rOnB an cos mm Sa ol (11) 
2 
œ œ ch Shin (+n) 
Aw= XE Kapama —j——— cos am (£432 1), (12) 
m=) k=0 ch Jimla 3h 2 


* The dependence on time ett is excluded in all formulas for the vector 
potential. 


48 


An= > Kf ch Om(X—8+6)ch ymb + 
m=0 


+e sh am(ôð—a-— b) chan(x—a) — 
Ym 


ch amô ’ 


-1 cham (a-+b—6) sh am(x—a)| sh mmo | E Omz 
m 


Ay= È Kn [en am(a+b—ò)chym(a+b—x)— 


æ 
—-" sh m(O— =p C E COS AmZ 
a ee 0) | Sra 


Ava > je Stan (= 8) coia, 
m=0 


ch amô 
where 
2r yi 
=z ™ am = y (2m+1); 
Ym =¥ Yam? + iopow; 
Pm= {en am (b — ò) ch ymb + Es sham(d—a—6)chana+ 
Y 
+—2 ch am(ò—a—b)sh ——; 
a, m( )s ana | sh ynd) cir aad? 
Ynm{sh am(6—8)ch vnb- | TE a=) E 
m 
Y ' l 
+= cham(ð—a—b)ch ——; 
am m( ye ama |sh vmb | aoa 
Omk= m? cos AM sin rep; 
1 4 
Bma= — A gat ga [0s am sin a (i, +E) + 
+sin x («4-7 '); 
3 
h 
er= z; (4+1). 
The coefficients Ka are determined from an infinite system of 
tions 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) /55 


(21) 


(22) 


equa- 
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œ 
2 KmOkm= 2pojolnXrs (23) 


where a 1 
Orm=m Pe (iaman th Avla + BymBva th > Ay n) z (24) 
-4 mY mõrm 3 
x= snan ; i 
saal 0; km, 
‘(1 k=m. (26) 


Thus, the determination of the vector potential may be reduced to sol- 
ving an infinite system of equations (23). Let us investigate certain spec- 
ial cases (Ref. 18). 


l. There is no conductive band IV in the clearance. In order to de- 
termine the vector potential in this case, we may set o = 0. Asa result, 
we have 


Qn= 1 ; Pm = th Omd, Ymn>= Gm, 
z N | 
Okm= > A fey m@yk th slat By mBve th oO Ayla + 
vol 
+5 adm th AmÔÔkm, 


ch Am (x—6) 


COS Am2. 
ch amd i 


œ 
Amn=4y=4Ay= 2, Kn 
m=0 


(27) 
II, In the case of the two-sided symmetrical inductor, the ex- {56 
pression for the vector potential may be written as follows: 


TE 
u 2 Ku ch am(x—a)ch ymb (28) 


Ym COS mZ 
—_—7 sh az m: 
a, SP am (x—a)sh vm] chad | 


A 5 COS AmZ (29) 
w= Kn, ch ym(8—x) amt 

m=0 i Ym( ) ch Gn 
The coefficients Kn are again determined from system (23), where 


(30) 


= T In z : 
Pm fen Oma ch Ym (ò —a) + an Sh ama sh Ym (ð ola amô’ 
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Ym=|—sh ama ch ym(d—a@) + in ch ama sh ym(6— a) 


ch amô’ (31) 
Formulas (28) -- (31) were obtained under the assumption that a + b = 
= 6 -- i.e., in the case of a two-sided symmetrical inductor, the horizontal 


plane of symmetry passing through the middle of the conductive band corres- 
ponds to the boundary of the medium with u = œ and o = O in the case of x = ô. 
The induction vector of the magnetic field B = rot A is perpendicular to this 
plane of symmetry. 


If there is no medium with u4 = œ and o = O above the conductive band, 
we have the following expression for the vector potential (Ref. 18): 


fe] 
Ain= J) Kme-em(2*6) {e-tnt —)ch ymb + 
m=0 


+f chan(x—a) — Yæ sh am (<—a)] sh vnb) COS On2, (32) 
Ym Qin 


æ 
Ay= }, Kme~amla+b) [ch yn(at+6—2) ~ 


m=0 
-$2 sh ym(x—a~6) | cos Amz, (33) 
Ym 
ao 
Ay= }, Kme-%m* cos amz. (34) 
mal 
The coefficients Ky as always, are determined from the system 57 
(23), where 
@m= e-tm'a+d) [esns ch ymb + ke ch ama + Ta sh ama] sh vnb) i (35) 
m m 


Pm = — e~2mla +b) fazne ch ymb + [fe sh ana + Ya ch ana] sh vmb} 
m m 


(36) 


Thus, it is possible to find the vector potential in other special 
cases. It must be pointed out that when the depth of the grooves ka strives 


to zero -- so that the total current in the groove does not change —- i.e., 
: I 
joln= 7 =const, 


the expression for the vector potential changes into a similar one for an 
inductor without grooves and projections (Ref. 6-10). 


Summing up the vector potential of three phases, we may obtain the 
traveling magnetic field, with allowance for the non-uniformity of the air 


51 


gap. 


In the final analysis, determination of the influence of grooves and 
projections may be reduced to an investigation of an infinite system of equa- 
tions (23). This system requires additional analysis, in order to establish 
the rapid decrease in the coefficients Ka when m increases, and the analysis 


requires additional computational work. It may be performed most readily on 
computers. 


APPENDIX 


4, Results of Analytical Computations 


The results derived in (Ref. 11) are employed to determine the distrib- 
ution of the force density and the Joule heat losses in the channel and in 
the walls of the model of the induction MHD machine (see Figure 4). A sol- 
ution of the equation for the vector potential A(0,A,0) was obtained in this 
study: 


A,=3 ), [a,Me-WF + ay Mevyt]ellwsyt—ayz’) 
y 
Ayy=3 5 [gyDe-W g, Oev] eot-a,2, 
v 
Ay=3 È [c,Me-Vit+-6,ev,2] elot—a,z), 
v 
where i UT 
2’=z-vl; sy=l——v; 
WT 
= Lid . —_ . 2 H . 
a=} ay av; Y= YVay?+ iopowsy; 
yy =Va,? tioo; v=6k+1; k=0, +1, +2,... 


<(Gyly—M,Ny) ; 
a, =K olor ch ay (h —d,) My + poloze’? ch a, (1a—d2) Gv. 


(M,N, —G,L,) 


gam 2 (1 + 2) a, (Det -yYy)? +(1-7 = a aeniea]; : 


y 
g= all z ijen ay Mey try) 4 +( 1+ ‘a voeni]; 


— 


* The z' axis is connected with the moving medium I. 
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cy) =} (: + 2) ay Verve) 4 ( < phoeni] ; 
9 


i% 
e= 2 ((1- J a, Devyty b4 ( 1+ 5) adet; —y,»»|; 
Le 
+ 
M,= Y apan [æv sha „ (b+ô2—l2) ak 


-f cha, (b+d:—h)]+ ŽZÝ oe eWay sha, (6+ 82—l2) + 
+7, cha, (b+82—t2)]; 
+ $ 
Gy -Sir e~ (¥yo+ vy 


y 


2) fa, shay (1; —b—8)) — 


-rcha (1-b—0:)]+ 2” ae [ev sh æy (l1 —b—8,) + 


p +e cha, (li —b— 8:)]; 
vory -vð 
Ny= yet vd [ash a, (f; -b—8)) —yscha, (fy —b—,)}+ 


gtr : 
FE KAW giat, [a, shay (fi — b — 1) +y5 ch a, (l1 ~b —8)]; 


Y-Y 
L =~ E ? @Vyer-Vy? [a sh &, (b +82— l2) — y5 ch ay (b +82—l2)]+ 


Kt _ e 
te WP F129) fa, sh as (b+ b2— l2) +y5 cha, (b +ô2—la)]. 


The winding coefficient Kyiv of the harmonics of the order v has the follow- 


ing form (Ref. 10): 
K_ .= Kyu Ku Kn, 


wiv 


where the coefficient of the winding distribution is 


tt, 
sin —! 


2t 
mvt,’ 
2t 


v= 
qsin—— 
the contraction coefficient of the winding step is 


TAv 


Kiv oT , 
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and the coefficient depending on the width of the inductor plates (winding 
coefficient "of the groove opening" [Ref. 15]) is 


27 . «vk 
Kin = xvh Sin gg 


The complex expressions of the magnetic field induction B and the curr- 
ent density j may be readily obtained according to the following formulas 


OA 


B= ane: 
rotA, j= os. 


According to (Ref. 11), for region I we have: 


B = 3 = , 
7 2 Yy [a, evs — a, Me v1] ellos taz 2, 


By= 3iy, a [a 2e- 4 aP evy] ellwsyt—a,2") 
Yy 
J £) Blea), Sv [ay De- pa Dev, x] ilosta, 


The expressions for the induction B and the current j in regions IV and /60 
V are similar to the corresponding expressions in region I. However, aang) 
stead of the coefficients a, and a; 2) we must substitute 8.) and 8, 
(region IV) or e; (region V) instead of z' - z, o - ooe YT YS . Alli sS; 
must be set equal to unity. 


Employing the expressions obtained for the magnetic field induction 
and current density, we may find the force density f which is averaged over 


time and the coordinates (Ref. 9): 


a= È faw 


j x= > F a(x, 
where 


_ 90w la ©’) sh (2b Rew) 
Faw y avsy [2 (DP ( a ' GER at 


eo sin (2 Im m yb) 


90w la? 
h= mg Yv” Sv Re Za)? om | e? Reve — 


* fj is averaged over time and the z-coordinate. 
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—en2Rey,x —2i Im(2 


Ih 


* is the sign of the complex conjugate. In regions IV and V, the expressions 
for the force same att are similar. As indicated (2 the appendix, we need 
only substitute a; 1) and ay (2) by S, (1) and 8, etc. 


Let us investigate several special cases (Ref. 2, 3, 7, 8, 10). The 
formulas for the force density and the, Joule heat losses remain the same in 
every case, only the coefficients a; and a; change. 


Thus, we must substitute 5) = 55 = 0, 5 =O. in the coefficients 


a (1) and a (2), in the case of a two-sided non-symmetrical inductor without 
conductive channel walls of the MHD machine. After this, we obtain the /61 
following expressions My» G,» N,» Ly for the coefficients a (1) and a 2): 


My =e" [ay sh æ, (b— l2) —y, ch æ, (b — l2)], 
G, =0 "a, shz, (h —b) —yy ch ay (l1 —b)], 
N, =e'v" [av sha, (l — b) +y, ch a, (h —b)], 
Ly =e fay sh ay (b — l2) + yy ch a, (b—l2)]. 


In the case of a one-sided inductor without conductive channel walls of the 


MHD machine (Ref. 8), we must set 91 = 0, 2i = 0, 8) = 85 = 0, ¢6= 0 in 


order to determine the coefficients a (1) and atts. In this case, we have 


D= Bola , cha, ae 
a=K is GL—MN, ’ 


@—K Pols ch aedi G, 
ENS ee” MANGE 


where M,» G? N,» Ly are the same as in the case of a two-sided non~symmetri- 


cal inductor. 


For a two-sided symmetrical inductor without conductive channel walis 
of the MHD machine (Ref. 9, 10), we must set 


Ion =lo.=fo, p=0, h=k=l, diı=d:=d, 8,=6.=0. 


After this, we obtain 


Holo, ch a, (1—d) 


ay) —@,(2) = E00 de 
. 2r ‘aychy,bsha, (1—b) +yshy,b cha (1—6) 
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The method indicated above may be employed in other special cases to 


obtain the coefficients a) which determine the distribution of the electro- 


magnetic field, the force density, and Joule heat losses in the channel of an 
induction MHD machine. 


10. 
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TRANSVERSE EDGE EFFECT IN PLANE INDUCTION MAGNETOHYDRODYNAMIC /63 
MACHINES 


A. Ya. Vilnitis 


1. Formulation of the Problem 


In the general case, a plane induction magnetohydrodynamic machine 
(pump, generator) consists of two plane stators (inductors) with m-phase 
winding, a plane channel with molten metal, and a thermal-insulation channel. 
Figure 1 presents a sketch of a plane induction pump (1 - molten metal in the 
pump channel; 2 - plane pump channel; 3 - diffuser; 4 - thermal-insulation 
channel; 5 - pump inductor; 6 - framework of the magnetic circuit; 7 - proj- 
ection of the magnetic circuit; 8 - groove of the magnetic circuit; 9-- groov- 
ed insulation; 10 - inductor winding; 11 - front portion of winding). 


FARENE AANA 
AETEIDIOU DL 


Figure 1 


The currents in the windings of the inductors flow in the transverse 
direction, forming a traveling wave, The induced electromotive force pro- 
ceeds in the same direction,which also forms a sinusoidal wave, which travels 
along the inductor, The sign of the electromotive force changes to the op- 
posite sign (half of the wave length) along the pole division t, The lines 
of the induced currents must be closed and -- in view of the limited width 
of the pump channel -- must be forced within the limits of the pole division, 
and closed in the longitudinal direction, in which there is no supporting 
electromotive force. This fact leads to attenuation of the induced curr- /64 
ents. If the fact is taken into account that only the transverse current 
component (in interaction with the magnetic induction component which is nor- 
mal to the channel plane) produces the effective force (in the longitudinal 
direction), then the so-called transverse edge effect becomes clear: the 
smaller the ratio between the channel width and the pole division, the smaller 
the section of the closed current line over which the supporting current of 
the electromotive force has an influence, and consequently the weaker is this 
current and the smaller the effective force for one and the same magnetic 
induction, 


It must be immediately pointed out that this is an extremely simplified 
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picture. The field of the induced currents is frequently comparable with the 
primary field of the inductors, which causes the corresponding re-distribution 
of these currents. The channel thickness of the MHD machine also influences 
the current distribution over its width. As a result, it is far from simple 
to make a quantitative determination of the transverse edge effect, and it 
requires very detailed research. 


The theory of induction MHD machines finds its source in the theory of 
electric machines, and is today based to a significant extent on a system of 
concepts developed by the latter theory. Therefore, we may draw a certain 
analogy between problems of electric machine theory and problems in- 
volved in the theory of induction MHD machines. If the theory of electric 
machines is regarded as one of the branches of Maxwell electrodynamics, this 
theory may be called the electrodynamics of linear circuits with currents -- 
currents which are closely connected with a system of conductors of previously 
specified directions, which are linear for the most part -- i.e., they have 
a zero cross-section area. The entire group of specific concepts and re- 
search methods provides a good orientation for the processes occurring in the 
electric machines which are used most extensively at the present time -- 
transformers, electric engines, and generators. The theory of electric mach- 
ines -- engineering science -~ and the requirement for computational simplici- 
ty which it entails plays a primary role, frequently moving a mathematically 
precise description of the process to second place. 


If a heavy conductor or an engine containing a heavy rotor is placed 
in the groove of an electric machine, the electrodynamics of the linear cir- 
cuits may become weak, and field theory must come to its assistance. 
Since we may assume that the direction of the currents in the system is known, 
no particular difficulties arise. This is the large rotor of infinite length, 
the cylindrical inductor of the traveling field, or the plane inductor of 
infinite width, Maxwell equations may be solved by tabulated functions, and 
we obtain a comparatively favorable picture of the one-dimensional skin-ef- 
fect with respect to the radius of the rotor or the thickness of the plane/65 
channel in pure form. Nevertheless, care must be taken in this comparatively 
simple case, when trying to express the results obtained in the language of 
the theory of electric circuits, which the theory of electric machines com- 
pletely utilizes [see, for example, (Ref. 1, page 447)]. In the general case 
of conductors having finite diemensions, it is difficult to introduce the 
concept of resistance and inductance in an un-ambiguous and rational manner. 


The following factor considerably complicates the problem of an engine 
with a large rotor having finite length or of a plane inductor having finite 
width: it is necessary to find not only the density of the current in the 
rotor (or in the channel of the plane pump, respectively) but also the curr- 
ent direction, which changes from point to point. The rotor in the engine is 
a solid body. In order to draw the corresponding analogy, it is necessary 
to construct a solid plate having infinite length, but finite width,between 
the plane inductors of the pump. Such a problem may be solved by the methods 
of mathematical physics. The theory of heat conductivity, the theory of 
elasticity, the theory of wave guides, etc. extensively employ tri- 
gonometric series -- which is a specialized method of mathematical physics 
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for solving partial differential equations in rectangular regions (for example, 
in exactly the same way that Bessel functions are characteristic for cylin- 


drical regions) in these branches of mathematical physics. 


This is the form of the problem regarding the transverse edge effect, 
in which we may today speak of its exact solution. The value of the exact 
solution thus obtained should not be over-estimated, since by no means all of 
the problem of the transverse edge effect in induction MHD machines has been 
solved. A very daring simplification of the problem is based on the assump- 
tion that a solid plate is placed between the inductors, and molten metal 
is the working medium of the MHD machine. When searching for technical solu- 
tions of the entire magnetohydrodynamic problem, it is very tempting to employ 
the comparatively simple solution for the solid plate. The trigonometric 
series thus obtained are orthogonal. Therefore, if we determine the averaged 
density of the body forces 


=- 1 
i= f [iB] av, i 


f is obtained in the form of a unary series, and not a binary series, if the 
current density j and the magnetic induction B are given by unary series. It 
is obtained in the form of a binary series, and not a quadruple series, if j 
and B are given by binary series. It is also possible to substitute the /66 
initial unperturbed field B, (which would occur in the absence of a conductive 
plate) in (1) and not the total field B, if this simplifies the computations 
[see (Ref.2,p.186)]. However, maximum care must be observed. For a solid 
body (1) has limited application. For a liquid working medium, 

there is no guarantee that the pressure developed by the pump will not differ 
significantly from that computed on the basis of (1). In particular, this 
pertains to each type of deviation from the mean parameters of the machine. 


In our opinion, averaging is a very risky procedure -- integration 
over volume. The dependence 


f =[jB] (2) 


with j and B computed for a solid medium, is of value for a liquid in many 
cases, whereas the integral characteristics f does not take several important 
factors into account. Therefore, it is valid to be interested in the total 
distribution of j and B both in theory, and in structural practice. 


In every report which has been published on the transverse edge effect, 
it is assumed that the molten metal moves like a solid body. For purposes 
of clarity, we shall consider the edge effect in induction pumps, although 
the main results may be extended to the case of the generator regime, 


In more recent studies (Ref, 3-7), the effect of attenuation of inn 
duction in the clearance and the transverse edge effect were each studied 
separately, Their purpose was to obtain two correction coefficients to be 
applied to the simple technical formula for computing pressure, Upon a 
closer analysis, both effects flowed together, since no separation was appar- 
ent in the phenomena, At the same time, the longitudinal edge effect, 
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produced by the finite length of the pump channel, previously occupied a 
separate position, and was studied independently of the two preceeding 
phenomena. At the present time, we must assume that this division is reason- 
able. If a purely sinusoidal (monochromatic) traveling field is produced by 
the inductors, then the finite dimensions of the pump channel, with respect 
to height and width, or of the inductor (with respect to width) are not 
caused by dispersion of the field wave. The amplitude and phase of induction 
change from point to point, but the form of the wave remains undistorted, 

and higher harmonics are not produced. At the same time, the periodic struc- 
ture of the channel or of the inductor in the direction in which the traveling 
wave moves (with respect to length) inevitably produces the appearance of 
higher harmonics and a non-uniform non-periodic structure (inductor having a 
beginning or an end, or both) produces the continuous spectrum of the travel- 
ing waves. Due to this fact, mathematical procedure designed for studying 
the non-uniformity over the pump length differs significantly from the 167 
procedure for studying the effect produced by attenuation of the field in 

the clearance and the transverse edge effect. It would appear that it is 
advantageous first to clarify the undiscovered patterns of both groups of 
phenomena separately. Therefore, from this point on we shall assume that the 
inductors and the pump channel are infinitely long. 


Figure 2 


When examining different solutions of the problem of the transverse 
edge effect, we shal. assume that we are dealing with a system of two in- 
ductors, the distance between which equals 26. These inductors carry on 
their surfaces (which face each other) an undistorted, traveling sinusoidal 
wave of the surface current (linear current loading). We shall direct the 
axes of the cooridinates as is shown in Figure 2. We shall let the x-axis 
coincide with the direction in which the traveling wave moves (along 
the pump), and the y-axis coincides'with the direction of the surface current 
and of the real primary currents in the inductor windings (along the pump 
width). The z-axis is perpendicular to the surfaces of both inductors (along 
the channel width). We shall assume that the magnetic permeability of the 
inductors is infinite (u = œ). In the clearance, u = Ho (both in the channel 


with the molten metal, and outside of it). The width of the inductors is 
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2c (or infinity), and the channel width is 2a. We shall place the origin 
of the coordinates in the middle of the clearance between the inductors, and 
also in the middle of the inductor channels, Let the linear current loading 
of the inductors be 


A,=Agcos (wt—ax) =Re Ap eflot-axi, (3) 


: T i PRT 
Here w is the angular frequency, and a = 8 where t is the pole division, 


All of the quantities are described in a system of SI units. If we [68 
stipulate the condition that the wave phase does not change 


ot—ax=const, 


we obtain the wave propagation velocity 
m= a (4) 


It is understood that the directions A coincide for both inductors for one 


and the same x. The following boundary conditions hold on the inductor sur- 
faces (in the case of z = + 6) (no matter whether induction currents develop 
in the system or not, and independently of the degree to which the intra- 
inductor (working) clearance is filled by the current-conducting medium) ; 


for z=§ Ay,=—Abo, 
(5) 
for 2=—6 Hx=Ap. 
We shall try to find the solution in a complex form, assuming the 

presence of a phase factor ei(wt -ax) everywhere. The amplitude of this 
phase factor equals unity (i is the imaginary unit). Let us examine the com- 
plex amplitudes H, B, j and E, which we shall designate by a point above the 
appropriate letter. Then the instantaneous value of the physical quantity 
Will equal the real part of the complex quantity 


B=Re B ellat~ax) (6) 
Its amplitude will equal the modulus of the complex amplitude: 
Br=|Ba| (k=x,y,2), O) 
and the phase shift will equal the argument of the complex amplitude: 
(Bx) =arg Bn. (8) 


We shall employ the phase of the surface current of the inductors as the 
a zero phase ~--~- i.e., we shall assume that the complex amplitude 
0 in (3) is in fact a real quantity, As long as we perform. the linear oper- 


ations with the components of the fields and currents, there is no necessity 
of dividing the real parts from the imaginary parts in the complex amplitudes, 
It is another matter if it is necessary to perform non-linear operations both 
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in formulas (1) and (2). We must then keep the fact in mind that it is /69 
only possible to multiply the real physical quantities -- i.e., the real 
parts of type (6) complex quantities.If we are only interested in the magni- 
tude of the product averaged over time (period), it may be computed more 
simply if no distinction is drawn between the real parts of the complex quan- 


tities: zoi > 
i=- Re [j ó*] (9) 


[see (Ref. 2, page 243)], where the asterisk (*) designates the complex con- 
jugate. 
2. Infinitely Wide Model 


Before examining works on the transverse edge effect, let us write the 
well known solutions for infinitely wide systems. 


1. The field in the clearance between two infinitely wide and infini+ 
tly long inductors, filled by a non-conducting medium, is 


A= Hao sh az Hao sh az 


shad ichad ’ 
Fyn Heachas Hacha, ao 
Ey= <8, =H, 
Ay=E,=E,= 
According to (5) we have 
Humo Ham an 


2. The field in the clearance which is uniformly filled with a medium 
having the conductivity o is 


7 _ Ha sh Bz _ B Hosh zZ 
Ha sh Bd iachpd ’ 


Å, = 19 Hoch pz Hwch pz 


B sh pd ch pd ; (12 
. Ww wW- 
E= 738:5 H, jy=oEy, 110 


H,=jz=j:=0. 
Here we have 
B?=0?+ipo00; (13) 
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Hx=— Ao, HO Hom A. (14) 


In formulas (12) -- (14) (the same pertains to the following formulas) 
we may divide the real and imaginary part, before commencing the numerical 
computation which, however, leads tocumbersome expressions. On the other hand, 
we may perform the entire computation from beginning to end with complex num- 
bers, as is sometimes done in the theory of circuits. It is our opinion that 


the second method is preferable. 


Let us study the Capers case 6 «<T, or ad << 7. Then th aô > 0, 


and according to (11), H z0 >Å x0" But a = a 0? and Ë o play the role of 


counter-electromotive force per unit of inductor width. There is no reason to 
assume that it strives to infinity, just like the applied voltage. For a 


sufficiently narrow clearance, Å > 0andå > 0, and Å = _% Ë does not 
x y z WË Y 


strive to infinity. In addition, we have 


0H, _aHy»shaz JÈ, 
a ha O A 


but 


OH, _achaz ay 
“Oz ichas p (15) 


and not to zero. Thus, the longitudinal field component strives to zero for 

a sufficiently small clearance. This does not occur for its derivative with 
respect to z, which strives tọ the finite limit (15). In its turn, the deriv- 
ative (with respect to z) of H strives to zero; the normal field component 


becomes independent of z. In order that A + 0 may be fulfilled in the case 


of a clearance filled by a conductive medium, we must set Bi) cea, which means 


aô << q and ow << m at the same time. The latter is fulfilled in all in- 


Ko 
duction pumps employed in practice (if Ho >>a2, such a pump is not feasible 
due to an exaggerated skin-effect). If we retain the linear current loading 
of the inductors (more rA the ratio of the surface current to the /71 
clearance thickness lim Ao , since the current itself strives to zero), we 


> S 
then obtain the following when there is no conductor in the clearance 


ÒH, jie Heo. tim 0 (16) 


30 Oz Q z> ò 30 6’ 


and we obtain the following when a conductor is present 
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i.e., we obtain one and the same result in both cases. 


Formulas (10) - (14) make it possible to determine the averaged force per 
unit of surface of solid plate cross-section between the inductors, which we 
shall designate as the pressure of the induction pump: 


a d 
L-1 1 sake 
p= ares, [dugg [azU,82. (17) 
-f —ð 
We would like to add a few words to (17). We obtain p from (1) by multiplying 
f by the pump length L and averaging over time, disregarding the longitidunal 


edge effects. For whatever reason, the product ihe always equals zero. Also 
ay 


v 
fy = f = 0. We employ Bt in (17) to designate the unperturbed field (the 


field in a hollow clearance), according to the condition mentioned above which 
simplifies the computation. If there is no dependence on y, we then have 
1 a 


Za f  dy=1. The same holds true for z. 
-a 


The pressure is calculated very simply for an infinitely thin clearance 
between infinitely wide inductors if the medium filling the clearance has such 
a low conductivity o that the field of the induced currents is much weaker 
than the primary field (ug ow<<a2). We may then employ formulas (10), in which 


6 > 0, and Jy = of (although it is small, o is entirely different from zero). 
In this case (17) yields 


pom WIB aL (18) 
3. The Theory Advanced by A. I. Vol'dek (72 


Let us now turn immediately to the problem of the transverse edge effect. 
In chronological terms, the first analytical study deals with the transverse 
edge effect “in pure form" -- i.e., when there is no attenuation of the field 
in the clearance -- which is achieved in the limit by the thin clearance be- 
tween the inductors. A. I. Vol'dek conducted this study (Ref. 3). In view of 
the importance of this solution, we shall present it in its entirety below, to- 
gether with certain comments. A similar problem was investigated in (Ref. 4) 
for rotational machines. 


A. I. Vol'dek divided the field in the clearance between the inductors in- 
to a primary and secondary field. The primary field is the field in the clear- 
ance filled by a non-conductive medium. (We shall designate it by H' and E'). 
In view of the fact that 4' = 0, we must have rot H' = 0. The secondary field 
is the field of the induced currents (we shall designate it by H" and E"). 

The total field equals the sum H = H' + H", E, respectively. The Maxwell equa- 
tions 


65 


0B 
rot E=- ori (19) 


rot H=j aes 


must be fulfilled for the total field. In view of rot H' = 0, we obtain the 
following for the secondary field 


Wee 0B” 
rot E” = OF? (21) 


rot H” =j=0 (E+ E”). (22) 


The reason for this division lies in the fact that E' is assumed to be a known 
function. In coordinate notation, (21) and (22) yield 


aE” aby ari," oÈ," 


“Oy de OBS Gy Ge OE FEL"), 
ce, sei oH” |. le A 
Es +iaE,''=—iwB,”, Jz +iaHH,""=0 (E,'+E,”), 

EF," oH," 


ia," + Gt =iwB,", -ih -E =i, 


We must simplify the system (21) - (22), based on the physical conditions 
in an extremely thin plate between the inductors. It is apparent that we must 
have dy = 0 in this plate -- i.e., the current lines lie on one plane. We must 


also assume that la = 0, since the plate has the same conductivity o along the 
z-axis. For infinitely wide inductors ET = E'S = 0. This effect may be /73 


expected for an extremely thin clearance between inductors of finite thickness, 
so long as regjons are examined which are not too close to the inductor edges. 
We then have dE 'y = 0 and Ey = 0. The occurrence of EY is caused by 


ay 


the finite, wide plate. Since rot A | A, then -- according to (22) -- H" | 4, 


which 
may be satisfied in the case of H" = e H" , although this is not the only pos- 
y Ht , y 


sibility. However, preference should be given to this, based on the following 
considerations: if we retain the surface current A when introducing a con- 


ductive plate, then for Bs we obtain zero boundary conditions [see (5)]. Thus, 
tt 


we also have 2H x = 0 [see (16)]. Even when os changes, we may write 
82 
3H" 
5 x = 0, since volume currents have no influence on the boundary con- 
y zZz = +6 
ditions, and aA, = 0 always holds. Therefore, it is natural to assume that 


ay 


e 
HY = 0. We also have zero boundary conditions for Hs since A. = 0. We must 
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thus assume that A" = 0. We may express this most clearly in the following 


way. The lines of the induced current field are closed between the inductors 
along the shortest path -- i.e., We an to the plane of the clearance. 
The zero boundary conditions for H" and H" designate the perpendicular nature 


of the lines to the inductors when they enter them. However, we than have 


a Esn i =0 we have ls =0, 


and equations (21) (22) yield the system which is used as the initial system of 
A. I. Vol'dek: 


iaE", 
P a 
Hea ob 5 (23) 
laH!,=3E"', +E’, 
e 
1 
Excluding E" and m and employing 9E y = 0, we obtain the following from (23) 
dy 
Ey pEi 4 awk", 
a. = PE", +1930 (24) 
where B is defined according to (13). (24) has the solution 74 
Eps 


' =A, yeh Rust Aysh Byang E's 
and then the current is 
jy=s(E", +E')= o (Ai ch By+ Aash pu +i È' a 


The natural boundary conditions are Jy = 0 in the case of y = + a [A.I. Vol'dek 
stipulates Bs = 0 in the case of y = + a which is one and the same for the 


system (23), see (23) -- third], from which it follows that Ay = 0, since i; 


1 
must be an even function of y and A, me yY . Correspondingly, we have 
B2chBRa 

» {ch By 

impo, (1 Spe); 

on _aoE’y _ ch By); (25) 

ite ip? (1 ch Ba/’ 
iacE", sh By 
Je RC Ba" 


We feel it is necessary to add certain remarks to the solution (25). As 
has already been indicated, we may assume that Èn = Bl = oe 0 for the 
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regions (which are central with respect to the width) of the clearance between 
the inductors. In this case, solution (25) must be explained by the assumption 
that the conductive plate is narrower than the inductors. Nevertheless, A. I. 
Vol'dek extends the solution (25) both to the case a = c (plate having the same 
width as the inductors) and to the case a>c (plate wider than the inductors). 

This is based on the fact that Ws decreases sharply at the outermost edge of 


the inductor, and 2c barely ee over the entire area, so that the Biss which 
is average with respect to the inductor width, differs little from maximum Br 


e 
This latter fact is indisputable, but we would like to note that oH", >+ œo 


ot 2n h OY 
close to the inductor edges, and consequently on y >+% and Jai x > tin œ. 
' az dydz 
We also have E y a = 
ETA > + ©, from which we have io ft - > + œ., Unfortunately, 
no conclusions may be drawn regarding the quantity att >» if we assume that 
a'z dy 


is given. No conclusions may be drawn with respect to the quantity /75 


ay 
f' either. In order to do this, we must examine the solution in the entire 


region of the field propagation. 


After these brief comments, let us turn to the solution (25). According 
to (9), we have 


È =J 5 Rel, (Ë, 4 Brey), 
However, it may be seen from (25) that jy and A" differ by the factor 77 


i.e., they are shifted with respect to phase by 90°. Consequently, iB? can- 


not provide the resulting force. According to (17), the pressure is 
_ 50! B2,o| L sal is [reg ree 
Ta Re; dy= (26) 
_9@| BIL a? th Ba 
~za Rll pa / 


A. I. Vol'dek recommends that the pressure obtained (26) be expressed by 
the pressure Po (18) in the clearance filled by a medium with low conductivity, 


and introduces the attenuation coefficient kI c’ 
P= k'o Po. (27) 


This method has certain advantages for technical designs of pumps. The corres- 
ponding attenuation coefficient (a# ~, o # 0, 6 + 0) is obtained as equal to 
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pee th 
Bye =Re&(1— wee). (28) 


It should be noted that ed is obtained for constant linear current 


loading of the inductors (more accurately, for constant lim _Y „since the cur- 
620 § 
rent loading A, and the clearance thickness 26 strive to zero). In other words, 


an operation is performed of inserting a highly conductive plate of finite 
width in the extremely narrow clearance (which was hollow before this) between 
two infinitely wide inductors, maintaining the current in the inductor windings, 
and not the supply voltage, which must thus be reduced since the total winding 
resistance has decreased. 


In the case of a + » -- i.e., when the entire, infinitely wide clear- /76 
ance is filled quite well by a conductive medium (retaining lim Ay ) -- we have 
6) — 
6+ 3 
2 4 
1 t= peta = 
Rol Rae Re p? \p4|° (29) 


From our point of view, the following expression is assumed for ae » which is 


obtained by expansion of th fa in elementary fractions: 


p2 N 1 1 
a >; 2 (a? + %,2) pp BPU?’ (30) 
where i (+ na)? 
tn = PRR (31) 


Under the condition that Ugo << a2 (a#m,o > 0, ô > 0), we obtain the fol- 


lowing from (28) ‘ 
thea 


aa 


ka"! =1 -7 ’ (32) 


which may be readily determined according to a table of functions. The dif- 
ference a may be written in the form of a rapidly converging series: 


k ik 1 2o’ 0? y EL EEEE E 
‘oc ‘oc a? TEETE pow? (33) 
n=l (a? + %,2)? 


In the case when a conductive plate thickness does not occupy 
the entire clearance 25, but only part of it, 2A, so that 6 > A, A. I. Vol'dek 
adheres to the previous hypothesis xs 0. In this case, the magnetic resist- 


ance to the flux produced by the induced currents increases by a factor of 6/A, 
which is equivalent to a decrease in the magnetic permeability (or conductivity) 
of the medium by as much as ("o > "25) ` 


6 
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It is a very simple process to pass from a plate at rest to a moving plate: 
it is sufficient to multiply the cyclic frequency of the field w by the slipping 
s (w > ws, wheres = 1 - v4 oe velocity of the plate; Ve -- velocity of the 

Yg 
field) [see (4)]. 


As was already pointed out, koc! may be derived for a constant linear cur- 


rent loading of the inductors. If we have a plate with high conductivity or if 
we are operating at a high frequency, then ug becomes larger, or even /77 


much larger than, a2. In these cases, kaas has very small numerical values. 


According to the terminology employed, we thus have a strong armature reaction 
(which is expressed as a drop in the total winding resistance). If it is not 
balanced -- i.e., if the previous linear current loading of the inductors is 
maintained -- then the effectiveness of the pump greatly decreases. However, in 
practice (if only the temperature regime of the winding permits it) it is fre- 
quently not the current which remains constant, but the supply voltage, and 

the armature reaction is balanced. It is difficult to make a precise calcula- 
tion of the effect of a constant supply voltage, due to the absence of a single 
theory for a non-magnetic clearance and for the inductor itself: the real in- 
ductor is replaced by a ferromagnetic body with a surface current loading for 
which there is no active inner resistance, nor inductive inner resistance. If 
we retain this idealization, we may only speak of counter-electromotive force 
to the supply voltage, caused by the variable magnetic field in the clearance. 
The constant supply voltage is then equivalent to a constant magnetic flux a 
for the pole division. In infinitely wide systems, a constant %, is equivalent 


to a constant Bo In systems of finite width, there is no point in speaking 
= fr pit - ou d z 
of a constant Bo B 20 + B z0? Since B 20 epends on y [see (25)]. According 


to this, we may introduce the attenuation coefficient for a constant supply 
voltage (in the sense of a constant magnetic flux 3), i.e., with compensation 


for the armature reaction. As has already been pointed out, the theory of A. I. 
Vol'dek encompasses the tacit assumption that the inductors are wider than the 
channel (to a sufficient extent). Let us clarify this assumption: c > a. It 
is assumed that the winding is of a diametrical step. We then have 


wy = iw®,° — EF ) 
T a 
U= iw®, = foe Ba'+ iw J oe dx J È 0 dy= 
0 a 


_ 40 | fend 2 bey, 
= Fg! [e+ et SE); 


pa 
e 
and stipulating that U = 00, we obtain , 
> Ba 
Bay’ = 
z0 1 4 Hoge 2(i- th Ba) 
ip c pa 
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and /78 


a? (1— t Ba) 
ka= Re - a) eae ; 
P- ipow t (1 — th in ba) (34) 


It may be readily seen that lim bo = 1 -- i.e., if the supply voltage 
a=c7" 
is maintained, the armature reaction in an infinitely wide plate is completely 
balanced. 


4. Other Theories with Plane-Parallel Field in the Clearance 
of the Machine 


L. V. Nitsetskiy (Ref. 5) has formulated a principle for drawing up a model 
of certain vortex fields in an electrolytic bath and on electroconductive paper, 
In particular, this principle may be applied to formulate a model of the trans- 
verse edge effect in an extremely thin plate. It has been verified by the 
author by comparison with the results derived from performing the calculation 
according to formula (32). The agreement was satisfactory. 


L. Ya. Ulmanis (Ref. 6) employed an analytical method to obtain the fol- 
lowing attenuation coefficient for an extremely thin finite width of the plate 


hae’ = lB Re& (1 - th Ba) (35) 


([Ref. 6], formula (1); the previous notation is employed). ca is determined 


as the ratio of the pressure in a plate of finite thickness to the pressure in 
an infinitely wide plate, with allowance for the armature reaction: 


p=k,."p'o, 


where P'o = A is the pressure in an infinitely wide, highly conductive 


oc 
plate, which is chosen as the initial pressure, instead of Po (18). As the 


plate narrows from an infinite width to the width 2a, a constant current load- 
ing of the inductors (sa ) is maintained. L, Ya. Ulmanis has checked (28), 
60 — 
§ 


(32) and (35) experimentally, and has provided the following experimental form- 
ula 
pisan on 
oc z 2° 
1+13(5) (36) 
Nothing more similar may be obtained, assuming o = 0 in the expansion of [79 
(30) and taking the first term of the series 
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Figure 3 Figure 4 


In an article by A. I. Vol'dek and Kh. A. Yanes (Ref. 7), the theory which 
was discussed for an extremely thin plate between inductors was generalized to 
the case of a plate with lateral short-circuiting bus-bars (Figure 3). Zone I 
is the plate which simulates the channel of a pump with molten metal (it is im- 
mobile here!), and the conductivity of the zone is o. Zone II is the short- 
circuiting bus-bars with the conductivity o3. The condition 6 > A is taken in- 
to account by the introduction of equivalent conductivity of the zone o' = o A, 

6 


which is similar for Jos It is assumed that the primary magnetic field is fully 


concentrated in the clearance: in the case of |y|<c Bt = const, and in the case 


>c P: = 0 (just like È). Two special cases were examined in this study. 


of |y 
1. a= c (short-circuiting bus-bars outside of the inductors). Equation 
(23) holds for zone I. BS = 0 in zone II, and the field of the induced cur- 


rents is also disregarded, in view of the fact that magnetic lines pass along a 
non-magnetic medium. Therefore, we may formally substitute Ho = O in (23). 


e 
The boundary conditions are: in the case of y=+t HS = 0, and in the case 


= P 34 15 rt 
ce ae Jy|z ~ Jy|I1; È x|1 z «IT. In addition, it is assumed /80 
that cae is an even function of y. The following result is obtained: 
A," eo . 
zi p? ao ch Bc+ oof tha(f—c)sh Bc }’ 
Jyh=iatt,” h; (37) 
k vizre% {1 — aoshe a 
oe p? Be lao ch Be + Boo tha(é—c) sh pc)l’ 


Substituting Jo = 0 from (37), we may obtain the formulas for a channel which 


is wider than the inductors, but without short-circuiting bus-bars. 
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2. t = c (short-circuiting bus-bars between the inductors). In this case, 
we have equations (23) for both regions; the boundary conditions are the same as 
before. The result is 


_. _ [(oB%2—o2f?) ch B2(c—a) + 0287] ch py } 
P2[02ß sh Bo(c—a) sh Ba+ofech B2(c—a) ch Bal)’ 


jy =lat" l, (38) 


2 
Rax! = Re F fı = 


= sh pa [(oB*2— 028°)ch B2 (c—a) +028% | 
Bae [02B sh B2(c—a) sh pa +oß2ch Be (c—a)ch Bal)" 


Here B57 = a2 + in goo. Cases are possible in which kator for small values 


of o and Sy » and also a. 


VII rT. 
koc and P PI are introduced under the assumption that the linear 
current loading of the inductors is constant when conductive plates are intro- 


duced into the clearance of the system. 


5. Transverse Edge Effect in a Clearance of Finite Thickness 


The following stage in an investigation of the transverse edge effect in 
induction pumps is related to the transition from an extremely thin plate be- 
tween the inductors to a plate of finite thickness. It is discussed in the 
publications of A. Ya. Vilnitis (Ref. 8), T, A. Veske (Ref. 9), and N. M, 
Okhremenko (Ref. 10, 11). In veiw of the fact that identical solutions are 
obtained in the studies (Ref. 8, 9), in spite of certain variations in the formu- 
lations of the problems, it is advantageous to analyze them concurrently, 


The following problem is postulated in the study (Ref. 8). There are /81 
two infinitely wide and long inductors with u = œ, which carry a linear current 
loading in the form of a sinusoidal traveling wave of the surface current A, 


(Figure 4). A conductive body having rectangular cross-section and the thick- 
ness 26 (which equals the distance between the inductors) and the width 2a is 
placed in the intra-inductor space. In the remaining intra-inductor space 

o = 0. It is assumed that the body is immobile, but -- as was already pointed 
out -- it is sufficient to multiply the cyclical frequency of the field w by 
the slipping s for the transition to the mobile body. When the problem is 
solved, the boundary conditions jy = 0 in the case of y = + a and Í, = 0 in 


the case of z = + 6 become apparent, in addition to the boundary condition (5). 
It is also necessary that the exciting action of the body may be reduced to 
zero at an infinite distance from this body, which is placed between the in- 
ductors -- i.e., in the case of y > +œ , the field between the inductors may 
be expressed according to (10). Finally, it is necessary that the magnetic 
field be constant in the case of y =+ a. 
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The solution of the problem which is given in (Ref. 8) follows very clearly 
from these conditions. Two qualitative consequences are clearly apparent in 
this problem, which were not apparent from the very beginning: 


(1) There is no transverse component of the magnetic field 2 over the 


entire clearance, both outside of the conductive body, and inside of it; 


(2) The field outside of the conductive body is an unperturbed field not 
only in the case of y > + œ, but also for any |y|/> a. 


By analogy with the formulation of the problem given by A. I. Vol'dek for 
an extremely thin plate, the author (Ref. 9) first assumed that the plate was 
thinner than the clearance between the inductors. However, when written out 
the solution everywhere becomes 6 = A -- i.e., the solution may be written for 
a plate occupying the entire width of the clearance. Secondly, the author /82 
did not define the problem of the inductor width more accurately, and only 
assumed that the primary field is described by expressions (10). This is 
equivalent to the assumption of infinitely wide inductors. Thirdly, it was 
postulated that the secondary field (the field of the induced currents) does 
not pass beyond the side boundaries of the conductive plate . It is not clearly 
apparent that this is an initial condition, but this is actually the case. 


If the solutions given in (Ref. 8, 9) are reduced to one and the same 
notations, the fact that they are completely identical becomes readily apparent 
-- i.e., expressions (Ref. 8) (42) are identical to (Ref. 9) (4), and expres- 
sions (Ref. 8) (41) are identical to (Ref. 9) (5) (with allowance for the fact 
that the total field is given in [Ref. 8], and that the primary field (Ref. 9) 
(3) must be added to (Ref. 9) (4), (5). This fact apparently indicates that the 
problem may be solved without errors under the given premises, and that with 
these solutions the next stage in terms of complexity is, in a certain sense, 
exhausted after the solution for the extremely thin plate, given by A. I. Vol'- 
dek. We cannot substantiate the fact that all of the qualitative patterns hid- 
den in the above-mentioned solutions have been clarified. Their clarification 
is of essential importance for orientation when proceeding to the next stage of 
the problem in terms of complexity. 


The expressions for the fields and the currents are given in six different 
representations in (Ref. 8) and (Ref. 9). We shall give the expressions for 
B_ and jy (i.e., for the components which are of interest when computing the 
pressure developed by the pump) in a form which is the most advantageous for 
the computations: 


> Jia chgz_ pooo, ch py 
B, = Baa f. shpd dp? ch pa” 


kee] 
2 pooma >; (—1)"™7!chymy 
FO Ly Vay m chyna “OS hm? } a?) 


40 
eB fe, ch gz _isw, ch py S 
k p shpd ap chfa 
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= t 
aey (—1)7-! ch Ymy ios dine! 
V'mchyma j (40) 
where 
jr Eu, = B24} 2-62 
Àm 5 ’ BHn = B+ Be ; (41) 
Ym? =O? 2,7, (42) 


and B is given by the expression (13). The series in (39) and (40), which can- 
not be summed up in tabulated functions, converge quite rapidly, pri- /83 
marily due to the factor Chymy , since ly |<. In the case of y = + a, the 


hyna 
series are summed up as follows: 


D a) 


(=1)"-1c0simz_ [fy _adchaz 
2 Vm Ym Dit ow oe EE shad, (43) 
1 — Boch Bz\) 
-p ~shpo jj’ 


It may be readily seen that (39) and (40) are valid both over the region, and 
at the boundaries of the conductive body. 


It is interesting to note that the first terms in (39) and (40) coincide 
with the formulas (2) for an infinitely wide channel, and the second terms may 
be written as follows 


ch BY limb —limé,, 


-Bao að? “chBa 2.07 330 
a>%  a@—0 

040 040 
oe, ch BY lim j — lim 
— Bxo ap? ch Ba ny z jy 


ao Žž a =e 


Thus, (39) and (40) may be written as the superposition of the special solutions 
(which are previously known) and of the special solution in the form of a series, 
which cannot be reduced to tabulated functions and is not amenable to 
separation of variables: 


B,=lim B: + lim B, - lim B,+ Fs (y, 2), 
0=40 3-30 o+0 
aco = ato aco (44) 
040 040 o40 
jy= lim jy + lim fy — lim jy +F; (y, 2). 
L 30 Iv Sardi 3—0 d : 
a= azz0o a-—o0o 


(44) is not the only possibility for separating the simpler special solutions 
from the general solution. Other examples are presented in (Ref. 8, 9). In 
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this connection, there is a good possibility of employing simpler special solu- 
tions of a different type, instead of the total solution, in the appropriate 
cases. However, for this purpose the magnitude of the discarded portion of the 
total solution must be determined. A general theoretical analysis of this /84 
problem has not yet been conducted. A calculation of the special case (a = 0.13 
m, ô = 0.015 m, t = 0.45 mn, liquid potassium at 500°C) on a digital computer 


has shown that the terms F and F; in (44) can be discarded in this case up to 


y = + 0.12 m (in terms of the modulus, they are no less than three orders smaller 
than the sum of the first three terms) and only in the immediate vicinity of the 
edge (y = + 0.13 m) do they become of the same order of magnitude as the remain- 
ing terms. 


According to (17) and (27), the corresponding attenuation coefficient 
equals 

fps aie ia? cth Bô th pa 

RIX, = th? aô Re Too cop - pad? 


22 SV thyma 
aò? Ym Vm 


m=l 


(45) 


It may be readily seen that, in the case of ô > 0, the kot of A. I. Vol'dek (28) 


is thus obtained. When verifying this, we must take the fact into account that 
th BS > BS 

_thymd o iat at 
330 Ô? Ym y2m and ooog Pi Regs 


Kose is introduced into (45) in exactly the same way as A. I. Vol'dek intro- 
duces his koc! (28): a conductive plate having a rectangular cross-section is 


introduced into the empty clearance between two inductors. The linear current 
loading of the inductors does not change. Since the clearance between the in- 
ductors has finite thickness 25, when a change is made from an extremely narrow 
clearance (p = Po ko Po is computed for a clearance which is infinitely 


narrow) to the clearance 26, the induction B'o remains constant just as before 
(the inductors are infinitely wide, and B' 0 does not change over their width). 


Just as previously, we may introduce the attenuation coefficient in the case of 
a constant supply voltage (in the case of a constant magnetic flux $). This 


yields 
2 
pX eRe = (46) 
S apo oo Y 
whe 
Yei(c~a) cth ad+ č cth pò — ne TY th pga — 
B ad p? (47) 


_2poooa N thyma, 
ô 7 ym Y 2m i 
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(34) is thus obtained with the limiting transition 6 > 0. 


85 
TABLE 1 
Nhe le pe Te 
e 0, 0.3 0,5 1,0 4,0 
D 

0 0,0317 0,2185 0,4151 0,6787 0,9120 
0,02 0316 2180 4140 6769 9096 
0.1 0,10 0298 2051 3894 6366 8554 
’ 0,25 0240 1563 2939 4781 6410 
0,50 0158 0917 1638 2577 3403 
1,00 0084 0462 0785 1471 1501 
0 0,0317 0,2150 0,3923 0,5928 0,7483 
0,02 0316 2144 3912 5912 7464 
05 0,10 0298 2018 3680 5560 7019 
" 0,25 0240 1539 2780 4181 5269 
0,50 0158 0904 1559 2282 2843 
1,00 0084 0456 0754 1068 13tt 
0 0,0316 0,1895 0,2700 0,2951 0,3045 
0,02 0315 1890 2693 2943 3037 
L5 0,10 0297 1778 2534 2769 2857 
: 0,25 0239 1360 1927 2106 2174 
0,50 0157 0812 1132 1253 1311 
1,00 0084 0416 0581 0667 0721 
0 0,0305 0,0820 0,0649 0,0515 0,0417 
0,02 0305 0818 0647 0513 0416 
50 0,10 0288 0771 0610 0484 0393 
>” | 0,25 0232 0609 0495 0404 0339 
0,50 0153 0408 0366 0331 0306 
1,00 0082 0221 0206 0191 0181 
0 0,0276 0,0305 0,0218 0,0158 0,0114 
0,02 0276 0304 0217 0158 0114 
10.0 0,10 0261 0288 0206 0150 0109 
: 0,25 0211 0245 0187 0148 0119 
0,50 0140 0183 0153 0134 0118 
1,00 0075 0101 0087 0077 0069 


Turning to (45), we should point out that ie may bè also represented in 


the following more symmetrical form: 


2 œ 
am TEE |C? È odn), (48) 
m=l 
where 
at thy„a\. 
(eda Real a ; 
af (18a) gs, G2 
(ku'o = Re Figs l= a =k; 


(koc D0 is obtained by substituting m = 0 in (k ePm; ar ce are deter- /86 
mined by the expressions (41) and (42). 
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The values of eae for several parameters are shown in Table 1. The val- 


ues of koc in the case of D = 8s 0 correspond to the attenuation coefficient of 


A. I. Vol'dek Kes As we may readily see, 


provides a good approximation up to D = ô = 0.02 (the error does not exceed 
T 


0.25%). However, even in the case of D 0.10 only the first sign is valid. 


The simplified dependence 


thaa 
a 
se l aa 
is valid within very narrow limits (in the case of A = Ž to 0.1, D -~ up to 
0.02, and e = HOLO == up to 1.5). This is clearly inadequate for practical applic- 
2 
a 


ations, since £, and particularily A, fall outside of these limits. 


We may obtain much better results by a comparatively simple method -- 
multiplying k__I by th2a6: 
oc a262 

k Xi thg Re& (1 th Ba 
oc a2? p ~ Se). 


In this case, the error does not exceed 1% up to D = 0.25 in the case of e < 1.5 
and A > 0.3, and in the case of D < 0.10 this holds for any A up to s = 10.0. 


Figures 5 and 6 present curves of B Jy and A which were calculated on 


a computer according to formulas (39) and (40). Figure 5 presents the distribu- 
tion of these quantities with respect to y in the middle of the channel (in the 
case of z = 0) for the special case mentioned above (á = 0.13 m, 6 = 0.015 m, 
t= 0.45 m. The entire clearance having a thickness of 0.03 m, is filled with 
potassium at 500°C). Curve 1 gives the amplitude B5 2 - phase; 3 and 4 - the 


same for j_; 5 -~ the mean force density (effective pressure). Figure 6 presents 


the distribution of these quantities in the case when only part of the clearance 
is filled with liquid potassium (a = 0.14 m, 26 = 0.024 m 2A = 0.014 m). In 
addition, transverse partitions made of stainless steel lower the effective con- 
ductivity of the body (Occ = 0.95 » 106 1/ohm - m). A re-calculation was /87 


performed on the recommendation of A. I. Vol'dek:e > e zs 


It may be seen in Figure 5 that, in the case of a highly conductive medium, 


the pressure PL between the inductors (e = 8.64) changes greatly over the 


channel width. It is comparatively simple to combat the decrease in ia down to 


zero along the channel edges, since this decrease occurs in a comparatively 
marrow zone (0.01 m ). However, the strongly expressed transverse skin -- 
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Figure 6 


effect -- the decrease in the pressure by almost a factor of two in the wide 
zone in the center of the channel -- requires special compensatory measures. 


The studies of N. M. Okhremenko (Ref. 10, 11) postulate a more complex 
problem, with the purpose of providing a better approximation of a real, plane 
pump of finite width: two inductors with a sinusoidal surface current load-/88 
ing, which produces a traveling wave; a non-conductive layer of thermal insul- 
ation on the surfaces of both inductors; a more conductive, immobile layer (walls 
of the pump channel) above, below, and in the middle -- a conductive layer /89 
which moves like a solid body with a velocity of v = (1 - s)v. in the direction 


of displacement of the traveling field wave (Figure 7). The average layer re- 
sembles a molten metal in the channel of the pump. 


If we check carefully the boundary conditions employed by the author, 

we shall see that they all apply only on different planes which are perpendicular 
to the z-axis (parallel to the inductor surfaces). On the flat sides (in the 
case of y = + a), no boundary conditions hold. Thus, it would appear that the 
solution obtained does not depend on the conditions beyond the flat sides of 

the channel and the inductor. This is an illusion. Based on the general form 

of the solutions of N. M. Okhremenko, we may readily establish these conditions 
which hold on the lateral surfaces of the channel y = + a, and beyond this, and 


79 


Figure 7 


we may formulate the model for which there is a solution of this type. 


We would now like to direct our attention to the fact that 8 ov ois: 0, 
since 8 = auth mt (the notation provided by N. M. Okhremenko). PEEN 


in the case of y = + a we have A = 0, A, = 0, Jy = 0 in every region (see 


sections 5-7, [Ref. 10]). These are the same boundary conditions which would 
hold if we used surfaces of ferromagnetic media without current loading in the 
case of y = + a. Thus, the first possible model which the solution of N. M. 
Okhremenko describes is as follows. There is a rectangular cavity in a ferro- 
magnetic inductor, in which a single-layered system is located. There is cur- 
rent loading in the form of a traveling wave on the horizontal surfaces of the 
cavity. There is no current loading on the vertical surfaces. A transverse 
field component ey must naturally arise in this model. 


However, this is not the only possible model. Such a solution is also 
obtained if the inductors, and also the channel, are infinitely wide, but the 
current loading over the width of the inductor is distributed according to the 


following series 
Qo 
„= Ae D4 ; cos Cyt I) (2v+I) ny 
T hol ~ 2a 


The sum of this series equals a in the case of -a<y<a, 3a<y<5a, 7a<y<9a, etc., 


and ZRN in the case of a<y<3a, 5a<y<7a, etc. The current loading on the in- 
ductor surfaces changes its sign over the width of the inductors with the period 
2a. An infinite succession of inductors arranged side by side is obtained. 

The opinion advanced by N.M. Okhremenko, to the effect that the sharp peaks 
(observed during the experiment) of the component 8 along the inductor edges 
substantiate the validity of the expressions for B and B, ([Ref. 10], page 23) 


as well as the similarity between the model under consideration and a real /90 
pump, is not convincincing. 
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In the solution of N. M. Okhremenko, a vertical current component (È = 

is lacking. z 
= 0) / For this reason, we may make the following statement: the fact that the 
obvious boundary condition j; = 0 is not employed in the case of z = A/2 means 


s 


simply that it is assumed that is 0O over the entire cross-section, as was the 
case from the very beginning. Otherwise, we would have to set Íz in the form of 


a series, which yields zeros in the sum in the case of z = A/2, but which may 
yield a non-zero sum for other z. Only if it may be shown that the, coefficients 
of the series Ï, equal zero, must we have 2 0. The presence of $y (due to 


ferromagnetic iei surfaces) makes it possible to assume that i, = 0. In 
this case, there are no apparent contradictions. 


Finally, attention must be called to the fact that the boundary conditions 
(section 2 [Ref. 10]) 


j! j" j" im 


are invalid. Let us try to determine what the boundary conditions must be for 
the electric field on a plane boundary of two conductive regions which slide 
over each other and which have complete electrical contact (when changing from 
one region to another, the current does not encounter any resistance). Region 
II is immobile (Figure 8); region III slides along it at a velocity of v = es 


= ey nit ~s) = Wee! - s). The boundary of the regions tesa through the 


ere z= b/2. ee the electric field in region II equal git in the case 

= b/2. The field in region III will be E*III = EIII + [vB], where EIT! is 
the field measured in an immobile coordinate system (connected with region ITI). 
E*III is the field in a system which moves together with region III (effective 
field; jIII = oIIIgsIII), This is the manner in which the Maxwell equations 
and the data of N. M. Okhremenko are interpreted (section 1 [Ref. 10]). 


Figure 8 Figure 9 


Unfortunately, B does not undergo a discontinuity on the boundary. /91 
We thus have 


Et = EM, 
E,“ n =E" — vB, 
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LW = EM 4 By, 


The first Maxwell equation is 


i 
iann — rot EM; 
m 
a =—rot E!!! = — rot E*! + tot [vB]. 


Let us draw the contour, as is shown in Figure 9, where h «2%, hg = Z, and 
let us multiply the left and right parts of the equations by a surface element 
ds. Let us then integrate over every È. We apparently have 


oB 
lim / — ds = 
lim | at 2550, 
= 
Since B is everywhere finite; 


i 
lim frotEds=tim fp Edl= f(E" —E") dl. 
Bats h>0 Pi Z 


In addition, n J rot[vB]ds depends on the mutual orientation of v and ds, 
h> 
z 


and equals zero in the case of v | ds. However, we do not even require this, 
since we have the dependence of E*III on EIII and B for the special case. We 
may then draw the conclusion that, independently of the orientation of v, we 
have the following on the surface: EHM = E IH or, in the given case, in the 


case of v = ev: 
E,*W= E11 = Eu, 


E, "=E, — vB; =Ey"'—0,B,. 


Changing to currents, we obtain 


fe" oJ j" J vB 


oll gil? “gill” “gil xoz (50) 


instead of section 2 (Ref, 10). The lack of agreement in the first of them/92 
is apparent. With respect to the second, we have 
1 OEM 


È L =! +1 dE as 
B, To POLE "CAR ETENE oy ; 


Ee Et — 2 (1—s) (C È+ 1 fat ) 


iw’ dy 
: i(i—s) 0E," 
msE,t p08) Se, 


And as a result 
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J" osd iUs), OE." 


“gil ~? Git PaaS tea 


gl a dy ’ 
which coincides with the boundary conditions of N. M. Okhremenko when, and only 
when, 


This occurs only in infinitely wide channels. 
We may draw the following conclusions from the above statements, 


At the present time, a sufficiently comprehensive solution has been pro- 
vided for the problem of the distribution of fields and currents in a rectangu- 
lar plate between two infinitely wide and long inductors which occupies the 
entire width of the clearance. The corresponding attenuation co- 
efficients of the pump efficiency are obtained (correction factors to the ele- 
mentary technical formula for computing the pressure Po = owl |B) 42 | ) The 


20 
formulas obtained make it possible to compute the desired quantities readily 


on a computer. In every case, terms with series of complex hyperbolic functions 
may be discarded (the sum of the series practically never exceeds the duplicate 
first term of the series). Then all of the requisite quantities may be ex- 
pressed in tabulated functions, which may be calculated by hand without any 
special effort. 


The results obtained find limited application in the construction of in- 
duction pumps, due to the great idealization of a real pump (solid conductive 
medium fills the entire clearance, inductors are wider than the medium and are 
replaced by ferromagnetic half-spaces with sinusoidal surface current loading). 
These results also find limited application due to the lack of rigorous 193 
solutions which provide an estimate of the error involved in thís idealization. 


One of the problems of paramount importance for future research is the 
development of a theory for a channel of finite width, whose thickness is less 
than the clearance between the inductors. In these solutions, the absence of 
perturbation of the primary field outside of the conductive plate produces a 
certain effect. This effect must vanish as soon as the plate thickness be- 
comes less than the clearance between the inductors, since the lines of the 
field will tend to pass around a highly conductive medium. In the limiting 
case £ + © complete displacement of the field from the conductive plate must 
occur. The caa&aculations must be conducted on the basis of the assumption of 
infinitely wide inductors. The physical meaning of the lack of latteral dis- 
tortion of the force lines must thus become apparent in the cases under con- 
sideration. 


The lack of a single theory for an induction pump (pump channel + in- 


ductor) reduces the accuracy with which we may calculate such quantities (which 
are so important in practice) as the active and inductive resistance of the 
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pump winding and its active and reactive power. It also reduces the accuracy 
with which we may calculate the changes in these parameters which are caused by 
changes in the secondary circuit of the pump. This problem is closely related 
to the problem of errors entailed when a real, tri-phase winding of the in- 
ductor is replaced by a system of an infinite number of phases (the higher harmonics 
are disregarded on the assumption of a sinusoidal form of the current loading). 


It is necessary to find a solution for inductors having a finite width, 
and it is desirable to take into account the influence of frontal sections. 
Narrow peaks of the transverse induction component have been observed experi- 
mentally along the inductor edges. There are theories which either ignore this 
phenomenon (Ref. 3, 7-9) or else they explain it incorrectly (Ref. 10, 11). 

The presence of B_ may lead to the re-distribution of By and Ez? which quickly 


has an influence on the pump efficiency. 


Finally, it is necessary to find an approximate hydrodynamic solution for 
a liquid-metal channel which is divided into several parts by longitudnal par- 
titions (this is a case which extremely important in practice). This makes it 
possible to correctly select structural methods for combating an excessive 
change in pressure along the channel width and to combat the danger of inverse 
flows of the molten metal, which greatly reduce the pump efficiency. 


The main emphasis must thus be placed on a clear representation of the 
concepts to be used, without struggling with complex computational formulas, 
The clear representation of the concepts is absolutely requisite for orienting 
pump construction in practice (thus, we do not advise an approximate machine 
solution of the existing differential equations for the given boundary con- 
ditions). The computational formulas contain, as a maximum, rapidly-con~ /94 
verging series of hyperbolic functions of a complex argument, These functions 
are apparently fundamental solutions in an investigation of induction 
processes in rectangular regions and may be readily calculated on computers. 
When one investigates semi~infinite regions (inductor of finite width), it is 
true that we may expect a solution in the form of Fourier integrals, which may 
lead to certain difficulties in the numerical calculations, 
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LONGITUDINAL EDGE EFFECT IN LINEAR INDUCTION MHD MACHINES /95 


Ya. Ya. Valdmanis 


1. Introduction 


As is well known, the development of a theory for linear induction MHD 


machines was initiated with research on the simplest idealized models, in which 
it was assumed that the dimensions of the inductor (primary circuit) and the 
molten metal (secondary circuit) were infinite in the longitudinal and trans- 
verse directions. Refinements related to the finiteness of real MHD machines, 
which were independent of each other, were introduced in the successive stage. 
The group of phenomena related to the finiteness of the dimensions in the 
longitudinal and transverse directions was designated as the longitudinal and 
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transverse edge effect, respectively. These effects were studied repeatedly 
both theoretically and experimentally. This article is devoted to an investi- 
gation of the present state of research on one of these problems -- the long- 
itudinal edge effect. We would like to point out that a brief summary of this 
problem may be found in the work by A. I. Vol'dek (Ref. 1). 


Before the development of linear MHD machines, the longitudinal edge effect 
was studied in connection with the development of asynchronous engines with 
arc and linear stators, but we shall investigate this problem from the view- 
point of MHD machines with allowance for their specific properties (unlimited 
secondary curcuit and a practically infinite magnetic permeability of the mag- 
netic circuit). Assuming that the channel of molten metal is infinite, from 
this point on we shall only relate the longitudinal effect to the finiteness of 
the inductor, which is represented by a smooth magnetic circuit (in theoretical 
computations) with linear current loading given on its surface, in the form of 
a traveling wave. 


The longitudinal effect of a finite inductor without a secondary circuit 
is manifested in the presence of additional pulsating fields, apart from the 
traveling field, in the working clearance. In real inductors of MHD machines, 
these pulsating fields are propagated over the entire length of the inductor 
with almost a constant amplitude. We should recall that, in a similar in- 
ductor of infinite length in the longitudinal direction, the field in the /96 
clearance has the form of an undistorted traveling wave. The formation of 
pulsating fields in the clearance is an undesirable phenomenon (the symmetry of 
the currents is distorted, the losses increase, etc.). Therefore, different 
methods have been advanced for inhibiting this field. When there is a secondary 
circuit, there is an unusual flow of the induced currents in the molten metal, 
beyond the limits of the inductor active zone, and the field is carried along 
in the direction of motion. 


Let us first investigate the field of a finite inductor (longitudinal edge 
effect in a primary circuit), and let us discuss changes in the field distrib- 
ution when there is a secondary circuit (longitudinal edge effect in a secondary 
circuit). Finally, we shall point out directions which may be pursued by further 


research, 


Let us first turn to a specific examination of the problem, and we would 
like to say a few words regarding the general computational method which is 
applied in almost every study devoted to the longitudinal edge effect, 


Only the electrodynamic portion of the computation is investigated in 
every report -~ i.e., the molten metal of the secondary circuit is replaced by 
a solid metal moving at a constant velocity. All the results are obtained 
from a solution of the Maxwell equations (in differential or integral forms) 
with the corresponding boundary conditions. 


2, Magnetic Field o; of a Finite Inductor 


ee CN eee eee 


The field of a finite inductor was first studied most extensively in a 
work by G, I, Shturman (Ref. 2), but the possibilities inherent in this study 
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found practical utilization by A. I.Vol'dek and co-workers (Ref. 3-5, 17). 


Let us discuss in greater detail the results derived in these studies. The 
following problem was studied theoretically in (Ref. 2) (see Figure 1). 


Figure 1 


We have a magnetic circuit with the height h and magnetic permeability uy. 
The current loading is given in the form of -J, cos (wt ~- az) [in the study 
(Ref. 2) Jo = A] on the active portion of its “length 2pt. In order to /97 


make allowance for the shunting currents at the ends of the magnetic circuit, 
portions of the length Y are left unwound, The lower magnetic circuit, which 
does not have winding, has the same height h, The air gap equals ô. 


Assuming that the magnetic induction in the clearance has only the com- 
ponent B, (only B in the magnetic circuit) and that the magnetic current at 


the ends of the magnetic circuit (in the case of z = + (Y + pt) equals zero, 
we obtain the following for the induction in the clearance in the region of the 


wound section: 
B=Bmsin (wf—az) — (— 1)’ Bm [ke ch pz sin wf — 1) 
—k, sh Bz cos ol], 


where 


B= ofo_. p ShBY . ,_  pehpyY ., (2) 
™ 8(1+B%/a?)’ “Ssh B(¥+pt)’ ~*~ ach B(Y+pr)’ 


oe 2p . aow 
B= Vins 5T (3) 


As may be seen from an analysis of expression (1), the first term re- 
presents the traveling field, and the last two terms represent, respectively, 
the portion of the pulsating field which is symmetrical over the inductor 
length (ach 8z) and non~symmetrical (ash gz). In the case of p which is odd 
(even), the phase of the symmetrical field coincides (shifted by 180°) with the 
phase of the traveling field in the middle of the inductor (z = 0), and the 
phase of the non-symmetrical field is shifted up to 90°, as compared with the 
symmetrical field. If the magnetic permeability of the magnetic circuit strives 
to infinity [B + 0 corresponds to this (3)], the non-symmetrical portion of the 
field disappears, and the symmetrical component is constant over the inductor 
length. Consequently, the non-symmetrical component of the pulsating current 
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is caused by the finiteness u of the magnetic circuit (saturation). Since the 
case mentioned above is observed in linear MHD machines, let us write expression 
(1) in the case of 8 > 0 


B=B,,[sin (ot—az) — (— 1)? k: sin ot], (4) 
Y 
k= Tape" (5) 


For an experimental verification of (4), let us calculate the distribution of 
the effective induction over the inductor length 


ü tjo B ta 
e V a f B'dt= 72 I= (CIP kos az Fk, (6) 


from which we readily obtain the maximum value Boff = Bia? and the mini- /98 
mum value B =B, k and B _ : 
eff min, c m 
B B 
Ba = 5 (1 +k.) Bain “Va (1 —k,), (7) 
Bmax — Beni B Bmart Bi 
Ro —max_ ming Sm _ Tmax min 
e 2 , y2 PA hd (8) 


Figure 2 presents the distribution of the field, which was measured experi- 
mentally, in the inductor clearance of one of the pumps constructed at the 
Institute of Physics of the Latvian SSR Academy of Sciences [similar curves are 
also given in (Ref. 5, 18)]. The curve corresponding to (6) is shown by the 
dotted line (this distribution is observed in the case of serrated pulsations). 


Figure 2 


Thus, we may experimentally introduce the coefficient ko which represents 


the ratio between the amplitudes of the pulsating and traveling components of 
the field. However, according to (5) a theoretical introduction of ke is far 


from settling this question. It is impossible to determine the quantity Y, 
introduced in (5), theoretically. When we are dealing with real inductors 
having unwound sections, we may take the length of this section as the provisional 
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quantity Y. The error which is introduced by the currents which are closed at 
the longitudinal ends of the inductor is thus retained. When there are abso- 
lutely no shunting sections at the ends of the inductor, the introduction of the 
coefficient ko completely loses any physical meaning, in terms of the calcula-.. 


tion indicated above. The study (Ref. 3) proposes that ko be expressed as 


follows 
I 
ke=TFR’ (9) 


where 1, is the ratio between the magnetic conductivity of the shunting /99 


sections of the magnetic circuit and the magnetic conductivity of the active 

zone of the air gap. However, this happens very little in theoretical ratios 
(the problem may be reduced to calculating currents which are closed at every 
possible secondary route). 


The pulsating component in MHD machines, which is non-symmetrical over the 
inductor length, is only slightly expressed according to [Ref. 1-3], although 
this problem not yet been studied sufficiently. The presence of this component 
may be established experimentally by employing the measurement method developed 
in (Ref. 21), or by studying the distribution of instantaneous values of the 
field induction (Ref. 1). 


The distribution of the field in a disconnected magnetic circuit was 
studied by A. A. Lebedev (Ref. 6-9). However, as was already indicated by A. I. 
Vol'dek (Ref. 1), the results which he obtained do not have practical value due 
to the errors entailed. We shall discuss below the method which he proposed 
for balancing the pulsating field. 


As was already indicated in the introduction, the formation of supplemen- 
tary pulsating fields of the edge effect is undesirable, and they must be 
eliminated as much as possible. Employing the terminology given in (Ref. 3), 
we shall call this process adjustment of the field, and we shall call the wind- 
ings with a supplementary pulsating component adjusted windings. Different 
methods were advanced in (Ref. 3) for balancing the pulsating fields. In the 
case of a single layered winding, this is achieved by means of a special coil 
with a current which reaches the core of the stator at the winding level, and 
which is switched on in phase symmetrically with respect to the inductor core, 
or with a similar short-circuit loop. In physical terms, this is based on the 
fact that the field of the coil with the current reaching the core is constant 
over the inductor length in the case of u = œ». We need only select the number 
of loops of this coil in order that its pulsating field may equal, and be oppo- 
site to,the pulsating component of the inductor edge effect. With respect to 
the short-circuit loop, when there is sufficiently small resistance a current 
automatically passes through it in such a way that the total induction current 
through the surface (which it reaches) will equal zero. 


For a m-phase single-layered winding, according to (Ref. 3), the number of 
loops of the compensating coil is 
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wy = Wohi 
T P 
where Wo is the number of winding loops; kyi is the winding coefficient; p is 


the number of pairs of poles. 


It is thus assumed that the same amount of current passes along the loops 
of the compensating coil as along the loops of the inductor winding. 


For two-layered windings, it is assumed that the half-empty grooves /100 
at the inductor ends are filled by the corresponding coils, whose free ends are 
closed around the edge of the inductor, or are located in special grooves at 
the ends of the inductor. If the clearance is small, in the region of inverse 
closure of the compensation windings it is proposed that the clearance be in- 
creased and the winding distribution along the length of the section be de- 
creased, in order that the influence of these sections upon the secondary cir- 
cuit be minimal, The use of a short circuit loop is also possible. 


With respect to the compensation method proposed by Blake (Ref. 10), and 
reiterated by A. A. Lebedev (Ref. 6-9) -- which amounts to a smooth change in 
the loading to zero at the inductor ends -~ it is our opinion that total de- 
struction of the pulsating component may occur, It is true that the question 
remains open as to how great the drop must be, There is certainly validity in 
the opinion of A, I, Vol'dek (Ref. 1) that this leads to incomplete utilization 
of the magnetic circuit, and that this method is barely justified in compara- 
tively short inductors. 


The studies (Ref. 1, 3) present arguments regarding the error entailed in 
this method, and argue that the author (Ref. 10) is only interested in the 
field in the magnetic circuit yoke. However, as was shown in (Ref. 1, 3), when 
there is a purely traveling field in the clearance, an additional pulsation 
current may exist in the yoke. Actually (see Figure 1) in the case of B> 0 
and Y > 0 we may obtain the following from the expression (4) for the magnetic 
current in the yoke 


z 
, Bp 
opg f Bn sin (of — az) dz =-7 [cos (of — az) — (10) 


— (—1)? cos of]. 


However, the error entailed in the method of Blake does not follow -- i.e., 
the fact that a pulsation field may exist in the clearance if it is a purely 


traveling field in the yoke: 
B ey OBro 
clear Oz , (11) 


We thus find that the field in the clearance will be a purely traveling field, 
if it is a purely traveling field in the yoke. 


An article by A. P. Rashchepkin (Ref. 15) recently appeared, in which he 


90 


studied windings with linear current loading which was not uniform over the in- 
ductor length (in other aspects, the problem was similar to that investigated 
previously in [Ref. 2], see Figure 1). The loading is represented in the form 
of a superposition of the windings with different pole divisions (t) and /101 
the same values of the frequencies (w) and amplitudes of the linear current 
loadings (Jo) The author (Ref. 15) obtained expressions for the magnetic in- 


duction in the clearance in the form of a series, and studied the possibility 
of balancing the pulsating fields. The method advanced does not differ basical- 
ly from the method proposed by A. I. Vol'dek which was just discussed. 


The distribution of the field in the clearance of a finite inductor, as a 
function of the inductor length, was studied in (Ref. 20). This article in- 
vestigated the problem, which was similar to that solved by G. I. Shturman 
(Ref. 2), if we set B and Y = 0 in the latter, and if it is assumed that p 
(length of the inductor in pole divisions) is an arbitrary number. Somewhat, 
unexpected results were obtained. In this model, the pulsating field was en- 
tirely determined by the inductor length. The generalization of expression (4) 
was also given for the field in the clearance in the region of the wound sec- 
tions: 


B= Bm [cos (wt—az)—cos mp cos ot], (12) 


where p, as was already indicated, is half of the inductor length in units of 
the pole division. 


In order to suppress the pulsating component, it is proposed that induc- 
tors be used whose active portion has an odd number of pole divisions. In this 
case cos TP _ 0 (p -- odd number) and the pulsating portion of the field dis- 


appears in(12). In the case of even 2p, we obtain cos tp = (-1)P, which co- 
incides with the results given in (Ref. 3-5, 17). The article (Ref. 20) pre- 
sents the results derived from experimental studies, which substantiate the 
validity of (12). 


Similar ideas were advanced in (Ref. 10, 19) regarding the possibility of 
a purely traveling field in the clearance of a finite inductor, However, a 
sufficiently comprehensive theoretical basis was not provided for these ideas, 


3. Field of a Finite Inductor with a Secondary Circuit 


All of the preceeding results pertain to no-load operation of the machine. 
However, the construction of a MHD machine must be optimal in the operational 
regime. Consequently, it is of greatest interest to study the distribtuion of 
the field in an inductor with a secondary circuit. If a comprehensive study of 
this type were carried out, it would be possible to refine the efficiency with 
which the field could be adjusted during no-load operation. However, as was 
already indicated in (Ref. 1), this problem has been studied to the least ex- 
tent, and there is no unified terminology for it. We shall deal with the study 
performed by A. I. Vol'dek (Ref. 12) in greater detail later on, retaining the 
terminology employed by him (Figure 3, previous notation). The difference from 
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the problem investigated above consists of the fact that the magnetic circuit 
is not limited in terms of height (h > œ) and in the longitudinal direction/102 
(y >œ). Its magnetic permeability if infinite (u > ©) and, in addition, the 
clearance is filled by a conductor having the conductivity o which moves like 
a solid body at the constant velocity u = ul. The study (Ref. 12) also takes 


into account the presence of a clearance between the inductor and the metal, 
like serration of the inductor, but this may be reduced to a change in the con- 
ductivity of the body which fills the entire clearance. We shall not deal with 
problem of whether this conductivity of the metal is real or effective, and 
shall take into account additional phenomena. 


Figure 3 Figure 4 


Assuming that the field is a plane-parallel field (i.e., it does not 
depend on the x-coordinate), for H = He we obtain the following equation 


CH OH ðH 
Ge TM (G7 + 032) =o. a 


The author (Ref. 12) writes the solution within the limits of the active 
zone as the sum of two fields: 


H=H,+H, (14) 
where 
Hy = Ay eiot-22) 4 H pelot, (15) 
Thus Hy is given and differs from zero only within the limits of the active 


zone (|z|<pt). 


In our opinion, it would be simpler to give the linear loading in the form 
of a traveling wave =- i.e., only the first term in (15), The pulsating com- 
ponent of the no-load operation [second term in (15)] would have to be obtained 
from Hy in (14), when o > 0. This problem is more general as compared with the 


study by G. I. Shturman (Ref. 2), if Y > œ, u >œ, h > æ in the latter, and the 
results given in (Ref. 2) must be obtained as a special case of the problem 
under consideration. Such a distinction would be valid if we wished to con- 
sider the finiteness of the magnetic circuit or the height of the yoke. How- 
ever, we would then have to define the pulsating fields beyond the limits of 
the active zone; these fields exist beyond this zone, as was shown by G. I. 
Shturman. 


For the given case (i.e., B > O and Y > œ), we obtain /103 


——— 


92 


Hin= — (—1) Aig. 


(16) 
Let us write the formulas derived in (Ref. 12), with allowance for the 
statements presented above: 
H, = Cye(*+B)2, (17) 
H = Ce- 0-0, (18) 
Hm= Cie %0 Cete Hogeta H ge-ta, (19) 
where 
zs iesHis. 2. aoc on, 
Hea a Fies’ x VB +ie=a +ib; 
p= Bert ; €= W00; (20) 
jats b= E; d=Vp+yp FE; 
y2 Y2d i 
gain tia (1Ha + AÊ (—1)e+ Ha | e-0-MPr; (21) 
21 20 
Cafè =E E (yn + AEP Hi| eat, eo 
C3= —2Cye(-B9* sh (4— B) pr; (23) 
C= —2Cye™+B)Pt sh (A+ B) pt. (24) 


The indices I - III designate the total fields in the corresponding regions 
(see Figure 3). 


Let us study in greater detail the structure of the field in the active 
section (|z|<pt). We shall only be interested in the additional terms in ex- 
pression (19), i.e., the first and second terms. The last two terms occur in 
infinite inductors, and we shall not discuss them. 


Thus, with allowance for dependence on time, we may write the field of the 


edge effect Hy sige in the following form 


Ho edge =Cie-la-Bzellot-bz) + Cyela+Przei(ot+b2), (25) 


As may be seen, the fields of the edge effect are fields which move in 
opposition from the ends of the inductor at a velocity of 


-+2 + d 
“edge gar a= ae co 


In (Ref. 12), the first field, which is propagated in a positive direction,/104 
is called the direct field of the longitudinal effect, and the second field is 
called the inverse field. As may be seen from (25), in the case of B # O (this 
corresponds to u # 0), the fields of the edge effect are damped in a dissimilar 
way: in the case of u>0 the direct field is damped more slowly. The field is 
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carried along by a flux of molten metal, which has been observed experimentally 
(Ref. 13). 


Such a phenomenon has been observed by the author beyond the limits of the 
active zone. A wave escapes from the right end, which is similar to the wave 
of the direct edge effect, and a wave which is similar to the wave of the in- 
verse edge effect escapes from the left end. In order to obtain a clearer pic- 
ture of this, let us examine Figure 4 given in (Ref. 12), which presents a dia- 
gram of the distribution and motion of the fields of the edge effect. The 
arrows designate the direction in which the waves of the edge effect are prop- 
agated; their velocity determines (26). Let us present two numerical examples 
from (Ref. 12): (1) induction pump for aluminum at 735°C; (2) induction pump for 
sodium at a temperature of 500°C. In both cases, it is assumed that f = 50 cps, 
the pole division is t = 0.15 m, and the number of pairs of poles is 2p = 6. 
The initial data for these examples are presented in Table 1. 


TABLE 1 


INITIAL DATA FOR THE EXAMPLES 
"Pump for U 


Quantities Notation Dimensions Aluminum Magnesium 
Layer Thickness | a) | mwm | 10 eee, RPS 
Equivalent Clearance ks mm 35 6.6 
Specific Resistance 1:0 10-80hm-m 21.3 18.44 
Parameter € 1/m2 528 1620 


Table 2 presents the quantities for the pumps under consideration which 
characterize the fields of the edge effect. 


TABLE 2 /105 
CHARACTERISTICS OF THE FIELDS OF THE EDGE EFFECT 


Quantities 


Dim- Pum for Aluminum! Pump for Sodium _ 
ensions} ser | s-o3 | s-005 | s=1 | s-03 | s~0.05_ 
$ i/u | 0 885 | 12,0 0 27,2 36,8 
a i/u 16,2 17,4 18,3 28,4 35,7 41,6 
b 1/4 16,2 15,2 14,4 28,4 22,7 19,5 
(1) 4p 2) aw/cex 19.4 20,6 | 21,8 | 11,0 13,8 | 16,1 
a-B | l/m 16,2 8,55 6,3 28,4 8,5 4,8 
a+ß 1/4 16,2 26,2 30,3 28,4 62,9 78,4 
Flay : Mig | — 0,772 0,34 0,06 0,965 0,743 | 0,182 
(1) - u ; (2) - m/sec 


edge 


94 


As may be seen from the table, the velocity of the edge effect fields in- 
creases with a decrease in the slipping and conductivity. Damping of the fields 
is characterized by _1 and _1 (corresponding to the distance at which the 

B+a a-8 
inverse and direct fields are damped by a factor of e). Table 2 also presents 
the ratio between the field of the secondary currents and the field of the pri- 
mary currents. 


It should be pointed out that this phenomenon is nothing else than diffu- 
sion of the magnetic field from each end of the inductor in both directions. 
Passing to the limit o > 0, u > 0 in expression (25) we obtain 

Hy edge Aaa iger (27) 
i.e., the field of the longitudinal edge effect of the primary circuit. This 
result. is obtained because the field of the direct and inverse longitudinal 
effect in a secondary circuit is not damped in the limit over the length of the 
inductor, and its velocity strives to infinity ~1_ in the case of o > 0. In 
o 


regions I and III, the fields disappear, since C3 and C, strive to zero when 


u > 0 and o > 0. 


The longitudinal edge effect in a finite inductor with a secondary circuit 
was also studied in (Ref. 14) (Figure 5), where -- in contrast to (Ref. 12) -- 
it is assumed that there is winding on both sides of the channel, the magnetic 
permeability of the magnetic circuit is u = const, the inductor length is 
arbitrary, and in addition damping of the field in the clearance is taken into 
account (surface effect). 


The problem is solved by means of the vector potential A, which is found 
from solving the equation /106 


AA — wos (34 24) o (28) 


with the corresponding boundary conditions. 


The author (Ref. 14) makes a 
Fourier transformation with respect to 
the z variable -- i.e., the solution 


Z for A is written in the form 


- - 1? 
Figure 5 A= zz f A, (£, x) ečzdĘ. (29) 
—o 


In view of the symmetry with respect to x, the Fourier component in the 
clearance may be written as 


A; (8, x) =Ao ch (V§?—iupo+iopa) x. on 


However, the assumption is also advanced that in an upper magnetic circuit the 
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given component of the potential is 
Ayy (È, x) = Boe-*. (31) 


It may be readily seen that, by employing the Fourier component (31) we 
obtain the following for the potential in region IV 


œ 
1 k 
Aw=5~ f Boo-tds, (32) 
—o 


Since x > 0 in the upper magnetic circuit, this integral diverges, because the 
integrand increases indefinitely in the case of —€ + -~, and consequently 
Aly (€,x) does not have the form of (31). 


Since the author (Ref. 14) employs the incorrect Fourier component to ob- 
tain all further results, the analysis of this study could be discontinued. 
The incorrectness of the results, in connection with the statements given above, 
pertains to all formulas and arguments when u of the magnetic circuit equals 
const. However, a large portion of the article(Ref. 14) is devoted to studying 
the case when u of the magnetic circuit equals infinity. In this case, the 
terms, which are obtained from the incorrect form of the potential in /107 
region IV, vanish, and a correct result is obtained. The magnetic field 
in the clearance may be expressed in the form of a series, whose first term 
coincides with the results derived in (Ref. 12). The author (Ref. 14) has also 
obtained analytical expressions with allowance for the first term of the series 
for Joule heat losses for the effective power and consumed power. In addition, 
the results were analyzed qualitatively for specific values of the parameters 
introduced, but we shall not present these expressions due to the fact that 
they are cumbersome. 


The author (Ref. 14) has also studied the problem of compensation for the 
longitudinal effect -- i.e., the possibility of obtaining a field in the con- 
ductive metal, in the form of a non-distorted traveling wave. In order to do 
this, it is proposed that winding be employed in the forms of the superposition 
of three windings with different t and different amplitudes of the current 
loadings. However, the use of such complex windings requires an additional 
theoretical basis or experimental verification. 


The field of a finite inductor with a secondary circuit was studied in 
the work (Ref. 11), but the results obtained (which was indicated in [Ref. 1]) 
have no practical importance. 


The longitudinal effect of a finite inductor without a secondary circuit, 
and with a secondary circuit, was only studied theoretically in a qualitative 
manner. To obtain a quantitative theory for the longitudinal effect of a 
finite conductor, we would have to find more precise methods for determining 
ko in expression (4). With respect to the secondary circuit, it would be desir- 


able to clarify the influence of a clearance between the inductor and the molten 
metal layer upon the decrease in the fieldsof the edge effect in the secondary 
circuit (25), since the results given in (Ref, 12, 14) were obtained upon the 
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assumption that the entire clearance was filled with a molten metal. 


It is desirable to provide experimental confirmation of employing inductors 
with an odd number of pole divisions in the active region, when there is a 
secondary circuit. 


Concurrently with these studies, an analysis must be performed of the 
longitudinal and transverse edge effects in MHD machines. This problem was 
studied in (Ref. 16), but -- as was pointed out in (Ref. 1) -- these studies are 
far from complete. The author (Ref. 16) has provided a very approximate solu- 
tion of the problem, and only relates the phenomena occurring in the molten 
metal, within the limits of the inductor active section, with the longitudinal 
effect. 


It would be interesting to study the influence of inductor projections up- 
on the longitudinal effect. 


The directions to be followed by future research, which we discussed above, 
are naturally only provisional, and are far from encompassing all possible paths 
of research on the longitudinal edge effect in MHD machines. 
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PONDEROMOTIVE FORCES ACTING UPON CONDUCTIVE BODIES IN THE 
TRAVELING MAGNETIC FIELD OF A CYLINDRICAL INDUCTOR 


Yu. K. Krumin' 


1. Introduction 


Up to the present time, the traveling magnetic field with axial symmetry 
has been studied less frequently than fields of the flat type. However, it is 
of great importance for the development of new designs for MHD machines. The 
advantage of fields with axial symmetry lies in the fact that inductors for pro- 
ducing such fields have no frontal sections, and there are no transverse edge 
effects in the working medium. This leads to an increase in the efficiency of 
the machine, as a whole. 


However, cylindrical inductors in electromagnetic pumps and other MHD mach- 
ines have been employed to a lesser extent than flat inductors. There are few 
articles in the literature on this problem (Ref. 19, 22). The use of a travel- 
ing magnetic field of cylindrical inductors is also well known for producing a 
different type of electromagnetic conveyers of solid conductive objects (Ref. 6, 
7). The reason for this may be found primarily in the fact that the distribu+ 
tion of the force density over the channel cross-section is different in a 
cylindrical pump (with a solid channel, without an inner core) than it is in a 
flat pump: it equals zero in the center of the channel. Due to this fact, a 
different type of closed flows may occur, which obstruct the normal operation 
of the pump. One opinion even holds that a cylindrical pump without a core can- 
not operate successfully, in general. Only recently have studies appeared (Ref. 
2, 17), in which it was shown that this opinion is incorrect, and that the con- 
struction of cylindrical pumps without cores is fully possible. 


A ferromagnetic core considerably improves the situation, and therefore it 
is employed in all existing cylindrical pumps. However, it is only employed 
under conditions in which the core still retains its ferromagnetic properties. 


One of the main problems which is encountered when designing any device 
utilizing ponderomotive forces, which influence the conductive media located in 
the traveling magnetic field, is the development of methods for calculating 
these forces, particularly their maximum values, as a function of other /110 
characteristics of the device. This chapter will be devoted to these problems. 


2. Cylindrical Inductor and its Electromagnetic Field 


As is known, the cylindrical inductor of a traveling magnetic field re- 
presents a system of coils located on a common axis, which are supplied by an 
alternating current which is usually a tri-phase current. The field produced 
by such an inductor is rather complex. Therefore, in theoretical designs sim- 
plifications are usually introduced which assume that the conductive medium 
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under consideration is located in the idealized field produced by an ideal in- 
ductor, and not in the field of a real inductor. A cylindrical surface having 
the radius R of infinite length, around whose circumference a surface current 
flows, is placed under an ideal inductor. Its density changes sinusoidally 
both in time and along the z-coordinate: A = Ajet(wt-az). In other words, this 


current represents a traveling wave, whose phase propagation velocity equals 
w/a. At such a velocity, the equal phase plane of the field excited by 
this current is displaced along the z-coordinate. 


The electromagnetic field in this inductor -- i.e., in the case of r < R -- 
has the following form if the inductor is filled with a uniform medium having 
the specific conductivity o and the magnetic permeability u: 


H, = aos, (Br) piot -az). 
B BR)" 


(la) 
H,=Hoo jee ellot -az); j 


= — fopHo , (Br) ol-az 
E ROA ae 


Here we have 


B= Ya? Fiopo, (2) 
and Hoo is the strength of the tangential (H,) component of the magnetic field 


on the inductor surface ~- i.e., in the case of r = R. 


If there are no conductive and ferromagnetic media within the inductor, 
the field has the following form 


A, = iH 4 Aca eilot-az); /111 
H= Hoo 7S. eilot-az)s (1b) 
= — iopoHoo (ar) sitot-az 
Ey sa To(aRy eilo ), 


Figure 1 presents a picture of the magnetic field of an empty inductor for 
two cases, when the ratio 2R/t equals 2 and 0.5. 


The electromagnetic field of a real inductor may be computed by sum- /112 
ming up the fields of all the loops along which the current passes. Such a 
calculation was performed in (Ref. 16), in which an expression was found for the 
radial (HL) component of the magnetic field for an inductor of finite length. 


It was also assumed that the inductor winding was infinitely narrow, and the in- 
ductor itself represented a group of solenoids located on a common axis and 
arranged close to each other. 


Since the magnetic field of a solenoid of finite length is well known (Ref. 
14), no particular difficulties are entailed in determining the total field of 
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Figure 1l 
Form of the Magnetic Field in a Cylindrical Inductor: 


a - For 2R/t = 2; b - For 2R/t = 0.5 


the inductor. Employing dimensionless notation, we may represent the radial 
component of the field of such an inductor in the following form 


~ tale -5V7 ((} -#) ken -2 zen |- 


= a Ae -g E (ba) } cos [oe (a= tye] l 


(3) 


Here E and K are the complete elliptic integrals of the first and second 
type; 
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(1+7)2+ Ej (4) 


(5) 


The values of HL and a are not given in the work (Ref. 16). 


Figure 2 presents a graphic illustration of the function (3) for the moment 
of time t = O at a value of r = 0.5 for certain T. 


The field of a real inductor differs from the field of an ideal inductor 
(if we do not assume that it is distorted close to the ends) in the fact that 
in a real inductor it is not a sinusoidal field along the direction in which it 
is propagated -- z-axis. It may be expanded in series of spatial harmonics, 
which all change at the same rate over a period of time, but which have differ- 
ent t. The larger the ratio of the pole step to the inductor radius, the larger 
is the amplitude of these harmonics. This is clearly illustrated in Figure 2. 
In the case of small Tt, the form of the curve is close to a sinusoid, but its 
amplitude is small. The conclusion may thus be reached that a certain optimum 
* optim must exist, at which the operation of the electromagnetic device will be 
the most efficient. 


The distribution of the magnetic field (averaged over time) in the {133 
real cylindrical inductors along the inductor length was studied in (Ref. 3). 
The influence of individual coils and projections of the inductor magnetic 
circuit upon the field may be clearly seen in the graphs. 


In addition to studies which calculate the electromagnetic field in an in- 
ductor filled with a unifrom substance, there are still many solutions of the 
problem when a conductive medium having a certain configuration is placed in 
the inductor. We shall investigate these solutions, together with the corre- 
sponding solutions for the ponderomotive forces, since they represent two sides 


of the same problem. 


3. General Formulation of the Problem and Method of Solution 


As has already been indicated, the problem of a cylinder having infinite 
length, located in the field of an ideal inductor, has been subjected to the 
most extensive theoretical treatment. Although this represents an idealization 
of reality, the results obtained nevertheless present a correct concept of the 
phenomenon -~ the nature of the dependence of force on frequency, conductivity, 
pole division, etc. The contribution furnished by a consideration of the real 
conditions amountsto a correction of the solution obtained. In addition, it may 
be stated that -~ if we know the spatial harmonics of the field -- then, in view 
of the superposition principle, we may obtain the complete solution by summing 
up all of the special solutions which correspond to each harmonics. All of this 
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Figure 2 Figure 3 


Dependence of H for a Real Inductor Cross-Section of Coaxial System 


Haying Finite Length Upon oa thë When the General Case if Calculated. 


Case of r = 0.5 and Certain Values 
of t: 1 - T=0.5; 2 - t=2.0; 
3 - t=4.0. 


speaks in favor of a detailed study of the idealized case, which may thus be 
assumed to be fundamental for a study of real cases, 


The solutions for ideal cases, which are presented in different stud- /114 
ies, differ primarily in the amount and order in which the layers are arranged 
(layers of coaxial cylinders), which are located in the inductor. For all of 
these problems, we may present a general method which may be employed to derive 
the solution, which is as follows. 


Let us assume the most general case, when the device consists of n coaxial 
cylinders, along which circular surface currents pass (we shall assume that they 
have the same w and a). These currents alternate with m uniform layers, whose 
specific conductivity and magnetic permeability equals Oy and Hk kS ly 2, ...; 


m) (Figure 3). The radii of the surface dividing the media may be designated 
by r. 
k 


Let us solve the Maxwell equations in cylindrical coordinates. Since the 
phenomenon does not depend on the angle p in this case, all of the derivatives 
with respect to ¢ in these equations must be set equal to zero. Due to this 
fact, the system of equations is broken down into two independent systems (each 
one consists of three equations with three unknowns). One of them contains HL» 


H and Ep? and the other contains EL Pz and H It may be readily seen that 


$° 
E = Ez = Hy = 0, and the solution of the first system has the following form 
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Ha= on [Cial (Bar) + ConK1 (Bar) J efet-a2); 


Ha= fe [Cinto (Bar) — Car Ko (Bar) ] efot-22); 7 


Egr =[Cirl (Bar) + ConKi (Bar) ] et-), 


where 
Pa= ya Fioro (k=1, 2,..., m). (7) 


Such a solution must be written for all m layers, substituting the corre- 
sponding values of Tk and Uys with the exception of the extreme values — of 


the inner layer containing points with r = 0 and the outer layer containing 
points with r = æ. In order that the solutions may be finite and unique, we 
must discard the functions K in the inner layer, and the functions I in 
the outer layer. 


The integration constants Cik and Cok are determined from the boundary/115 


conditions, which hold on each dividing surface, and which have the following 


form 
pall, rk = urn H Ts k+l 


Hir — Az, k= AR for r=fp. (8) 
Eọr=Ep, k+1 


It may be readily seen that the first and third conditions of (8) yield 
identical equations, and we thus obtain a system of equations from 2k - 2 equa- 
tions with 2k - 2 unknowns. With respect to the second condition of (8), if no 
surface current flows along the dividing surface, then zero must be substituted 
instead of AL: 


The solution obtained completely determines all of the electromagnetic 
processes in the system under consideration. However, the complexity of the 
computations increases greatly with an increase in the number of layers, and 
therefore only problems with several layers have been investigated more or less 
extensively at the present time. 


Cylindrical pumps with the number of layers reaching 6 may be frequently 
encountered in practice (Ref. 1). In this connection, E. K. Yankop (Ref. 20) 
proposed a method for solving this system of equations in the case of a large 
number of equations, 


4. Methods of Representing the Results 


The direct purpose of the computation is to determine the ponderomotive 
forces influencing the working layer in the inductor. 


In this case, there are two force components which differ from zero ~- the 
radial component fy which contracts the cylinder, and the tangential component 
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f whích operates in the direction of motion of the field, Their density 


averaged over time may be computed according to the well known formulas: 


f= ReligoB.o"]= + Re [EqoBu" 
5) 2 (9) 
f= —4 Re [jooBn"]= — £Re[EqoB.0*l. 


The total force may be found by integrating the force density over the 
entire volume of the conductive substance, Thus, the radial forces are mutually 
balanced, and therefore produce no effects if the cylinder under consideration 
is solid. When they influence a liquid, different radial flows may arise. 


The force component operating in the direction of motion of the field pro- 
duces the useful effect of displacing the conductive body, Therefore, /116 
when one speaks of the forces in the traveling magnetic field, it is usually 
this component which one has in mind, unless otherwise stipulated, 


We must now turn to the method of defining the boundary conditions. It 
must be stated that (8) is not the only possible method for defining them, In- 
stead of defining the current values in the inductor, we may define the strength 
of the magnetic field on any dividing surface, This is the customary procedure, 
especially when there is only one surface along which the surface currents pass, 
Thus, the strength CH,» HL) is given, both on the surface where the currents 


pass, and on the surface of the conductive cylinder which is located in the in- 
ductor, In other words, one and the same problem may be represented in a com- 
pletely different manner, and the expressions for the force will also be differ- 
ent. When we are computing the force, we are only interested in the processes 
in the conductive medium under consideration, Therefore, we shall assume that 
it is more advantageous to divide the problem into two parts: to define the 
strength of the magnetic field (one of its components) on the surface of the 
conductive medium (if it is solid), with which the electromagnetic processes 
within it will be clearly determined, and to calculate this strength by the 
method indicated above as a function of the currents flowing in the inductor 
windings, as well as other characteristics of the device, It is possible to 
define the strength of the field on the surface of the medium, because the 
distribution of the electromagnetic fields does not change as a function of the 
surface on which we define the boundary conditions, 


If the medium which we are considering is a hollow body, it will be in- 
sufficient to define the magnetic field only in the form of a single boundary 
condition on its durface, We must here define two conditions, This may be 
done by defining the values of both field components from one side of the body, 
or with respect to one component from the outer and inner side. However, such 
a method is not employed in practice. By employing this method, for a specific 
form of the body we shall only have one formula for the force; we shall have to 
substitute a specific value of the magnetic field strength in each case in this 
formula. This makes it possible to compile the general curves or tables for 
the calculations, which may be employed in every case, and to analyze the phe- 
nomenon, abstracting from the foreign influences introduced by other elements 
of the device. 
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Naturally, all of these statements only pertain to problems regarding in- 
finitely long cylinders in an ideal inductor. 


In order to decrease the number of variables, the computational results 
may be represented in dimensionless form. Thus, the following dimension- /117 
less numbers are introduced as the arguments: 


@=opoa?, G=aa, (10) 


where a is a certain characteristic dimension. Customarily, if a solid cylinder 
is being considered, its radius is taken as a. For an empty cylinder, it is 
natural to use its wall thickness as a. In several articles w is called the 


relative frequency. 


In addition, it is advantageous to intorduce one dimensionless condition -- 
the complex relative frequency wj: 


z Aree itt ae 
aim pa= |/ tate tyta piyee talot y tiya (11) 


The quantities w and a have a simple physical meaning: they indicate that 


a variable field penetrates to a portion of the characteristic dimension a. 


Thus, w characterizes the penetration of a variable, non-traveling field, and 


the real portion wi characterizes the traveling magnetic field (Ref. 8): 


where hy and hy represent the effective penetration depth, respectively, of the 


variable non-traveling field and the traveling field in half-space with a plane 
boundary dividing the media. 


As regards representation of the force in dimensionless form, the situation 
is more complex, since different variants are possible. In this connection, 
different determinations of the dimensionless force are encountered in different 
articles, 


We shall assume that we have calculated the total force on a segment of a 
cylinder having the length 2. Intorducing the dimensionless parameters w and 
a which have already been indicated, we find that this force may be represented 
as follows 


re BH F(a, a). (13) 


If we divide both sides of this equation by the factor in front of £(w,a) 
and if we designate the left side by F (as, for example, was done in [Ref. 12]), 
we then find that F depends on a ~~ a quantity on which a also depends /118 


(in addition, w and a depend on a). This renders an interpretation of the 
results more difficult. 


106 


We may select a dimensionless condition F which does not include any of 
the quantities contained in w and a. This condition is as follows*: 
F 


P= H, (14) 
o 


However, Foy does not have the dimension of force here, although in physical 


terms it may be regarded as a force acting upon a certain volume of substance 
(this will be demonstrated in specific problems). Far has the dimension of 


pressure, but at the same time it is not pressure which may be produced by an 
apparatus. 


In this section, we shall adhere to the opinion that it is better to pro- 
ceed with a certain complication of the physical meaning of F, than it is to 
apply interdependent conditions, whose relationship is difficult to interpret 
in physical terms, 


The separation of the problem into two stages, as was mentioned previously, 
is also advantageous when one wishes to present a graphic illustration of the 
results, The relationship between no more than three quantities may be depicted 
in the form of a set of curves on a plane, but the solution of the general 
problem, even in the simplest cases, usually depends on no less than four quan- 
tities, Nomography would be of considerable help, but nomograms have not as 
yet been compiled for the results investigated in our article. A graphic illus- 
tration is completely impossible, since all the formulas are very complex for 
computations by hand. 


We would like to point out the following. If we use w to designate the 
angular frequency of a current which is supplied to an inductor, all the for- 
mulas of this article reflect the processes occurring when the conductive media 
are not in motion, When the conductive medium moves with respect to the in- 
ductor, all the formulas remain in force, but -- instead of w -- we must sub- 
stitute another yalue of the frequency w in these formulas. This value de- 


pends on the slipping S: 
O1=0S. (15) 


In the case of S = 0, when the medium and the traveling field move at the 
same velocity, all the forces are equal to zero. In the case of S < 0, these 
formulas will also be real, but the value of the forces will have to be assigned 
the opposite (negative) sign. f i 


5. Solid Conductive Cylinder in the Field of an Ideal Inductor /119 


An attempt was made to solve this problem in a report by A. I. Tyutin 
(Ref. 18). However, he only calculated the force density, and did not integrate 
it over the cylinder volume. This problem was studied in greater detail in 
* Strictly speaking, F and w depend on u. However, since we usually only deal 
with substances in which u = uo? in this context they are independent. 
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(Ref. 12, 13), where the authors calculated the force acting upon a solid con- 
ductive cylinder with the radius p, located in the traveling magnetic field of 
a coaxial inductor. If we assume that the z~component (having the amplitude 

Ho? of the magnetic field strength is given on the cylinder surface, we obtain 


the following result 


pa For im Bh (GN EM 


pH? joilo(0:) [2 s (16) 


We shall employ Eor to designate the force acting upon a volume of substance 


of 2p units of volume. Figure 4 presents a graphic illustration of dependence 
(16), where values of the condition 2pọp/t = 20/7 which are more suitable for 
practice are given, instead of a. This dependence represents curves which have 
one maximum. 


In the case of w <<a, the force 


increases in direct proportion to w, 
according to the approximate formula 


r Hi(@) Loa) +7 -OE 0- 
a?l? (a) (17 y 


and in the case of w >a it strives 
to zero in inverse proportion to ys: 


=. yia /120 
Fru ie f (18) 
The maxima of these curves are of 

special interest. However, since we 
still do not know how to derive. the ex- 
pressions for them analytically, we 
shall confine ourselves to analyzing 
their behavior in a figure. With an in- 
crease in the ratio between the cylinder 
diameter 2p and the pole step of the inductor t, the maxima F strive to the 
largest possible value, equalling approximately 0.85. For small ratios of 
2p/t, all the maxima occur for one specific value of w, equalling 6.33. For 
small values of w, the maximum force is reached at one specific value of 2p/t, 
equalling 0.92 (in the case of w £15; 2p/t < 1.5). 


Figure 4 


Dependence of F on w for Different 
2e/t for a Solid Cylinder 


The study (Ref. 13) presents formulas and curves for calculating Ho with 


respect to the inductor characteristics for the following case: the diameter of 
the inductor is R; there is a non-magnetic, non-conductive clearance between the 
conductive cylinder and the inductor; the space outside of the surface current 
(in the case of r > R) is filled with a substance having the specific conduct- 


ivity o' and the magnetic permeability u'. We then have 
E EEA E ON E: -A AEE A EERE 
A afem Cofes -sS | 6 Ee CORI } we 
“| ilo (oi) | aK (0s) í pa’ Ki (o's) a 
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where S,=1o(a) Kı (€) + Ko(a) l(a’); 
S2= [lo (€) Ko(a’) — Ko (@) lo (a’); 
Ss=I (a) Ki (@) — Kı (a) I (@); 


= = ett ooo (20) 
S4=/, (a) Ko(a’) + Ki (a) /o (a’). 
Here we have 
D’ = pR] rrita +i1/Yo?+a"— a" ; 
_ (21) 
o=o WoR? a/=aR; 
sate ce 
p=— =—, 
m "Tio 
Substituting the value of Hy from (19) in the expression for F, we must 
select the modulus EA . In practice, we are primarily interested in two cases: 
when the space r > R is filled by a substance with u' = œ and w' = 0. /122 
This occurs in an approximate manner when there is an iron magnetic circuit in 
the inductor, and when u' = Ho? o' = 0 (inductor without a magnetic circuit). 


In these special cases, the number of variables is reduced to 3, and it is pos- 
sible to express them graphically. 


First case. u' =œ, o' = 0 outside of the inductor. Formula (19) changes 
into 
Hy _ 1 
Ay S uhle al (22) 
Ao z|- BN (i) 54 5, | 
oilo (01) 
Second case. wu' = Ug» o' = 0 outside of the inductor. Formula (19) 
changes into 
TE hn nwt whut ie Tach 
Ap -fal (%:) [Ko(a’)S, Ko(a’) Ss i (23) 
ai= m = ~S2}|—-|— = — -S, 
Oilo(wi) | Ki(a’) Ky (a’) 


Introducing the_condition R = R/p, which is more convenient in practice, 
instead of a' (then a' = aR), (22) and (23) may be graphically illustrated as 
is shown in Figure 5, a and b. Due to the fact that these functions depend 


slightly on w, it was possible to plot all of the curves on two graphs. We 
should note that the values of Hoto in the case of w = 1 and w = O differ by 


less than 0.1%. The study (Ref. 12) compares formula (16) with the results ob- 
tained experimentally. Measurements were performed with cylinders having a 
different length. It was found that the force acting upon a unit length is 
greater for a short cylinder than it is for a long cylinder, in the limiting 
case. Extrapolating the results obtained to an infinitely long cylinder, there 
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was satisfactory agreement with the theoretical value of (16). 
6. Hollow Conductive Cylinder in the Field of an Ideal Inductor 


A problem of this type was first studied in the article by I. A. Tyutin 
(Ref. 18), which was already mentioned. A fairly complex problem was solved: 
it was assumed that the cylinder is located in the field of a composite in- 
ductor, consisting of two coaxial surfaces. Currents flow along the surfaces, and 
one of them is located outside of the conductive cylinder, while the other is 
located within the conductive cylinder. However, this solution was limited to 
only describing general formulas. 


The study (Ref. 12) derived a formula for computing the force acting /123 
upon a hollow cylinder with a wall thickness d (outer radius Tos inner radius 
r4). The inner space of the cylinder is filled with a non-conductive, and non- 


magnetic medium. Taking into account the given z-component of the magnetic 
field strength on the outer cylinder surface, we obtain the expression for the 
force: 


Fe Fer_ 
= ASETET E — peH? 
— 2a Im f0: [a plo (x) Ti — il (x) Ta) [a plo (x) Ta— Dily (x) Tal") o 
(2—d) |o; [a plo(x)Ti— o: (x) Ta] |? i 


where x = a(l - d); 
Ts =1o(w1) Kı [oi (1 —@)]+ Ko (1) [os (1 —) ]; 
T2=19 [wi (1 —4)]Ko (01) — Kolo: (1—@)] Jo (©:); 
Ta=1, (0:) K, fa: (1 —8)]— K, (0) h [o:(1 —d)); 
T,= ow: (1 —d)] Kı (ai) + Kolo: (1—3)] 1, @:); 


(25) 


wi is expressed by means of w and o with the same relationship (11), if 


the outer cylinder radius is substituted in w and a: 


O= opore?; @=are; a=“. (26) 


2 


In this case, Foy designates the force acting upon a volume of substance having 


2d units of volume. If we set rS 0, r3 = P, it is natural that these formulas 


change into a formula for a solid cylinder. 

The study (Ref. 12) alse presented formulas for computing the field 
strengths, if it is assumed that the tangential component of the magnetic field 
strength on the inductor surface is given. They are quite complex, and we shall 
not present them here. 


The works (Ref. 5, 15) presented the results derived from a small number 


111 


of experiments, in which they measured the force acting upon a hollow conductive 
cylinder. They did not compare the results with theory. 


The results derived from a small series of measurements with respect to 
the force acting upon hollow cylinders were presented in (Ref. 4), but they were 
compared with a theoretical formula for a plane inductor. 


Finally, if the thickness of the cylinder wall is small as compared with 
its radius, it is possible to employ the formulas for a plane layer (Ref. 21), 
instead of the above formulas, without entailing any large amount of error (see 
the article by A. Veze and L. Ulmanis in the present collection). 


7. Solid Cylinder Having Finite Length in the Field of an /124 
Ideal Inductor 


The theory presented above pertained to cases when it was assumed that the 
conductive cylinder had an infinite length -- i.e., the influence of the finite 
length of the cylinder was disregarded. 


G. Kh. Kirshteyn (Ref. 9) discovered a method for solving a similar prob- 
lem for a solid cylinder having finite length, on the assumption that there is 
no air gap between the cylinder and the inductor and that the space behind the 
inductor has infinite magnetic permeability. 


This method is as follows. 


Let the cylinder have the length 2 (Figure 6). Let us divide all of the 
space within the inductor into three regions: z < -& (region I); -L« ze 2 
(region II); z 2 & (region III). 


In each of these regions separately, 
we shall try to find the magnetic field 
in the form of the sum of two fields 


Hy and H,» where HL is the field of a 


hollow inductor (lb) in the regions 

I and III, while in region II it is the 
field of the inductor filled with a 
uniform, conductive substance having 
the characteristics of the finite cy- 
linder under consideration, which is 
described by expressions (la). 


Figure 6 


Arrangement of Conductive Cylinder 
Having Finite Length in Inductor 


H, is the perturbation of the field 


H close to the ends of the cylinder. It thus follows that in the case of 


z > +œ the field H, must strive to zero. It follows from the condition that 


u = © outside of the inductor that the z-component of the magnetic field is 
the same over the entire surface of the inductor, no matter whether there is a 
cylinder there or not. As is customary, it is assumed that this component 
equals es = Apei(ut-az). However, since this condition is fulfilled 
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by the field Hs it follows that the tangential component of H, on the inductor 


surface equals zero. This is the boundary condition for computing H,. 


The following solution satisfies these requirements: 


In region I (z < -£) 125 
i æ 
X, 
H, = -aR QyCinbonl ( a ePan zelat. 
, @ 
La t ` m Kont lot. 
2 opoR £ Cintonto| R J tenzet (27a) 
æ 
Egy! = `, Cindi (3) ePanzelot. 
=I 


in region II (-2« z < 2) 


ke] 
i Xan 
ie > n 
Hn wyp L bind R 


(Coneoin? — C3ne7~ >in?) elot, 


co 
i Xon! l : (27b) 
H,"t=—— > Xeno") (Coneoin? + Cyne~bin2) elot; 
opR n=l R 
oo 
Ege"! = X ka ) (Canettin? 4- Cane tin?) eit, 
n=l 
in region III (z > £) 
i oo 
Haim t- X} Canbar c) e—tan zelot; 
n=l 
Ha!!! = i ‘ > CinXondo (x) e~*anzgiot . (27c) 
@poR n=l R 
œ 
Ep! = pe Canhi 5) e—banzeiot, 
n=l 
Here 
y 
bin=V/ 285 iono; ban =, (28) 


where Xon are the roots of the Bessel function of the first kind of zero order 


olon?) = 0); R is the radius of the inductor and the cylinder; Cia? t> Cyn {126 
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are the integration constants, which may be determined from the boundary con- 
ditions (for the total field) 


poll,’ = pH" uH"! = pof!!! 
H= H," for z=—i;; A, Ht t for z=l. (29) 
Eq! = Eg! E! = Eqit! 


The first and third conditions of (29) yield identical equations. Thus, 
four equations are obtained with four unknowns for each n, from which the con- 
stants may be determined. 


The study (Ref. 9) found the expression for the perturbation of the force 
caused by the presence of the cylinder ends: 


pane co 
T 2a2@? 1 
AF=— imf 5 a 
k — (Xon?+ a’) (Pn?+ a?) A 


T xon sin ak +a cos ak Aee 7 (30) 
z ch pak + xon sh Pak (Pn ch pak sin ak—a sh Pnk cos ak) + 
Xon COS ak—a sin ak zy in: sas 
Pn sh Pnk +Xon ch Pak (Pn sh Pnk cos ak+ach Pnk sin ak) D. 
Figure 7 
Dependence of AF on k in the Case of a = 1.5 for Certain Values of w 
Here /128 
wat fo. l aAF Ere. 
in Vito, Peles k=: = 
Pa eV Se Re Be OE ae aagl” (31) 
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AF designates the difference between the total force acting upon a finite 
cylinder, and the force acting upon a segment having the same length of an in- 
finitely long cylinder. 


In the case of k>>0.4, this formula may be simplified as follows: 


n— Xon + 02 272 +202 x Xon 24-a? 
e278 $ BK i ae axon “Ylu, |. 
mes] (xon +02) [ (xon +02)? + 0 á (32a) 


Un = VVO? txon +x on. 


Figures 7 and 8, a and b, present the graphic dependence of AF for certain 
values of w, a and k computed according to formula (30). It should also be 
noted that, for a sufficiently large value of a, AF may be negative for certain 
values of k. 


A comparison of the results obtained with the experimental data existing 
previously (Ref. 12) shows satisfactory agreement. 


b 


Figure 8 


Dependence of AF on a for Certain Values of k: a - In the case 
of w = 10; b - In the case of w = 15 [1 


However, it must be stated that this problem must stíll be investigated 
in greater detail. 


8. Conductive Cylinder in the Field of a Real Inductor 


U. A. Saulite, A. E. Mikel'son et. al. have studied the problem of the 
behavior of an infinitely long cylinder in the field of a real inductor, when 
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this inductor approximated the model which was discussed in Chapter 2 (Ref. 16, 
17). In this formulation, the problem quite accurately reflects the case in 
which the cylinder ends project sufficiently far from the inductor ends, where 
the field of the inductor has decreased almost to zero. 


It was shown in (Ref. 15) that, if the cylinder is length is small as com- 
pared with the inductor length, the force does not remain constant when the 
cylinder is displaced along the axis of the inductor. The force can even change 
sign at the end of it (Figure 9). The experiments were carried out with alumin- 
um cylinders. 


A calculation was made in (Ref. 17) of the radial component of the magnetic 
field strength in a system consisting of a real inductor -- whose model is de- 
scribed in Chapter 2 ~- and an infinitely long cylinder placed in it (Figure 10). 
It was found that He is expressed in dimensionless form as follows, with- /129 


in the conductive cylinder 


= _ 2nHn Es TKA) 
Ha= poo J (kr) 1, (Āri) + 


E T EAJAK Ar) p+ R Jo(kr)Ki Ari) 
A a ba n—An 
x > feosx (22—z} -cos [C722 shy, gae a 7 ‘lax, 


outside of the conductive cylinder in the case of p< r< R /130 


ane x Si (Ri)? Io(Ar1) w—R Jo(kriyh (Ari) 
H, 2e 1 
a efko ) T EDEKA a tEh (Fri) Ki (Ar) 


Ji (FF) Elo (A) a Engg ona), Ki (ar my |x ee 
a 


S > i - a—1)s0 
x>, barf ne A OSE —al a rye ba Say 


Syllis] 


n=l 


= (2 — ba) K éa) 2 a E (ka) |? -i [pt 23] 


The notation employed here coincides with the notation in formula (3). In 
addition, 


(33b) 


' (34) 


where P is the radius of the conductive cylinder. 


The report in question does not present the expressions for the tangential 
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Figure 9 Figure 10 
Change in the Force Acting Upon a ' Mutual Arrangement of an Inductor 
Cylinder Along the Inductor Axis; (F - Having Finite Length and a Conductive 


in newtons; x - in meters, t = 0.273 m) Cylinder 


component of the magnetic field, nor the expression for the electric field 
strength. The forces acting upon the cylinder have still not been computed. 


9. Conductive Sphere in the Field of a Cylindrical Inductor 


A sphere, as a body with the greatest symmetry, occupies a special position 
among all bodies having finite dimensions. Therefore, it is definite interest 
to determine the ponderomotive forces on both a solid and hollow sphere. The 
author of this article attempted to calculate the forces acting upon a solid 
sphere located in the field of an ideal, cylindrical inductor, under the assump- 
tion that u = © outside of the inductor (Ref. 10, 11). However, in formal terms 
the solution obtained is not suitable for the computations. 


Experiments measuring the forces acting upon spheres have shown that the 
dependence of the force change upon w Lw = opwp2 „p is the radius of the sphere 
if the sphere is solid, and w = ouwb? (b is the thickness of the sphere wall) 
if the sphere is hollow] is exactly the same as in all other cases. For small 


w the force is directly proportional to w. A maximum is then reached, and when 
w increases, it begins to decrease. 


Figure 11, a, b and c, presents the results derived from measurements /132 
performed with spheres having a different radius and made of a different materi- 
al -- copper, aluminum, tin, and lead. These spheres were placed in an inductor 
without a magnetic circuit (the quantity F/Ap2p2, which is a dimensional quan- 


tity, is plotted along the ordinate axis. However, in this case this is of no 
value). In the case when a = 0.31, it was found that the ordinates of all four 
curves were proportional to the radius of the sphere. Therefore, if we intro- 
duce the condition F/Ay 203, they may all be combined on one curve [see /133 
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Dependence of F/A9202 on w for Solid 


Spheres Having Certain Radii: a - in 
the case of a = 
a = 0.62; c - in the case of a = 0.31. 
The Inductor Radius is R = 4.0 cm; F is 
Given in Newtons; A, is Given in a/m; 

p is Given in Meters. 
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0.93; b - in the case of 


Figure 12 


Dependence of the Maximum Force on T 
and t for Solid Spheres. 


- in a/n; 
p/R. 


F is Given in Newtons; Ay 


max a 
p - in Meters; t = T/R; r= 


also (Ref. 11)]. However, for other 
values of a, this did not occur. This 
indicates that the above dependence is 
not universally valid. 


The maximum force is of great 
importance in designing different 
electromagnetic carriers. The study 
(Ref. 5) presents a graph which makes 
it possible to determine the value of 
Fal agr ee from tT, r and w, where 


T = T/R, and r = p/R (R is the inductor 
radius) (Figure 12). This graph was 
compiled on the basis of experimental 
results, and is incomplete in the sense 
that it does not encompass a large 
enough range of possible parameter 
changes. 


The results derived from measuring 
the ponderomotive forces acting upon a 
hollow sphere (Figure 13) were also 
presented in (Ref. 5). F' = F/Ay’b*, 


OG = onwb2; b is the thickness of 
the sphere wall. 


Figure 13 


Dependence of F' on w' for Certain Hollow Spheres 


We shall refrain from discussing the results in greater detail, since it 


is expected that a suitable method for computing the forces acting upon a sphere 
will be found in the very near future. 
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THEORY FOR THE PROPAGATION OF PULSED ELECTROMAGNETIC FIELDS 
IN MOVING CONDUCTIVE MEDIA 


G. Ya. Sermons 


1. Introduction 


The necessity of developing a theory for the propagation of electromagnetic 
pulsed fields in conductive media has arisen in the last twenty years due to the 
utilization of relaxation processes in ferromagnets (Ref. 1). This also pertains 
to the necessity of studying the influence of transitional processes which occur 
in electrotransmission lines when the current is switched on or off, as well as 
in telephone and telegraph lines (Ref. 2-3). A theory for the propagation of 
electromagnetic pulsed fields has been developed intensely during the last 50-60 
years due to the use of pulsed fields for geophysical research (Ref. 4-22). The 
essence of the so-called method for establishing a current (Ref. 23), which is 
used for geological exploration, consists of establishing an electric or magnetic 
field on the surface of the earth when a constant current is switched to a trans- 
mitting antenna. A horizontal, grounded conductor or loop is used as the trans- 
mitting antenna. In many cases, the distance from the receiving antenna to the 
transmitting antenna is large, as compared with the dimensions of the transmit- 
ting antenna. The transmitting antenna may be regarded as an electric (hori- 
zontal conductor) or magnetic (loop) dipole. 


Studies on the propagation of an electromagnetic field impulse in conductive 
media may be reduced to investigating the formation of the field of an electric 
or magnetic dipole which is located in a conductive medium (Ref. 2, 6), ina 
uniform conductive half-space (Ref. 4, 13, 20, 21), or in a non-uniform, strat- 
ified medium (Ref. 7-17). 


It was found long ago that electromagnetic pulsed fields may be used to 
develop new methods of measuring the flow rates of an electroconductive liquid 
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(Ref. 24-26). In these methods, pulsed currents are induced in a moving, con- 
ductive liquid by switching on the current in the primary coil. The flow /136 
rate of this liquid may be determined by changes in the time intervals of the 
maximum shift (Ref. 24) or the electromotive force passing through zero (Ref. 
25) in the receiving coil. This is caused by the fact that the conductive 
liquid removes the pulsed currents. Research on the propagation of electro- 
magnetic pulsed fields in moving, conductive media was carried out due to the 
necessity of making an exact determination of the dependence between the 
measured intervals and the velocity of the liquid (Ref. 27-29). 


The review by V. V. Novikov (Ref. 22) presents an extensive bibliography 
and a survey of research on the propagation of pulsed signals in conductive 
media and above the earth's surface when there is no motion. In this article, 
we shall confine ourselves to examining the methods and results given in these 
articles, which may be directly applied for studying the propagation of pulsed 
fields in moving, electroconductive media. We shall not deal with the problem 
of considering the displacement currents, since they do not play a significant 
role in this formulation. Thus, the phenomena which we shall investigate are 
limited to problems which include different solutions of the Maxwell equations 
for moving media (Ref. 27) 


ou 
rot E=—po7;; 
Ho Of (1) 
rot H=oE+op,[vH], div H=0. 
We previously disregarded the magnetic properties of the medium (Up = 4r 1077 
h/m). Introducing the electromagnetic potentials 
1 OA 
= j rot AandE=grad p= TE’ 
we obtain JA 
2A — = pels 
VW 7A— opo (VW A) =opo OE (2) 
and 


l 
P= 5), div A—VA. 


2. Propagation of Pulsed Electromagnetic Fields in Conductive Media 

The investigation of the propagation of an electromagnetic field impulse 
in conductive media is considerably simplified when the medium is not in motion. 
In a few of the simplest cases (the medium moves like one entire body, there is 
only one component of the vector potential) there is no necessity of resorting 
to the solution of equations for moving media (2). We may thus confine /137 
ourselves to employing the Galilean transformation (Ref. 28, 29). In this 
connection, it is pertinent to investigate the most characteristic methods and 
results of solving the problem when there is no motion -- i.e., the group of 
problems which include different solutions of the equations 
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(3) 


with the corresponding boundary conditions. 


Let us begin by examining the simplest case, when an electric dipole having 
the length dS is placed in a medium with uniform conductivity o. The current 
passing through it has the form of a single inclusion function 


where S(t) =S3y(t), 
0, <0; 
vo- 1, t>0 tag 


is the unit step function of time. 

The solution of this problem has been obtained by Ollendorf (Ref. 2, 3) by 
means of an operational method, òn the basis of the well known relationship of 
the Hertz vector (Ref. 2) for a dipole supplied by an alternating current 

Sds enor (5) 


n(o) =- F nor 


It is known from the theory of the Laplace transformation (Ref. 30) that, by 
applying the operator L~+ to the constant 9, we obtain 


atio 


j t 
L719= 51 fs © dp=S (0). (6) 
a—{00 


Consequently, substituting p= iw in (5) and applying the operator L to the 
expression (5), with allowance for the well known relationship (Ref. 30), we 
obtain the expression for the Hertz potential 


no-i |!- (Fa) |v. 
where 
2 x 
z 4t? f u? 
t= , O(x) =—= | edu. 
opo? (=) Ve, 
For a dipole directed along the Z axis, the relationship between the /138 


single component of the vector potential AL and the Hertz potential has the 
following form 
A,=opoll, (8) 


and the electric field components may be expressed by means of the Hertz poten- 
tial as 


0p 0” oll 
£,=——t = -a = 
Ey dy’ P Ot , (9) 


PI oil 


E.- Gat oy (9) 


Potential II satisfies the equation 


oll 
Vl =opo 5; - (10) 


For an infinitely long conductor directed along the Z axis (Ref. 2), ọ = 0 and 


nab ffeil a 


“Tno r 
where r=Yp?+¢2 3; 9 is the distance from the observation point up to the axis; 


5S. is the distance from the observation point to the dipole. In this case, 
only one component of the electric field remains 


l oo fð 1 ds 
nolan faoa 


za (12) 
_ 3 1} 
— ngo pE 
; Ē np? = 
Figure 1 presents the dependence of £,= 3 E ont. 


A. N. Tikhonov (Ref. 4) first obtained the most complete solution of the 
problem for the same dipole placed on the surface of a uniform, conductive half- 
space. A. N. Tikhonov and O. A. Skugarevskayaalso solved the problem for the 
dipole placed in a non-uniform, stratified medium (Ref. 7-12, 15-17). The 
solution of all these problems is based on the solution of the thermal con- 
ductivity equation for the second and third boundary value problems. The 
initial equations with partial separation of the variables are reduced to the 
one-dimensional case. The use of the two methods -- the reflection and /139 
the Fourier method -- yielded the dependence of the electromagnetic field over 
a wide range of changes in time t. 


We shall try to present the main results of these studies. Let us first 
examine the case when an electric dipole, which is placed at the origin, is 
directed along the Y axis and is located at the boundary z = 0 of a uniform, 
conductivity half-space. Just as previously, the current changes accord- 
ing to the unit step function (4), which leads to zero initial conditions. This 
orientation of the dipole makes it possible to set A, = 0, and the boundary 


conditions for the remaining components of the vector potential and the scalar 
potential may be reduced to the following equations of continuity: 


Ay le<o= tyll = Ay [no = Act lz=0 ’ 


(13) 
dAy| Ayu 5 
gal n a alaare 
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Figure 1 


Curve Showing the Establishment of an Electric Field of an 
Infinitely Long Conductor in a Conductive Medium 


where Al» 4 are the potentials for the region z>0 and o # 0, and Ary dry are 


the potentials for z<0 and o = 0. In addition to the boundary conditions (13), 
the vector potential must satisfy the condition of being limited at infinity, 


and must have a singularity of order W3 at the origin (Ref. 31). 
We shall try to find the solution ane problem which has been /140 
posited by the method of separation of variables. We shall assume the following 
Aes ek To(r) Y (z, à, t) dd (14) 
and 
A,= nod dy , ip | ONZen Nd (15) 


where r = Ven. It may be readily seen that the functions Y and Z 
satisfy the equation 


PU sry oy, 2U 
Oz? Ot ` (16) 


In view of the well known relationship (Ref. 31) 


j Jo(Ar)e*|2 |= 7 


it is clear that, in order that the vector potential have the requisite singu- 
larity at the origin, the following condition must be fulfilled 
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Yilzz0™ Ynl:= St- Oz 7T (17) 


z=0 


In view of the fact that Yir satisfies the following equation in the case of 


z<0 
Yy 


Jz? =0, 


whose solution is given by the following function 
Y (z, à, t) =u (à, t) e, 
we find that Yn 


Dz lem l 


lz=0 z=0 


Consequently, in order that the first and second conditions of (17) be satisfied, 
it is necessary that 


Yil=0= (à, £), (18) 


oY; 


Dee N 


jz=0 


rae (19) 


Thus, this problem may be reduced to the one-dimensional case of the /141 
third boundary value problem for an equation of the parabolic type. As is 
known (Ref. 32), the solution of the latter problem has the following form 


Seek P ag olto) | 
2 [|e (9H | @ ED ds! 
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With allowance for (14), the following expression is obtained for By 


ate | (ar) 7 iG eo = 
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The boundary conditions for the component A, are obtained from div Arr = 0 and 


(13) or 
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OA, ZS Apu 
oz |s=0 Oy 'z=0 (22) 
and 
Ay \ao= Ay, l= . 


With allowance for (15), expressions (22) lead to the following equation 


2 ðZ 


and the following boundary conditions 


Ls 


Z,| 


jz==0 


=~ jut, N=-4%, 


In order to find A,» following the procedure given by A. N. Tikhonov (Ref. 4), 


let us determine the magnitude of the scalar potential as follows: 


/142 
Sl ana 34.9 fF y OZ 
= a “Tae | Pan Y+ dì = 
0 
_Sdy (25) 
=7 af rons (z, à, t) dh, 
where 
ðZ 
S (z, à, t) =Y (2, A, t) + 5 A, t). (26) 


It is apparent that S satisfies equation (16), the initial zero condition, and 
the following boundary condition 


Shara an =2. (27) 


Flat ri 


Employing the well known solution of the first boundary value problem for an 
equation of the parabolic type (Ref. 31), we obtain 
t 


1 z 2 2 
a ape exp [oO A 
Vropo J ah xP | 4(t—t) om (¢— »] dr. (28) 


Taking into account the following relationship (Ref. 35) 


Jrone ane Tae Toe ES, be) 


z =r 8(i—1) 
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we find that | i ) 
am | 2?-+— 73 

-34.2 aia E p h (29) 

~ 8no Ey (¢— +)? °18 (E= 


Omitting the very cumbersome calculations (Ref. 4), we shall only present certain 
results which will be of interest to us in the future. 


The electric dipole field on the surface of the half-space (z = 0) is 
determined by the following expression: 


d a 
than Fallon!) -Jr J co 
where £143 
1 2 1 eF - 4t 
F(t =0|—~=]|— T = 
2 fae ya YTS a opy? ’ 


where y is the unit vector directed along the Y axis. 


If the source of the field is not a dipole, but any line L which connects 
the points A and B, then the electric field of the line L is determined by the 
integral 


E= f E,(x,4,0,t; $ n)ds= 
L 


(31) 
-gay (raduni) )- fr Bas}, 
where 
E, (x, y, 0, £; & 7 =~ PE {grad yl Os (+)-$r@} 


is the electric field of the dipole located at a point with the coordinates 
&, n(z = 0) in the s-direction. 


When the conductor is located on the Y axis and is infinitely long, a 
simple formula may be obtained from (31) for a single component of the electric 
field 


E,=— Ea (=e), (32) 


and the scalar potential is ọ = 0. We shall assume that the electric dipole 
directed along the Y axis is located at the origin on the surface of a con- 
ductive layer having the thickness 2. Assuming that A,. = 0, we find that the 


continuity conditions must be satisfied in the case of z = 2, in addition to 
the boundary conditions (13), i.e., 


(33). 


Aule = YII l= An low = Fatt PE 
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— Ayn 
z=! dz 


Log Fu leat = Put leat? (33) 


where AGT? oar are the potentials for 0 $ z < 2 and A are the potentials 


III’ rrr 
for z >. 


Thus, as follows from (19), the influence of the region z<0 may be replaced 
by the boundary condition 


ðY 


= - ay! =-—2i. 3 
Oz \2=20 |z=0 a (32) 
Taking (33) into account, we readily find that the influence of the 144 
region z > & may be replaced by the condition 
oY 
ae =0. 35 
oz z=! z=" ( ) 
In a similar way, we have 
1 (36) 
Z =—~Y 
and 2=0 à 1z=0 
z aeiy (37) 
jz=1 À iz=l 
It is clear from conditions (34) - (37) that the solution for a stratified 


medium may be reduced to the solution of the third boundary value problem of a 
one-dimensional equation of the parabolic type (16) with the boundary conditions 
(34) - (37). The solution of this problem was provided by A. N. Tikhonov and 

O. A. Skugarevskaya (Ref. 8-10), who employed two methods which yielded the 
value of the electric field over a wide range of time changes. We shall not 

go into a detailed examination of the calculations for solving this problem. 

We shall only try to illustrate by means of individual examples the general 
pattern of the methods employed by A. N. Tikhonov and 0.A. Skugarevskaya, which 
provided a comprehensive representation of field propagation in stratified media. 


In order to obtain the solution for the initial stage during which the 
field is established, it is advantageous to employ the reflection method (Ref. 
8). As was shown by A. N. Tikhonov (Ref. 8) (for variable t = _t_), in order 


oy 
to obtain the solution of the equation 2 
PY., ay (38) 
gat Yoho oF 
with and initial zero condition and the boundary conditions (34) - (37), the 


function Y may be represented in the form of a series for small values of t 


Y (z, t) = Y,[Yon(2nl+z, A, T) — Yon41(2(n+1)l—z, a, t)]= (39) 
a=0 
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œ 
= Yo (z,%, 1) + D (Yent (2(0+1)1+z, d, t) — Yone (2 (0+1) i- 
n=0 


(39) 
—2, À, t)J]=Yo(z, Ay 1) an Y,(2 l-z, À, t) + Y2(2 {+z, À, 1) ~e. 
All the functions included in the right part of equality (39) are deter- 145 
mined so that they satisfy the boundary conditions (34) - (37) as follows: 
— By HA Yo=2A C0) 
and 
Yang oY. 
BP tA Yon = aot +A Yon for 2=0 (41) 
and 
OY: 
oe +A Yon = OF Yan for z=l. (42) 
All the functions Ya may be completely determined with these conditions. 


The function Yo may be found from the solution for a semi-limited line. The 


remaining functions may be found by means of the operator 


Lilg)= | eg (t) de (43) 
so that 
Yn (z, à, 1) =20(—14+24L,)"(1—A La) vo. (44) 


Vo is determined by means of a fundamental solution of the thermal conductivity 


equation, i.e., 
T 
vo=2 [Giz t)di= 
0 


: (45) 
1 l 2 
e ie | -ayeo a. 


As was shown by A. N. Tikhonov (Ref. 8), the convergence of the series (39) is 
determined by the condition 


| Yan (2nl-+2, t) — Yonn (2(n-+1)1—z, t) | <9" ce- E. (46) 


As may be seen directly from (46), the series converges particularly rapidly 
for small t. For the initial stage, it is absolutely permissible to confine 
oneself to the first three terms of series (39), whose values expressed by 
means of the operator L, are as follows: 


À 
Yo=2Av9—2A2L (vo); 
Y, = ~—2A09+ 6A7L (vp) — 4XL? (vo); (47) 
Y2=2Avo— 10ML (v0) +1643L? (u0) — 84L? (vo). 
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Following the procedure of A. N. Tikhonov, in order to determine the /146 
function S (26) let us introduce the auxiliary function W (z, à, T), so that 


Z(z, mt) =—2V(2, dt) + W (2, Ay). (48) 
According to conditions (36) (37), we find that W (z, à, T) satisfies the same 
equation as Y, with initial zero conditions and with the following conditions 

W (0,2,1) =0; W(l4,7) =Y (L 4,7). (49) 


We are interested in the magnitude of the electric field on the surface of a 
layer, and consequently the quantity we desire is 


A 1 ðY 2 aW 
E E 
; (50) 
g A Oz Ig 


According to condition (27), the first component represents the value of S for 
2 =o (in the case of half-space), and the second represents the correction for 
the finiteness of the conductive layer thickness. Thus, the computation of S 


may be reduced to calculating the function W. Representing it in the form 
of series æ 
W (z, A, T) = $ [Won (z, À, t) = Wont (z, À, t)], (51) 
n=0 


according to conditions (48) we find that 
W,(z,4,0)=0, Wn (0, A, 1t) =0, 
Wan (1, A, t) = Yon [(2n +1), A, x], 


Won (4, A, t) = Vonsil (22+ 1), À, t). 


Consequently, the functions Won and Wontl are solutions of equation (38) 
with the conditions (34) - (37). It is apparent that they equal 
oo 
Won (z, bt) = }, (Yon [2(n +2+1)l—z, a, 1] — (52) 
k= 
— Ym [2 (n+k +1) +2, A, t)} 
and /147 
o0 


Wanti (2, à, 1) =}, (Yonu [ntk+1)l—z, A, sJ— 
k=0 
— Yann [2(n +2 +1)l+z, A, tH}. 
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As the first approximation, we may take 
W (z, à, t) = ¥o(2/—z, À, t) —Yo(2l+z, A, t) — 


(54) 
—Y,(2l—z, A, t) + ¥1(2l+2z, A, 1) 
and 
OW) __ {0% _ oY, 
Jz lT? a a] > a 
For the desired function S, we obtain 
= OY, - OY, 
S|,.9=2— ae z e: (56) 


Thus, the approximate value of S in the case of z = 0 may be expressed by 

means of Y, and Y,. 
0 1 

The study (Ref. 9) presents detailed computations of the functions Yoo 
which are employed to express the dependence of the electric field in which we 
are interested. In view of the cumbersome nature of the final formulas, we 
shall not present the end results here, which have no special interest for our 
purposes. We shall next investigate the solution of the problem for the final 
stage during which the field is established, which was given in (Ref. 10). 


In order to find the function Y (z, à, t) of our preceding problem in a 
form which is suitable for the investigation in case of large t, it may be 
conveniently represented in the following form 


Y (z, A, t) =Y (z, A) +F (z, a, t), (57) 
where y® is the stationary solution, and Y is the deviation from the stationary 
solution. In order to determine these functions, according to (34) - (38) we 
have the following condition 

ZYD BY =O; 0<2<l; 
Ot ayos —2, for z=0; (58) 

(0) 

oe +aAYO=0 for z=l 
and ay F 148 
av —i1Y=0 for 2=0; 
pa (59) 
Jz -s—+AY=0 for z=% 


F(z, à, 0) =— YO (z, a). 
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It is apparent that 

YO (z, 4) =e-%, 
Assuming 

Y(z, A, t) = Yn (Z, A) Tn (A t), 
we find that 


ar +hn2¥n=0, (60) 
ots =0 for 2=0; ce) 
OFa Yq = =0 for z=! 
and 
Ta=Ce-",t, 
where 


t 
= 2 2 = 
pn Rn +A and T ou” 


The general solution of equation (60) is 
Yn=C, cos knz + Cosin knz. 
We find the following from the boundary conditions (61) 


2knl 
C= Ec and tg kn l= pa ae 7 yr 


Separating the right and left part of the latter equality for n = 0 in powers 
of A and retaining terms up to the third order, we obtain 


1. 2 
2 -2 -HÉ IR. (62) 
The function Y may be ee in the following form /149 
Y(z, 4,4) = os CanYn (Z, A) e-n tA, (63) 


The coefficients C_ are determined as Fourier coefficients of the initial func- 


6), 


tion 


JV (2,4) Ya (2, 4) dz 


C,=— 2 (64) 


Ca 


i [Yn (2, à)]?dz 
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(0) 


Taking into account the values of Y 
l 


f Y0 (2,4) Ya (z, A)dz= 2h 
0 


and Ya’ we readily find that 


(65) 
Rn? FAR" 


With allowance for the expansion of (62) we obtain 
t 
fire, A)Fdzwt+ 1 yyy 2 e (66) 
ò 3 45 g 


Consequently, 


an 2a 
DT pes ea ae Oe (67) 
(ko? +42) (1+ Fate 5 PA?) 


and we obtain the asymptotic formula 


24 [cos koz +e sin Roz] 
Y(z, 4, t) =— 


EARS the (68) 


(k? +9?) var! = 2) 


of 
which may be applied to the case/small values of à, which we shall be interested 
in later on. In order to illustrate the use of the formula obtained, let us 


calculate the derivative aA, in the case of z = 0: /150 
at 
dA, sE 
DE ha froni Y (0, A, t)dìi = 
(69) 
ueSidy 24 exp |- G A+ 2 4 zo) | 
= He f non dà. 
l+ +a te =P) 
The presence of the exponential factor exp |- 2 At] leads to the fact that the 


value of the integral for large t is determined by the value of the integrand 
for small à, which substantiates, the use of asymptotic formula (68). Decompos- 
ing the factor J Qr) exp - = àt] under the integral sign in series with 


& . 
respect to à and confining ourselves to three terms of the expansion we obtain 


OA, wSdy T 2, 2., 2/1, 2 _2 
2A, „MoS Y fhiar) [m -32+3 (5!-77)* Ie T^ d) = (70) 
0 
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woSdy [2 20 =f A _ I2% 
Ua, 3| in, en aC 
cee Ga a, 


: 70 
+3(5-7") 18x 12073 (70) 


= 
dt 


since 
Jo(Ar) edd = salio 
Y+? 
and 
JolAr)are-a - ais 
[i )are-ada = (—1)9 2. ae 
The function W may be determined as the solution of the boundary value /151 
problem 
Cw ow 
age VW = one Gps 
W|,,=0, WI_,=0, 
z=0 t=0 (11) 
œo 
War =Y lme > CaYa(l, A)e-#, 5, 
a=0 
which leads to the asymptotic formula (Ref. 10) 
(72) 


W (z, A, t) ~e-™ shiz Mg +Co(A) Yoan È sin na pent 


We shall not continue the calculation for the electric field components 
on the surface of a conductive layer, and we shall thus conclude the examination 
of the methods developed by A. N. Tikhonov and 0. A. Skugarevskaya. These 
methods were employed to obtain the relationship for the initial and final stage 
during which the field was established in a conductive layer lying on an ideally 
conducting base (Ref. 11-12), as well as to determine the process during which 
the field is established in a three-layered medium (Ref. 15). Figure 2 presents 
curves showing the dependence of E= zor? p on t = 4t taken from (Ref. 

Sdy ~¥ ouor? 

7), which were computed according to formulas for the electric field of a dipole 
lying on the surface of a conductive layer. Calculations according to formulas 
for the initial stage were performed for t < 0.5, and for the final stage -~ for 
t 20.1. The dashed curve presents the results derived from calculations for 
half-space. As may be seen from the curves, the presence of a limited layer 
leads to an increase in the field amplitude. This may be explained by reflection 
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of the field from the layer surface. 


We may obtain the electric field o 
on the Y axis, which is given 
ductive medium and half-space 
the dipole directed along the 
expressions, it is not always 
and then to integrate over an 


by integr 
Y axis. 

advantage 
infinitel 


f an infinitely long conductor, located 


by the simplest formulas, in the case of a con- 


ation of the expressions obtained for 
In view of the complexity of the dipole 
ous to solve the problem for the dipole, 
y long conductor. Since we have only 


one component of the vector potential Ay in the last case and since all quantities 


depend on the two coordinates x and y, we may regard the problem as a two-dim- 


ensional problem. 


A, of the vector potential vanish. The 


mined by expression (14), may be determ 
may be readily seen: 
ates fr 
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Figure 2 


Curves Showing the Establishment of the 
Electric Dipole Field on a Conductive 

Layer Surface for Different 2 = P 
where 2 is the Layer Thickness; r - is 
the Distance from the Dipole Along the 


X Axis. 
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In view of this fact, 


the scalar potential and the component 
component 4, which may be deter- /152 


ined by the following transformation, as 


Z, A, t)cos Ax x: 


(73) 


The problem of determining Ay may be 


reduced to applying the methods which 
we investigated above for determining 
the function Y. 


The methods developed by A. N. 
Tikhonov and 0. A. Skugarevskaya, which 
are based on an analysis of the solution 
of the thermal conductivity equation 
and which are rather complex in nature, 
make it possible to analyze the solution 
graphically over a wide range of changes 
in time. In order to examine similar 
problems, other authors have employed 
an operational method (Ref. l-2, 5, 6, 
13, 14, 18-21, 27-29). As we shall sub- 
sequently see, the operational method 
may be visualized more clearly, as com- 
pared with the methods of A. N. Tikhonov. 
Difficulty is entailed here when we de- 
termine the inverse transform and ana- 
lyze the solution. D. N. Chetayev (Ref. 
13) has provided a graphic application 
of the Laplace transformation for solv- 
ing the problem regarding the /153 
establishment of the field in stratified 
media, when the conductivity of the two 
layers differs only slightly. 


Let us now investigate the propaga- 
tion of electromagnetic pulsed fields 
in moving media, where the operational 


method will be illustrated at the same time. 


3. Propagation of Pulsed Electromagnetic Fields in Moving 
Electroconductive Media 


Let us begin the examination of the propagation of pulsed electromagnetit¢ 
fields in moving media with the problem in which the source of the field is an 
electric dipole located at the origin. The vector potential of this dipole 
satisfies the following equation (Ref. 29) 


OA A 
y’A—opo(vy)A— opo z = — Hobo (74) 
and the initial condition 
Aļ:=0=0, (75) 


where Jz = QSdrd(r) and (r) = 6(x)6(y)6(z) is the three-dimensional delta-func- 
tion. 


In order to simplify subsequent computations, we shall assume that the 


medium moves in the direction of the X axis. In this case, it is natural to 
assume that vV 7V” Yy he a 0, A = A, = 0 and iy satisfy the following 


equation 


ðA 0A 
VAy—opov 3” — opoz] = — woSdy6(r). (76) 


Applying the Laplace transformation 


[- 2] 
A= | Ayer'dt=LA, (77) 
6 
to equation (76), with allowance for the initial conditions (75), we obtain 
0A, 
AA — opo g TPA, =— 5 mo3dyð (r). (78) 


The solution of equation (78) may be veitt ene an3 of the Green's function for 
a point source (Ref. 34): 


A, = boSdy Í dk, J d z 79 
A= np ho J ketikan tomp” ae 
where kr = k ee tke z and ‘Pe kx 4k, 417 . The inverse transform of 154 


the integrand may be Sonenn B. by applying the multiplication 
theorem for the Laplace transformation Eer: 30): 


Lai 
aait fa fafaj aneo 


Xcos kyy cos k:zdr. 


(80) 
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Integrating over k? E and kK» we obtain the following expression 


A 8H [beso fg 4r qe (E= o) tty] | de, GY 
o 


which, as may be readily seen, may be reduced to the expression for the Hertz 
potential (7) in the case of v= 0, with allowance for (8). 


The same result is obtained for the case when the dipole is directed along 
Z, and the motion takes place in the X direction, just as previously. For a 
dipole directed along X (i.e., in the direction of motion of the medium) we have 


Ant 


A -2t fire Fexp ae; =42 ((x- —or)* +2424 |e (22) 
o 


The scalar potential may be determined by the following relationship 
; 
H div A—v, Ax. (83) 


In order to find the field of an arbitrary circuit or conductor, we must inte- 
grate the dipole potential over this circuit or conductor, assuming beforehand 
that it is formed of elementary dipoles, whose vector potential is determined 

by formulas (81) and (82). 


By way of an example, let us calculate the electric field of a circular 
circuit in the plane (Y, Z) with the center at the origin. The vector potential 
of the dipole located in the plane (Y, Z), with the coordinates ,g and the 
length ds, has the Hon DE form: 


4-38 se jig Fexp| — FE eot ya) -0h ar. ee 


Introducing the cylindrical coordinates /155 


y=rcos), z=rsiny, 
N=acosa, C=asina, 
ds=ada, 


we obtain the expression for the single component Ay of the vector potential 


of the circular circuit: 


Sa flop o 
hen ( Si exp — SP a24 r24 (xo a) 
4 7) "5X 
o J 4t (85) 


x feof ar cos(a— Y) } cos adadr. 


-x 
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In the case of cylindrical symmetry, we may set y = 0 and, employing the integral 
representation for the Bessel function of a purely imaginary arguement (Ref. 31), 
we obtain 


t 
2 flori 
E j af- sto [a? + 2+. (ory h(a ae (86) 
0 


The scalar potential ¢ = 0, and the following expression is obtained for the 
electric field 


OAq - 2na3 (2¥o\% opo 

sce e SS sf — ‘entice ee. 2 eee Cg ar 

Sana ag (a) exp Po fat 4 re- (x on) (ae ). (87) 
The electromotive force, directed in the circular circuit of radius r which 


is located at the distance x from the first circuit, may be expressed by the 
integral 


2x 
E= f Eyrdg=2arE, . 


(88) 
Consequently, we have 
= 1 a SSS = 
= -(z) e{- þu +r?+(x-v Dan (F) (89) 
where 
g oa? 34t. -_ OS, .>7_1 
€ NS Ç; b= a r= a’ x= gard v= 4 Ouoav 


are dimensionless variables. The latter results in the case of r = a coincide 
with those obtained in (Ref. 27). 


Let us present still another method for solving the problem for mov- /156 
ing media when a solution of equation (2) is not required, and when we may 
confine ourselves to employing the Galilean transformation (Ref. 28, 29). This 
method is sometimes simpler, as compared with the solution of the problem with 
an equation for moving media (Ref. 27). In order to examine this method, let 
us introduce the four-dimensional potential % with the components (Ref. 33) 


DAs D=Ay OAs W=Ly, (90) 


and the four-dimensional radius vector x with the components 
Xi=X, XQ=Yy, xy=z, xq=ict. (91) 


When changing from an inmobile system of reference to a system which moves in 
the x direction with the velocity v with respect to the first, according to the 
theory of relativity, the four-dimensional potential and the coordinates may 
satisfy the Lorentz transformation 


D=, D=O; (92) 
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D= — 
1- (92) 
a 
and 
xX +i— x, 
x= a x = X2, x’ =X. 
y! wv 2 2 3 3» 
-F 
(93) 
x4—i— x, 
x= F 
y'a 
Examining the case of slowly moving media and disregarding terms containing 
&/c2, let us turn to the Galilean transformation 
A’n=Ax, AY=Ay AL=A, (94) 
g =9-0A, 
and {157 
x'=x—vt, y'=y; z =z; ľ=t. (95) 


Employing the transformation formulas obtained for electromagnetic potentials, 
we may readily obtain the solution for moving media. For this purpose, we must 
first solve the problem in a coordinate system moving together with the con- 
ductive medium, by transforming the potential and the coordinates of the im- 
mobile source of the field according to (94) (95). Employing the inverse trans- 
formation, we may then change to an immobile coordinate system. 


Let us illustrate the statements presented above with the example of a 
two-dimensional problem (Ref. 28). Two infinitely long conductors, with 
different currents moving in the opposite direction and having the magnitude 
3 , are arranged on the surface of moving half-space at the height h in the 
Y direction. The conductors are located symmetrically with respect to the YOZ 
plane. The XOY plane coincides with the surface of the half-space. The currents 
change according to a unit step function of time (4). In this case, the scalar 
potential vanishes, and only one component of the vector potential E = A re- 


mains, which satisfies the following equation in a system which is rigidly 
connected with the moving half-space: 


PÄ PA 0A’ 
Ox? T zit TO OFF oe 
with the initial condition 


A’(x’, 2’, 0) =0. (97) 
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The continuity conditions at the boundaries z' = 0 and z' = h are the boundary 

conditions, just as previously. The derivative 9A' on the surface z' = h 
az' 

undergoes a discontinuity, which equals the density of the surface current j': 


A'i (x’, h, t) =A’) (x, h, t’) , 


0A’) 0A’; | =p ’ 
D 5 -izr a = Hoj » 
Oz' ‘s'=n OZ zsh (8) 
'n (x’, 0, t) =A’ m (x’, 0, t) , 
duj _2A'm 
OZ’ ‘s= OZ |er=0" 
where A'i is the solution of equation (96) for the region z' >h, Ve peas 


solution for 0 « z' x h, and A’ is the solution for z' < 0. In addition to 


IIL 
the boundarv conditions, we must take into account the limiting condition /158 
for z' > +, 


The density of the surface currents with respect to the immobile coordinate 
system may be expressed by means of the delta function 


œ 
j=S[d(x—a) -oxta)}= 3 f sin Ax sin adh. 
0 


In the mobile coordinate system, we obtain the following by means of the trans- 
formation (95) 


co 
Pata’ . 
273 sin A(x’-++ ut) sin adh. 


0 


Applying the transformation (77) to equation (96) and boundary conditions (98), 
we obtain the equation 


GA’ PÄ ; 
gt + gga = OHA (99) 
and the boundary conditions 
A(x’, h, p) =A’ n(x’, h, p), 
0A'y 0A’, ae 
Oz’ laran OZ’ laan CO? (100) 


A'n(x’, 0, p) =A’ m(x, 0, p), 


04‘) Am 


OZ’ ‘gn  OzZ' 


’ 
2’=h 
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where 
23 f 
f= f F sin dada 
0 
and 


F=? sin Ax’ + Av cos Ax’ 
pits 


Separating the variables of the last equation, we may readily obtain the solu- 


tion satisfying the boundary conditions (100) and the limiting condition for 
z' + +o, 


After all the constants are determined, the solutions for all three {159 
regions assume the following form: 


oo š 
At A Ar 2 X r dh 
A'=A n=A =S f Fe-z'+h) TE a 
0 


x (101) 
da woS f dd 
t ra P 
13 “fy eren ES f Fene -MS, 
0 0 
Ale 2m j Fe-hhrnr i (102) 


where x “Vi! + SUgp i 


The inverse transform of the integrand may be found according to the well 
known formulas (Ref. 35) by employing the multiplication theorem: 


t © 
A’= us far f eM +M (A, t) sin A(x’ + of — vt) sin Kadi — 
0 


(103) 
wi a 
7 He f e-Mz'+h) sin A(x’ + vt) sin KST 
0 
5 co) 
+ He J e-Mz'-M sin A(x’ + uf) sin Aa a, 
u 
(104 


œo 


A= au f e-MN (2’, A, t) sin A(x’ + vf—vt) sin adh, 


where 
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and Opoz’? AT 
o(a- 


Vxopot 
ae er ee opoz”? ga 
žar [i o(y/ re +} II. 


Changing to the stationary coordinate system,according to (94), (95) we 
obtain 


N (Z’, a 1) = 


t co 


A= 2w3 a f e-z+h) M (à, x) sin À (x— ut) sin Aada— 
ò o 


—y(t) A +y (As 
œo 


t 
An= 2S f ds f e-AN (z, à, t) sin.A(x— ot) sin Aad, 
o 0 
where 


o0 
A= B e-^z+h) sin Ax sin Aa a 


g a 


is the expansion of the vector potential of the reflected currents and 


3 Ff dd 
w3 f e™^z-h) sin Ax sin Aa a 2h; 
0 


A=} | 


3 F dh 
BS f eMe-M) sin Ax sinha Z, z<h 
o 


is the expansion of the vector potential of the external currents. 
The strength of the electric field for z 20 is 


A 
Ey= — pf Et Ep = Ef”, 


where 


o 
E= LS f e-MetMM (a, £)sind(x—vl)sin hadh, 
o 


Ey =0(t)A s; Ey 8(t)Ay™ 


and ô (t) is the impulse function. 


(106) 


(107) 


(108) 


(109) 


(110) 


We shall omit the very cumbersome calculations of the integral (110) ,-‘and 
we shall only present the final result for the field of the induced currents: 


143 


w (et) -1 w (=) 1 w (42) -1 
yi yi 


Fx yt 
=- ut (ge ere OT 
w (=225) a ie 
ac ee oe ee 
(Ee + if)? 7a (aree Ex? +? }’ /161 
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w (z) =e77 ( +7 f edu) 
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is the probability integral of the complex argument which is a tabulated func- 
tion (Ref. 36); 
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Figure 3 


Dependence of the Function EH on t and Vy in the Case of 


x = 2 and ¢ = 0.5 
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For the special when ¢=0 we have /162 


S BA 1l A 
E; =2|zz(- ŽE) — perp -$ (112) 


We may obtain the latter formula by means of (32), summing up the field of the 
two conductors and employing the substitution x > x - ut. Figure 3 presents 
curves showing the dependence on t and u, computed according to formula (111) 
for x= 2 and 5 = 0.5. 


As the last examples have shown, when allowance is made for the transforma- 
tions (94) (95), we may readily change from computations which are performed 
for a conductive layer to computations for a moving conductive layer. For the 
case of an infinitely long conductor, these results are only achieved by re- 
placing x > x - ut. 


In conclusion, we would like to point out that the transition of the func- 
tion (111) through zero between two main maxima (see Figure 3) is determined 
by the condition x = vt. It may be readily seen that the electromotive force 
induced in a rectangular coil, which is arranged vertically with respect to the 
surface of half-space or the moving layer, satisfies this condition. This fact 
has been employed to compile a new method for measuring the velocity of con- 
ductive media (Ref. 25); this method is based on the transition through zero of 
the voltage in the measuring coil. 
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